code: def intensity_at_radius(self, radius):
          
          return self.intensity * np.exp(
              -self.sersic_constant * (((radius / self.effective_radius) ** (1. / self.sersic_index)) - 1))
docs: Compute the intensity of the profile at a given radius.
  
          Parameters
          ----------
          radius : float
              The distance from the centre of the profile.