function y=glag_PeronaMalik(x)

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