[SciPy-user] Differential Algebraic Equation Solvers

Johann Rohwer jr at sun.ac.za
Wed Jul 30 06:43:12 EDT 2008


On Wednesday, 30 July 2008, Nils Wagner wrote:
> Hi all,
>
> I am looking for a python tool to solve a set
> of DAE's, e.g.
>
> m \ddot{x} = 2 \lambda x
> m \ddot{y} = 2 \lambda y
> x^2 + y^2 - l^2 = 0

Our group has developed PySUNDIALS 
(http://pysundials.sourceforge.net), which is a Python package 
providing bindings for the SUNDIALS suite of solvers (using ctypes). 
SUNDIALS has a DAE solver (IDA), which should be able to handle your 
kind of problem. The PySUNDIALS distributrion comes with examples for 
each of the wrapped solvers.

Johann



More information about the SciPy-User mailing list