Thanks for spotting this ... actually there was a missing minus sign in `eval`: `np.exp(x / x_cutoff)` should have read `np.exp(-x / x_cutoff)`.