import numpy as np
import matplotlib.pyplot as plt
from spectractor.extractor.psf import *
Ny = 50
y = np.arange(Ny)
amplitude = 10
norm = evaluate_moffat1d_normalisation(gamma, alpha)
a = evaluate_moffat1d(y, amplitude=amplitude, y_c=Ny/2, gamma=5, alpha=2, norm=norm)
plt.plot(a)
plt.grid()
plt.xlabel("y")
plt.ylabel("Moffat")
plt.show()