code: def intensity_prime(self):
          
          return self.intensity_break * (2.0 ** (-self.gamma / self.alpha)) * np.exp(
              self.sersic_constant * (((2.0 ** (1.0 / self.alpha)) * self.radius_break) / self.effective_radius) ** (
                      1.0 / self.sersic_index))
docs: Overall intensity normalisation in the rescaled Core-Sersic light profiles (electrons per second)