import numpy as np
import matplotlib.pyplot as plt
from spectractor.extractor.psf import *
Ny = 50
y = np.arange(Ny)
amplitude = 10
gamma1 = 5
alpha1 = 2
eta = 2
gamma2 = 4
alpha2 = 3
norm1 = evaluate_moffat1d_normalisation(gamma1, alpha1)
norm2 = evaluate_moffat1d_normalisation(gamma2, alpha2)
a = evaluate_doublemoffat1d(y, amplitude=amplitude, y_c=Ny/2, gamma1=gamma1, alpha1=alpha1, eta=eta, gamma2=gamma2, alpha2=alpha2, norm1=norm1, norm2=norm2)
plt.plot(a)
plt.grid()
plt.xlabel("y")
plt.ylabel("DoubleMoffat")
plt.show()