[SciPy-user] odeint with digital data

Ryan Krauss ryanlists at gmail.com
Mon Jun 4 11:24:30 EDT 2007


I have a vector of experimental data that I need to use as part of a
system of ode's.  I would like to solve this system using
integrate.odeint.  Can odeint be forced to solve only at the discrete
points in time where the experimentally measured signal is available?
Or do I need to set up some interpolation function to find the
measured signal at any time?  If that signal is the only thing that
explictly depends on time and I set up a digital look-up table that
returns the same constant value for the range from t to t+dt, would I
effectively force odeint to do what I want?  Am I making any sense?
Is there a better way?

Thanks,

Ryan



More information about the SciPy-User mailing list