function y=Gauss_2D(x1,x2)

y=exp(-x1.^2-x2.^2);
