在python中使用质量矩阵的DAE

Question

像

这样的系统

M * y'= y

其中M是单数，可以用Matlabs ode23轻松解决。在我的情况下，M是konstant。

我找到了PyDStools，应该可以在FORTRAN或C ++中使用DASSL。所有这些都非常复杂，不太直观。 （至少对我来说）

让某人遇到同样的挑战，可以指出一个好的开始点吗？

我想我做不了类似的事情： y'= np.linalg.inv（M）* y

并使用scipy.integrate和vode求解器这个等式？