% g3000    2-D Gaussian
%
%       function A = g2(sz,s)
%
% sz : size of returned matrix, [s1 s2]: variance
 
%       M Jagersand 2000
 
 
function A = g3000(sz,s)
  A = exp(-(((1:sz(1))' -(sz(1)+1)/2).^2/(s(1)*2))*ones(1,sz(2))...
          -ones(sz(1),1)*(((1:sz(2))-(sz(2)+1)/2)).^2/(s(2)*2));


