clear all;close all
amp=5;x0=80;Rx=20;dx=1;y0=32;Ry=15;dy=1
nx=160;ny=64;
for j=1:ny;
    for i=1:nx;
        c(i,j)=amp*exp(-((i*dx-x0)/Rx)^2)*exp(-((j*dy-y0)/Ry)^2);
    end
end
figure(2);
surf(c');
save 'h_init.dat' c -ascii;