function z=dblgauss(x,y) z=exp(-x.^2-y.^2);