File Exchange

image thumbnail

Manuscript of 'Solving Index-1 DAEs in MATLAB and Simulink'

version (308 KB) by Jacek Kierzenka
Manuscript of the SIAM Review paper + M-files for all examples.


Updated 01 Sep 2016

View License

This is a manuscript of the paper
Lawrence F. Shampine, Mark W. Reichelt, and Jacek A. Kierzenka,

Solving Index-1 DAEs in MATLAB and Simulink,
SIAM Review, Vol. 18, No. 3, pp. 538-552, 1999.

The paper describes mathematical and software developments needed for effective solution of differential algebraic equations of index 1. Modifications introduced to ODE15S and ODE23T in MATLAB 5.3 allow for solving linearly-implicit DAEs of index 1.

The M-files for all the original examples described in the paper are available in the directory DAE_examples.

The directory DAE_examples_65 contains the files for the examples updated to take advantage of features available in MATLAB 6.5 (R13). Correspondingly, DAE_examples_70 contains files modified for use with MATLAB 7.0 (R14).

The paper describes solving linearly-implicit DAEs with time-dependent, singular mass matrix
M(t)*y' = f(t,y)

In MATLAB 6.0, the solvers were extended to handle DAEs with state-dependent mass matrices

Cite As

Jacek Kierzenka (2020). Manuscript of 'Solving Index-1 DAEs in MATLAB and Simulink' (, MATLAB Central File Exchange. Retrieved .

Comments and Ratings (5)


Luis Enrique Avendaño

valialah zareei

Dear sir
I am simulating packed tower steady-state
by jacobian but don converge and need to
dynamic simulation ,Please help me
Best regard

George Chamoun

the examples only described systems of the form M(t)*y' = f(t,y).

there were no state-dependent mass matrix examples of the form M(t,y)*y'=f(t,y)


Updated license

Copyright update

MATLAB Release Compatibility
Created with R11
Compatible with any release
Platform Compatibility
Windows macOS Linux