Comparison of biochemical simulations using integrators derived from "Adomian" decomposition with traditional numerical methods.

Comparisons are made between some traditional numerical integrators and integration using "Adomian" power series solutions to the ordinary differential equations. These are initial investigations to determine the viability of their application to the simulation of large complex metabolic pathways. A small set of test equations was employed to represent the types of problems encountered in biochemical applications. It was found that the "Adomian" method is as accurate as the numerical methods and, for 'nonstiff' equations or for small simulation times, the "Adomian" method is often more efficient. The results suggest that it may be worthwhile refining this method for biochemical simulations for situations where the traditional numerical methods fail.