import matplotlib.pyplot as plt
import numpy as np
from spectractor.extractor.psf import DoubleMoffat
p = [2,20,30,4,2,2,5,1.5,20]
psf = DoubleMoffat(p)
yy, xx = np.mgrid[:50, :60]
out = psf.evaluate(pixels=np.array([xx, yy]), values=p)
fig = plt.figure(figsize=(5,5))
plt.imshow(out, origin="lower")
plt.xlabel("X [pixels]")
plt.ylabel("Y [pixels]")
plt.show()