gauss2D_draw

 p(1) = peak height
 p(2) = mean x
 p(3) = sigma x
 p(4) = mean y
 p(5) = sigma y
 p(6) = offset gradient x
 p(7) = offset intercept x
 p(8) = offset gradient y
 p(9) = offset intercept y

0001 function z = gauss2D_draw(p,x,y)
0002 
0003 % p(1) = peak height
0004 % p(2) = mean x
0005 % p(3) = sigma x
0006 % p(4) = mean y
0007 % p(5) = sigma y
0008 % p(6) = offset gradient x
0009 % p(7) = offset intercept x
0010 % p(8) = offset gradient y
0011 % p(9) = offset intercept y
0012 
0013 
0014 if length(x) < length(y)
0015    diff = length(y) - length(x);
0016    added = zeros(1,diff);
0017    x = cat(2, x, added);
0018 elseif length(y) < length(x)
0019     diff = length(x) - length(y);
0020     added = zeros(1, diff);
0021     y = cat(2, y, added);
0022 end 
0023 
0024 xmesh = x;
0025 ymesh = y;
0026 
0027 x_slant =p(6)*xmesh;
0028 y_slant =p(7)*ymesh;
0029 
0030 z=p(1)*(exp(-0.5*(xmesh-p(2)).^2./(p(3)^2)-0.5*(ymesh-p(4)).^2./(p(5)^2)))+ x_slant + y_slant + p(8);
0031 
0032 
0033 end
0034 
0035