The approximate solutions (24)–(26) are generally, as shall be shown in the numerical experiments of this Letter, not valid for large t. A simple way of ensuring validity of the approximations for large t is to treat (17)–(22) as an algorithm for approximating the solutions of (1)–(4) in a sequence of intervals choosing the initial approximations as