First-order stiff ordinary differential equation solver

Runs 20 implicit and semi-implicit methods for first-order initial value stiff ODEs.

Implicit methods are important in numerical solution of stiff ODEs due to the instable nature of explicit methods. We wrote a program of a general one-dimensional implicit-solver including diagonally implicit case.