Abstract
Algorithms based on the numerical inversion of Laplace transforms are developed for computing the solution x(λ) of a homogeneous ordinary linear differential system of the form dx/ = Ax. The local truncation error and the stability of the algorithms are analysed. It is shown how the algorithms can be made unconditionally stable and how the local truncation error can be minimised.

This publication has 1 reference indexed in Scilit: