Abstract
In Part I of this paper [SIAMJ. Numer. Anal., 25 (1988), pp. 1055–1073], a Riccati transformation method for solving linear boundary value problems for ordinary differential equations is presented. Here the implementation aspects of this method are discussed. In particular a double sweep integration process which includes a strategy for solving a matrix Riccati differential equation and forming a continuous representation of its solution is described. The efficacy of the method is demonstrated on a number of singular perturbation problems. Particular attention is paid to the choice of scaling for the variables and its effect on the solution method. The computational experience, while not based upon a refined implementation, indicates that the Riccati method may well prove to be competitive with the other general methods.
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call