import numpy as np
import matplotlib.pyplot as plt
from spectractor.extractor.psf import *
Ny = 50
y = np.arange(Ny)
amplitude = 10
a = evaluate_gauss1d(y, amplitude=amplitude, y_c=Ny/2, sigma=2)
plt.plot(a)
plt.grid()
plt.xlabel("y")
plt.ylabel("Gauss")
plt.show()