function y=glag_gauss(x)

x=x/sqrt(256);  % sqrt(h)*x where h is the spatial stepsize
y=exp(-(x.^2));