This link describes the algorithms used by lsqcurvefit and the other least squares optimizers
BTW, I believe the trust-region method is the default for LSQCURVEFIT, not Levenberg-Marquardt.
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!