[Help-gsl] Fitting data in ODE system/model

[Christian Dudek]

Hello everybody!

I've got an Octave program that fits experimental data into a model of 
different ODEs. The Octave program uses the leasqr() function to fit the data, 
which uses the Levenberg-Marquardt-Algorithm.

My task is to reimplement this in C++, but I haven't found any library yet, 
which provides a similar function. I've seen that GSL provides functions to 
simulate an ODE system without experimental data, or fit data into a single ODE.

The ODE system look like that:

dA/dt = -k1*A
dB/dt = k1*A-k2*B-k3*B
dC/dt = k3*B-k4*C

Has anyone an idea how to solve this problem with GSL? I can't imagine that 
this is possible in Octave (or R), but nobody has accomplished that before in 
C++.

Thanks for your help!

Christian