MathGroup Archive 1995

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: DAE Solver in Mathematica?

In <3npp74$ohp at> achillgc at (Geoff C. Achilles) writes:

>I am trying to solve a system of differential and algebraic equations
>using Mathematica.  NDSolve seems like the logical function, but it
>doesn't seem to like algebraic equations.  I post a trivial example
>below to illustrate my problem.  Converting the algebraic equation to a
>differential equation is an alternative (but is difficult in my
>particular case).  Is there a package out there for solving DAE's?

>Example calculation:

>NDSolve[ { T'[t] == -c[t],
>           c[t]  == 1,
>           T[0]  == 400,
>           c[0]  == 1 },
>         { T[t], c[t] }, { t, 0, 10 } ]

>Obviously my case is more complicated and the algebraic equations
>(c[t]) are time dependent, but the basic setup is the same.  This
>particular format gives and error message saying that c[t] isn't in the
>form of an ODE - because it isn't an ODE.  Any suggestions?

>Geoff (achillgc at

Yes, there are packages for solving DAEs and you can get them from
netlib. The most famous one is the one written by Linda Petzold.
Send mail to netlib at for information.

Mathematica will have DAEs implemented in some future version as well
but not in the current one (V 2.2).

Sirpa Saarinen
saarinen at

  • Prev by Date: More on "Extracting data points from Plot[]"
  • Next by Date: Re: How to slow down plotting of curves?
  • Previous by thread: More on "Extracting data points from Plot[]"
  • Next by thread: DAE Solver in Mathematica?