[SciPy-user] integrate.odeint , stiff chemical equations and mass conservation -any hint?

ms devicerandom@gmail....
Wed Jun 3 09:27:05 CDT 2009


Rob Clewley ha scritto:
> FYI if you conserve something explicitly in your system using an
> algebraic constraint then you have to remember to reduce your list of
> differential equations accordingly. Imposing one constraint will
> eliminate one D.E. when you substitute it's effect into them. You'll
> have to work that out symbolically in advance of setting it up in any
> integrator. 

Right, trying to work it out.

> You will need to use Radau in pydstool to solve the
> resulting differential-algebraic eqn (DAE), 

Uh? Why a DAE?

thanks!
M.


More information about the SciPy-user mailing list