Abstract:
The method for solving stiff systems of ordinary differential equations based on calculation of matrix exponent is introduced. In addition, the method of accuracy control for linear problems is developed. Computational results for linear and non-linear problems and their dependencies on problem and method parameters are shown.
Keywords:chemical kinetics, stiff systems, matrix exponent, exponential integrator, Van der Pol oscillator.