On the numerical solution of stiff systems Nuran Guzel, Mustafa Bayram * Yildiz Teknik Universitesi, Fen-Edebiyat Fakultesi, Matematik Bolumu, 34210 Davutpasa-Esenler-Istanbul, Turkey Abstract In this paper, we use power series method to solve stiff ordinary differential equations of the first order and an ordinary differential equation of any order by converting it into a system of differential of the order one. Theoretical considerations has been discussed and some examples were presented to show the ability of the method for linear and non- linear systems of differential equations. We use MAPLE computer algebra systems for numerical calculations [G. Frank, MAPLE V, CRC Press Inc., Florida, 1996]. Ó 2005 Elsevier Inc. All rights reserved. Keywords: Stiff system; Power series; Numerical solution of the ordinary differential equations; MAPLE 1. Introduction A system of first order differential equations can be considered as: y 0 1 ¼ f 1 ðx; y 1 ; ... ; y n Þ y 0 2 ¼ f 2 ðx; y 1 ; ... ; y n Þ . . . y 0 n ¼ f n ðx; y 1 ; ... ; y n Þ 8 > > > < > > > : ð1:1Þ 0096-3003/$ - see front matter Ó 2005 Elsevier Inc. All rights reserved. doi:10.1016/j.amc.2004.11.035 * Corresponding author. E-mail address: msbayram@yildiz.edu.tr (M. Bayram). Applied Mathematics and Computation 170 (2005) 230–236 www.elsevier.com/locate/amc