[SciPy-user] Passing DAE mass matrix to integrate.odeint and/or ode

Scott Askey scotta_2002@yahoo....
Wed Jan 28 11:28:57 CST 2009


I am trying to solve and semi-explicit DAE system (index 1).
For a simple pendulum I 5 equations, there 1st order sytem and the explicit contraint.

0=x1^2+x2^2
x1'=x3
x2'=x4
x3'=-lam*x1
x4'=-lam*x2-g

Can I do this with scipy or must I try pydstools?  The problem is doable with  matlab ode15s which is simlar to  scipy.integrate.ode(f).\
set_integrator('vode', method='bdf', order=15) 

http://www.scipy.org/NumPy_for_Matlab_Users

V/R

Scott
 


      


More information about the SciPy-user mailing list