[SciPy-user] Parameter estimation / Data fitting in scipy

[Nils Wagner]

Hi,

Suppose it is desired to fit a set of data y_i to a known model 
[given in form of an IVP (initial value problem)]

y' = f(y,p,t) , y(0) = y_0(p),       y' = dy/dt

where p is a vector of parameters for the model that need to be found.
y denotes the state-variables. The initial conditions y(0) may also
depend on 
parameters.

How can I solve this problem using scipy's optimization and ode tools ?

A small example would be appreciated.

Thanks in advance

                      Nils