Novel Padé‐based algorithms for numerical integration of ODEs

The paper aims to propose several new approaches for the discrete‐time integration of stiff non‐linear differential equations.

The proposed approaches build on a method developed by the authors involving Padé approximates about each function sample. Both single‐ and multi‐step methods are suggested. The use of Richardson extrapolation is recommended for increasing efficiency.

The efficacy of the methods is shown using two examples and results are compared to a standard integration technique.

The paper shows that the methods are suitable for application in any field of science requiring efficient and accurate numerical solution of stiff differential equations.