function f = bell(x,a)

f = exp(-0.5*a*x.^2)/sqrt(2*pi);

plot(x,f)

