code: def flux(self, emin, emax, params=None):
          
          params = self.params if params is None else params
          return np.squeeze(self.eval_flux(emin, emax, params, self.scale,
                                           self.extra_params))
docs: Evaluate the integral flux.