File Exchange

image thumbnail

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

version 1.1.0.1 (308 KB) by

Manuscript of the SIAM Review paper + M-files for all examples.

3.66667
3 Ratings

14 Downloads

Updated

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)

Comments and Ratings (4)

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)

Updates

1.1.0.1

Updated license

1.1

Copyright update

MATLAB Release
MATLAB 5.3 (R11)

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video