An explicit two-step method for solving stiff systems of ordinary differential equations

In this paper we propose a numerical method to integrate stiff ordinary differential systems of the form Y′ = f(t Y)t ∊ [t 0 t N ]Y ∊ R m m positive integer, with Y(t 0) = Y 0. The method is an explicit two-step scheme, with variable coefficients, depending on a stability parameter. We prove that the scheme is of the first order in accuracy and that it shows good stability properties. We conclude by giving some numerical results obtained solving known stiff systems of differential equations.