function gauss = Gaussian_2d ()

x_start=-3.0;
x_end=3.0;
x_sub=0.1;
Nx=(x_end-x_start)/x_sub;  %number of points along x
y_start=-3.0;
y_end=3.0;
y_sub=0.1;
Ny=(y_end-y_start)/y_sub;  %number of  points along y

for i=[1:Nx+1]
  for j=[1:Ny+1]
       x(i)=x_start+i*x_sub;
	   y(j)=y_start+j*y_sub;
	   f(i,j)=exp(-(power(x(i),2)+power(y(j),2)));
  end
end

gauss=f;
return
end
