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