File Exchange

image thumbnail

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

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

21 Downloads

Updated 01 Sep 2016

View Version History

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).

NOTE:
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
M(t,y)*y'=f(t,y)

Cite As

Jacek Kierzenka (2021). Manuscript of 'Solving Index-1 DAEs in MATLAB and Simulink' (https://www.mathworks.com/matlabcentral/fileexchange/7481-manuscript-of-solving-index-1-daes-in-matlab-and-simulink), MATLAB Central File Exchange. Retrieved .

Comments and Ratings (5)

Tucker McClure

mdg

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)

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

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!