function gg=ggauss(x,y,x0,y0,a,b,amp)
dx = (x-x0)/a;
dy = (y-y0)/b;
gg = amp*exp(-(dx.*dx+dy.*dy));