function  y=normal(a,x);
% returns values of normalized Gaussian function with parameters a
 y=1/(sqrt(2*pi)*a(2))*exp(-(x-a(1)).^2./(2*a(2).^2));
return
