I've looked at the examples given in PDL-levmar documentation. Question: if I need to do curve fitting on function with 2 coordinates, how do I do this? For example: x = p0 * exp(-t1*t1 * p1) + p2 * exp(-t2*t2 * p3); So the above function has 2 coordinates: t1 and t2. Thanks! Bobt