amp=10

for j=1:ny;
    for i=1:nx;
        c(i,j)=amp*exp(-((i-i0)/Rx)^2);
save 'h_init.dat' c -ascii;