import numpy as np
import matplotlib.pyplot as plt
from astropy.modeling import models
X = np.arange(100)
PSF = models.Moffat1D()
p = (50, 50, 5, 2)
Y = PSF.evaluate(X, *p)
plt.plot(X, Y)
plt.show()