Fit 2D Gaussian with Optimization Toolbox

Fit a 2D rotated gaussian. http://en.wikipedia.org/wiki/Gaussian_function

FMGAUSSFIT performs a gaussian fit on 3D data (x,y,z).
[fitresult,..., rr] = fmgaussfit(xx,yy,zz) uses ZZ for the surface
height. XX and YY are vectors or matrices defining the x and y
components of a surface. If XX and YY are vectors, length(XX) = n and
length(YY) = m, where [m,n] = size(Z). In this case, the vertices of the
surface faces are (XX(j), YY(i), ZZ(i,j)) triples. To create XX and YY
matrices for arbitrary domains, use the meshgrid function. FMGAUSSFIT
uses the lsqcurvefit tool, and the OPTIMZATION TOOLBOX. The initial
guess for the gaussian is places at the maxima in the ZZ plane. The fit
is restricted to be in the span of XX and YY.
See:
http://en.wikipedia.org/wiki/Gaussian_function

Examples:
To fit a 2D gaussian:
[fitresult, zfit, fiterr, zerr, resnorm, rr] =
fmgaussfit(xx,yy,zz);
See also SURF, OPTIMSET, LSQCURVEFIT, NLPARCI, NLPREDCI.