Research Article

A Quantitative Comparison of Numerical Method for Solving Stiff Ordinary Differential Equations

Table 3

Numerical results for problem (2.2).

TOLMethodTSsAverage errorMAXETime

10−25oBBDF24 3 . 1 6 7 6 2 𝑒 0 0 5 2 . 1 5 6 8 𝑒 0 0 4 0.01150
ode15s34 1 . 0 9 0 𝑒 0 0 2 1 . 0 9 0 𝑒 0 0 2 0.0313
ode23s21 1 . 0 0 4 𝑒 0 0 3 6 . 4 0 0 𝑒 0 0 3 0.0313
10−45oBBDF48 4 . 3 9 6 3 4 𝑒 0 0 7 1 . 8 6 5 2 𝑒 0 0 6 0.033346
ode15s71 5 . 7 5 0 1 𝑒 0 0 5 2 . 0 2 7 4 𝑒 0 0 4 0.0469
ode23s55 1 . 3 9 5 3 𝑒 0 0 4 3 . 5 3 7 5 𝑒 0 0 4 0.0156
10−65oBBDF127 1 . 4 1 7 1 1 𝑒 0 0 8 2 . 0 6 2 9 𝑒 0 0 9 0.009460
ode15s140 8 . 3 6 8 8 𝑒 0 0 7 4 . 0 0 5 9 𝑒 0 0 6 0.0469
ode23s233 7 . 2 8 2 4 𝑒 0 0 6 1 . 7 0 2 3 𝑒 0 0 5 0.1250