code: def intensities_from_grid_radii(self, grid_radii):
          
          np.seterr(all='ignore')
          return np.multiply(self.intensity, np.exp(
              np.multiply(-self.sersic_constant,
                          np.add(np.power(np.divide(grid_radii, self.effective_radius), 1. / self.sersic_index), -1))))
docs: Calculate the intensity of the Sersic light profile on a grid of radial coordinates.
  
          Parameters
          ----------
          grid_radii : float
              The radial distance from the centre of the profile. for each coordinate on the grid.