Picard Iteration method, Chebyshev Polynomial Approximation, and Global Numerical Integration of Dynamical Motions

Abstract

view Abstract Citations (12) References (5) Co-Reads Similar Papers Volume Content Graphics Metrics Export Citation NASA/ADS Picard Iteration method, Chebyshev Polynomial Approximation, and Global Numerical Integration of Dynamical Motions Fukushima, Toshio Abstract The Picard iteration method and the Chebyshev polynomial approximation were combined to obtain numerically a global solution of ordinary differential equations. The method solves both the initial and boundary value problems. The method directly provides not the tabulated values of the solution but the polynomials interpolating the solution in the integration interval given. In the case of scalar computation, the method is a few to several times as fast as the multistep method when (1) a good approximation of the solution is known beforehand, (2) the right hand members of differential equations are weakly dependent on the solution, and/or (3) the magnitude of the right hand members are small. Thus the method is suitable for (1) orbital improvements, (2) the integration of almost uniform rotations, and (3) perturbed dynamics in general. The method will be greatly accelerated by using vector/parallel computers since its main part is the numerical quadrature of known function of time. Publication: The Astronomical Journal Pub Date: May 1997 DOI: 10.1086/118404 Bibcode: 1997AJ....113.1909F full text sources ADS |