docs: def eval(e, amplitude, e_0, alpha, e_cutoff, beta):
          """One dimensional power law with an exponential cutoff model function
          """
  
          xx = e / e_0
          return amplitude * xx ** (-alpha) * np.exp(-(e / e_cutoff) ** beta)