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