To run the test problems type 

	>> main_problem

in MATLAB, choosing from the following list
of test problems (10 stiff, 1 non-stiff):

------------------------------------------------------------
1. main_testequation

The simple stiff test problem: u' = - l*u

a  = 6.1267
a0 = 1909.8424

gain: 311.7245

k0 = 1.0472e-03
------------------------------------------------------------
2. main_testsystem

Diagonal stiff test system with l1 = 100, l2 = 1000

a  = 18.3575
a0 = 1910.8165

gain: 104.0891

k0 = 1.0467e-03
------------------------------------------------------------
3. main_threescale

Diagonal stiff test system with three scales

a  = 17.8751
a0 = 1910.8166

gain: 106.8982

k0 = 1.0467e-03
------------------------------------------------------------
4. nonnormal

Non-normal stiff 2x2 system

a  = 16.8788
a0 = 3042.8564

gain: 180.2768

k0 = 6.5728e-04
------------------------------------------------------------
5. main_reaction

Stiff problem from Hairer-Wanner

a  = 594.1775
a0 = 2887.1203

gain: 4.8590

k0 = 6.9273e-04
------------------------------------------------------------
6. main_hires

Stiff test problem from the ODE test set (Lioen/Swart)

a  = 8.471
a0 = 283.6674

gain: 33.4869

k0 = 7.0505e-03
------------------------------------------------------------
7. main_akzonobel

Stiff test problem from the ODE test set (Lioen/Swart)

a  = 2.1074
a0 = 18.1225

gain: 8.5995

k0 = 1.1036e-01
------------------------------------------------------------
8.	main_petzold

First stiff problem from the book by Ascher and Petzold

a  = 97.1425
a0 = 137.1608

gain: 1.4120

k0 = 1.4581e-02
------------------------------------------------------------
9.	main_vanderpol

Van der Pol's equation

a  = 137.4475
a0 = 10240

gain: 74.5012

k0 = 1.9531e-04
------------------------------------------------------------
10. main_heat

The heat equation in 1d

a  = 3855.5399
a0 = 64742.5466

gain: 16.7921

k0 = 3.0892e-05
------------------------------------------------------------
11. main_sincos

A non-stiff to show that the method is ok.

a  = 70.5049
a0 = 70.5049

gain: 1
