code: def create_from_flux(cls, params, emin, emax, flux, scale=1.0):
          
          params = params.copy()
          params[0] = 1.0
          params[0] = flux / cls.eval_flux(emin, emax, params, scale=scale)
          return cls(params, scale)
docs: Create a spectral function instance given its flux.