PROCEEDINGS OF THE RUSSIAN HIGHER SCHOOL
ACADEMY OF SCIENCES

Abstract
In basis of  L-stabile (2,2)-method with second order of accuracy constructed algorithm for numerical solving Cauchy problem for stiff systems of ordinary differential equations. Intermediate computations realized by L-stable numerical formula. In this algorithm used freezing of Jacobi matrix . Produced results of numerical comparison of constructed algorithm with known methods.

References

1. Хайрер Э., Ваннер Г. Решение обыкновенных дифференциальных уравнений. Жесткие и дифференциально-алгебраические задачи. – М.: Мир, 1999. – 685 c.
2. Byrne G.D., Hindmarsh A.C. ODE solvers: a review of current and coming attractions // J. of Comput. Physics. – 1987. – № 70. – P. 1–62.
3. Rosenbrock H.H. Some general implicit processes for the numerical solution of differential equations // Computer. – 1963. – № 5. – P. 329–330.
4. Новиков Е.А., Шитов Ю.А., Шокин Ю.И. Одношаговые методы решения жестких систем // ДАН СССР. – 1988. – Т. 301. – № 6. – C. 1310–1314.
5. Новиков Е.А., Шитов Ю.А. Алгоритм интегрирования жестких систем на основе (m,k)-метода второго порядка точности с численным вычислением матрицы Якоби // Препринт №20, Красноярск: ВЦ СО АН СССР, 1988, 23c.
6. Kampowski W., Rentrop P., Schmidt W. Classication and numerical simulation of electric circuits // Surveys on Mathematics for Industry, 2(1). – 1992. – P. 23–65.
7. Mazzia F. and Magherini C. Test Set for Initial Value Problem Solvers, release 2.4. Department Mathematics, University of Bari and INdAM, Research Unit of Bari, February 2008. Available at http://pitagora.dm.uniba.it/~testset/problems/ringmod.php
8. Iavernaro F. Test Set for Initial Value Problem Solvers // Department of Mathematics, University of Bari, August 2003. Available at http://www.dm.uniba.it/ testset.