function y=gr(x,r)
y=sqrt(r/pi)*exp(-r*x.^2);