Behind and beyond the MATLAB ODE suite. (English) Zbl 0955.65055

Summary: The paper explains the concepts of order and absolute stability of numerical methods for solving systems of first-order ordinary differential equations (ODE) of the form \[ y'= f(t,y),\quad y(t_0)= y_0,\quad\text{where }f: \mathbb{R}\times \mathbb{R}^n\to \mathbb{R}^n, \] describes the phenomenon of problem stiffness, and reviews explicit Runge-Kutta methods, and explicit and implicit linear multistep methods. It surveys the five numerical methods contained in the MATLAB ODE suitable (three for nonstiff problems and two for stiff problems) to solve the above system, lists the available options, and uses the odedemo command to demonstrate the methods. One stiff ode code in MATLAB can solve more general equations of the form \(M(t)y'= f(t,y)\) provided the Mass option is on.


65L20 Stability and convergence of numerical methods for ordinary differential equations
34-04 Software, source code, etc. for problems pertaining to ordinary differential equations
65L05 Numerical methods for initial value problems involving ordinary differential equations
65L06 Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations
34A34 Nonlinear ordinary differential equations and systems
68W30 Symbolic computation and algebraic computation
Full Text: DOI


[1] Hairer, E.; Wanner, G., Solving ordinary differential equations II, stiff and differential-algebraic problems, (), 5-8
[2] Lambert, J.D., Numerical methods for ordinary differential equations, () · Zbl 0745.65049
[3] Dormand, J.R.; Prince, P.J., A family of embedded Runge-Kutta formulae, J. computational and applied mathematics, 6, 2, 19-26, (1980) · Zbl 0448.65045
[4] Hairer, E.; Wanner, G., On the instability of the BDF formulas, SIAM J. numer. anal., 20, 6, 1206-1209, (1983) · Zbl 0538.65051
[5] Shampine, L.F.; Reichelt, M.W., The Matlab ODE suite, SIAM J. sci. comput., 18, 1, 1-22, (1997) · Zbl 0868.65040
[6] Ashino, R.; Vaillancourt, R., Hayawakari Matlab (introduction to Matlab), (1997), Kyoritsu Shuppan Tokyo, 6^{th} printing, (1999) (in Japanese); Korean translation, (1998)
[7] Using MATLAB, (1997), The MathWorks Natick, MA, Chapter 8
[8] Shampine, L.F.; Gordon, M.K., Computer solution of ordinary differential equations, (1975), W.H. Freeman San Francisco · Zbl 0347.65001
[9] Hull, T.E.; Enright, W.H.; Fellen, B.M.; Sedgwick, A.E., Comparing numerical methods for ordinary differential equations, SIAM J. numer. anal., 9, 4, 603-637, (1972) · Zbl 0221.65115
[10] Shampine, L.F., Numerical solution of ordinary differential equations, (1994), Chapman & Hall New York · Zbl 0826.65082
[11] Enright, W.H.; Hull, T.E.; Lindberg, B., Comparing numerical methods for stiff systems of odes, Bit, 15, 1, 10-48, (1975) · Zbl 0301.65040
[12] Shampine, L.F., Measuring stiffness, Appl. numer. math., 1, 2, 107-119, (1985) · Zbl 0555.65050
[13] Enright, W.H.; Hull, T.E., Comparing numerical methods for the solution of stiff systems of ODEs arising in chemistry, (), 45-67
[14] ()
[15] van der Houwen, P.J., Construction of integration formulas for initial value problems, (1977), North-Holland Amsterdam · Zbl 0359.65057
[16] Hindmarsh, A.C.; Byrne, G.D., Applications of EPISODE: an experimental package for the integration of ordinary differential equations, (), 147-166
[17] Kahaner, D.; Moler, C.; Nash, S., Numerical methods and software, (1989), Prentice-Hall Englewood Cliffs, NJ · Zbl 0744.65002
[18] Shampine, L.F., Evaluation of a test set for stiff ODE solvers, ACM trans. math. soft., 7, 4, 409-420, (1981)
[19] Lambert, J.D., Computational methods in ordinary differential equations, (1973), Wiley London, Chapter 5 · Zbl 0258.65069
[20] Butcher, J.C., The numerical analysis of ordinary differential equations. Runge-Kutta and general linear methods, (1987), Wiley Chichester, Chapter 4 · Zbl 0616.65072
[21] Harier, E.; Nørsett, S.P.; Wanner, G., Solving ordinary differential equations I, nonstiff problems, (1987), Springer-Verlag Berlin, Section III.8 · Zbl 0638.65058
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.