[SciPy-user] Specific Odeint Problem
Ted Kord
teddy.kord at googlemail.com
Sun Feb 15 11:22:04 EST 2009
Hi
I have a specific problem for which I am unsure how to use odeint.
I'll try to outline the problem in steps so please bear with me.
1. I have a set of data in a file (one column only), say resistances.
2. I have a set of 5 time-dependent ODEs to solve.
3. At each time step, a specific resistance value from the
aforementioned data set must be used in the solution of the ODE set
e.g dydt = 4.56 * exp(resistance)/F
4. I am able to solve the ODEs using odeint using a time range from
'a' to 'b' according to the requirements of odeint, i.e.,
odeint(f, y0, t)
5. However, it's proving difficult to incorporate 'at each time step',
the specific resistances required.
I'd appreciate any help on this issue.
Thanks.
Ted Kord
P.S: The resistance values cannot be hard-coded in. There are 5001
values each for a specific time step.
More information about the SciPy-User
mailing list