where
is a vector of dependent
variables and . Differential-algebraic equations can be solved numerically
in the Wolfram Language using the
command NDSolve,
and some can be solved exactly with DSolve.
For suitable systems, differentiating the algebraic constraints can yield relations used to reduce the system to ordinary
differential equations. The original constraints must still be imposed, since
differentiation can introduce additional solutions.
Ascher, U. and Petzold, L. Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations.
Philadelphia, PA: SIAM, 1998.Brenan, K.; Campbell, S.; and Petzold,
L. Numerical Solutions of Initial-Value Problems in Differential-Algebraic Equations.
New York: Elsevier, 1989.Brown, P. N.; Hindmarsh, A. C.; and
Petzold, L. R. "Using Krylov Methods in the Solution of Large-Scale Differential-Algebraic
Systems." SIAM J. Sci. Comput.15, 1467-1488, 1994.Brown,
P. N.; Hindmarsh, A. C.; and Petzold, L. R. "Consistent Initial
Condition Calculation for Differential-Algebraic Systems." SIAM J. Sci. Comput.19,
1495-1512, 1998.Deuflhard, P.; Hairer, E. and Zugck, J. "One-Step
and Extrapolation Methods for Differential-Algebraic Systems." Numer. Math.51,
501-516, 1987.Hairer, E. and Lubich, C. On Extrapolation Methods
for Stiff and Differential-Algebraic Equations. Leipzig, Germany: Teubner, 1988.Hairer,
E. and Wanner, G. Solving
Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems, 2nd
ed. Berlin, Germany: Springer-Verlag, 1996.Hindmarsh, A. and
Taylor, A. "User Documentation for IDA: A Differential-Algebraic Equation Solver
for Sequential and Parallel Computers." Lawrence Livermore National Laboratory
report, UCRL-MA-136910, December 1999.