Solving linear ordinary differential equations by exponentials of iterated commutators

Summary A sequence of transformations of a linear system of ordinary differential equations is investigated. It is shown that these transformations produce new systems which represent progressively smaller perturbations of the original set of equations.The transformations are implemented as a basis of a numerical method. Order, stability and error control of this method are analyzed. Numerical examples demonstrate the potential of this approach.