Hi Mathieu, I think there is also a Levenberg-Marquardt implementation in ITK which may be better suited to perform non-linear least squares fits. I don't know how to do this with vnl_powell, I am currently investigating Levenberg-Marquardt for parameter fitting. Kind regards, Markus