Hi all,

I am trying to use Levenberg-Marquardt method to fit an incomplete gamma function, but the routines included from NR (mrqmin etc.) donot work at all, i.e., cannot converge. I have tried to let the code write "da", it was so large.

Not sure where I can find a workable and stable source code in Fortran 77 to do this fit?

Thanks.

Yun