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

[Scott Askey]

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