The numerical solution of unstable ordinary differential equations

摘要

Invariant imbedding, or the Riccati transformation, has been used to solve unstable ordinary differential equations for a few years. This paper compares the above method with parallel or multiple shooting and a method using Chebyshev series. Parallel shooting gives a solution as accurate as that obtained using the Riccati transformation, in a comparable time.