docs: def eval(e, amplitude, e_0, alpha, beta):
          """One dimenional log parabola model function"""
  
          ee = e / e_0
          eeponent = -alpha - beta * np.log(ee)
          return amplitude * ee ** eeponent