[x y]=meshgrid(-100:0.5:100,-100:0.5:100);
y1=exp(-(x.^2+y.^2)/(2*s1^2));y2=exp(-(x.^2+y.^2)/(2*s2^2)); image((y1-y2)*500);
