code: def intensities_from_grid_radii(self, grid_radii):
          
          return np.multiply(np.divide(self.intensity, self.sigma * np.sqrt(2.0 * np.pi)),
                             np.exp(-0.5 * np.square(np.divide(grid_radii, self.sigma))))
docs: Calculate the intensity of the Gaussian 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.