A NEW APPROXIMATE METHOD FOR NUMERICALLY SOLVING OF ODE SYSTEMS IN CHEMICAL-KINETICS

In a recent paper we described an algorithm to calculate analytical fu nctions for the time-dependence of concentrations based on a system of ODEs by Laplace transforms. In this paper, we extend the procedure to the case in which the analytical solutions are finished with an addit ional change function. The method is suitable for models containing fi rst and second order steps. No differences between stiff and nonstiff ODEs are made. We show the application of our method to the case of pr opane conversion, a mechanism with 32 steps and 20 species. Comparison s with LARKIN (DEUFLHARD et al.) are made.