I have a coupled system of differential equations, which I defined in the following way:

dx= zeros(2,1);

dx(1)=K1*K2*x(1)*x(2)./(K3+x(2));

dx(2)= -K1*x(1)*x(2)./(K3+x(2));

I have some experimental data for x(1) at a given time grid. (I have about 30 values).

I would like to find the values of the parameters K1,K2,K3 such that the solution of the ODE fits these experimental data.

How can i program this in matlab?? I'd like to use a least square minimisation to get my parameters… how can I use the function ''lsqcurvefit'' in the optimisation toolbox starting from a coupled system of differential equations??

thank you very much in advance.

simona

## Best Answer