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

Table 1

Numerical results for problem (1.1).

TOL	Method	TS	AVEE	MAXE	Time
10^−2	5oBBDF	21	1.09720e−005	7.3154e−005	0.01067
	ode15s	29	1.400e−003	8.700e−003	0.0313
	ode23s	19	1.200e−003	4.800e−003	0.2813
10^−4	5oBBDF	37	1.88948e−007	9.1173e−007	0.02234
	ode15s	61	3.4444e−005	1.7046e−004	0.0313
	ode23s	43	1.2509e−004	2.7400e−004	0.0156
10^−6	5oBBDF	79	1.72338e−009	8.8279e−009	0.009139
	ode15s	96	9.0543e−007	2.7175e−006	0.0469
	ode23s	148	8.0523e−006	1.3309e−005	0.0313