function y=gauss(x)
% function y=gauss(x)
% gauss(x) = exp(-(x/epsilon)^2)/(sqrt(pi)*epsilon)
epsilon = .01;
a=1/(sqrt(pi)*epsilon); 
b=x/epsilon;
c=b.*b; 
y = a*exp(-c);


