Adaptive nonlinear boundary value

Solution of the equations in (2.8-2.18) proceeds with an adaptive nonlinear boundary value method. The solution procedure has been discussed in detail elsewhere (10) and we outline only the essential features here. Our goal is to obtain a discrete solution of the governing equations on the mesh Af... 

Twopnt, solves systems of nonlinear, stiff, boundary-value problems [158]. It implements a hybrid, damped, modified Newton algorithm [159]. Local error is controlled using adaptive placement of mesh intervals. 

---