Code covered by the BSD License

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

• ex1 This is a singular perturbation problem involving Michaelis-Menten
• ex1 This is a singular perturbation problem involving Michaelis-Menten
• ex10 This is the DAS 2 test problem of I. Cameron, Solution of differential-
• ex10 This is the DAS 2 test problem of I. Cameron, Solution of differential-
• ex10ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
• ex11 This is a mechanics problem studied in section 4 of Chapter XII of
• ex11 This is a mechanics problem studied in section 4 of Chapter XII of
• ex11ode(t,u,flag) Copyright 1999, The MathWorks, Inc.
• ex12
• ex12
• ex12ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
• ex13 This is the discharge pressure control problem from E. Hairer, Ch. Lubich,
• ex13 This is the discharge pressure control problem from E. Hairer, Ch. Lubich,
• ex13ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
• ex14 DAE examples of index greater than 1
• ex14 DAE examples of index greater than 1
• ex14ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
• ex1ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
• ex2 This is an artificial problem with two sets of consistent initial
• ex2 This is an artificial problem with two sets of consistent initial
• ex2ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
• ex3 This is the pendulum problem posed as an index 1 problem. Formulated as
• ex3 This is the pendulum problem posed as an index 1 problem. Formulated as
• ex3ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
• ex4 The Robertson problem is a classic test problem for codes that solve
• ex4 The Robertson problem is a classic test problem for codes that solve
• ex4ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
• ex5 This is the two phase plug flow problem from E. Hairer, Ch. Lubich,
• ex5 This is the two phase plug flow problem from E. Hairer, Ch. Lubich,
• ex5ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
• ex6 This is the two phase plug flow problem from E. Hairer, Ch. Lubich,
• ex6 This is the two phase plug flow problem from E. Hairer, Ch. Lubich,
• ex6ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
• ex7 This is the one transistor amplifier problem displayed on p. 377 of
• ex7 This is the one transistor amplifier problem displayed on p. 377 of
• ex7ode(t,u,flag) Copyright 1999, The MathWorks, Inc.
• ex8 This is the two transistor amplifier problem of E. Hairer, C. Lubich,
• ex8 This is the two transistor amplifier problem of E. Hairer, C. Lubich,
• ex8ode(t,u,flag) Copyright 1999, The MathWorks, Inc.
• ex9 This is the DAS 1 test problem of I. Cameron, Solution of differential-
• ex9 This is the DAS 1 test problem of I. Cameron, Solution of differential-
• ex9ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
• ex1.mThis is a singular perturbation problem involving Michaelis-Menten
• ex10.mThis is the DAS 2 test problem of I. Cameron, Solution of differential-
• ex11.mThis is a mechanics problem studied in section 4 of Chapter XII of
• ex12.m
• ex13.mThis is the discharge pressure control problem from E. Hairer, Ch. Lubich,
• ex14.mDAE examples of index greater than 1
• ex2.mThis is an artificial problem with two sets of consistent initial
• ex3.mThis is the pendulum problem posed as an index 1 problem. Formulated as
• ex4.mThe Robertson problem is a classic test problem for codes that solve
• ex5.mThis is the two phase plug flow problem from E. Hairer, Ch. Lubich,
• ex6.mThis is the two phase plug flow problem from E. Hairer, Ch. Lubich,
• ex7.mThis is the one transistor amplifier problem displayed on p. 377 of
• ex8.mThis is the two transistor amplifier problem of E. Hairer, C. Lubich,
• ex9.mThis is the DAS 1 test problem of I. Cameron, Solution of differential-
• View all files

3.66667

3.7 | 3 ratings Rate this file 24 Downloads (last 30 days) File Size: 308 KB File ID: #7481

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

19 Apr 2005 (Updated )

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

File Information
Description

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)

MATLAB release MATLAB 5.3 (R11)
27 Feb 2013
27 Oct 2007
21 Nov 2005

Dear sir
I am simulating packed tower steady-state
by jacobian but don converge and need to
Best regard

08 Aug 2005

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)