a = 5; // demi-largeur du domaine de définition
// en: half width of the domain of the function
h = 1; // hauteur du pic
// en: peak height
p = 0.02; // pente de la ligne de fond
// en: slope of the background line

theta = linspace(-0.5*%pi, 0.5*%pi, 51);
foo= cos(theta); // pour avoir plus de valeurs au sommet qu'aux bords
// en: to have more values close to the summit than at the rim
clean(foo);
x = gsort(a*[1-foo, foo(2:$)-1], "c", "i"); // domaine de définition symétrique entre -a et a
// en: symmetrical domain of the function between -a and a

y = h*(exp(-x.^2) + 0.5*exp(-(x + 1.5).^2) + p*x); // deux pics gaussiens
// en: two gaussian peaks

plot(x, y);