... Press), pp. 2 57 27 5 ; 1 979 ,Applied Mathematics and Computation,vol.5,pp. 93– 121 . [3]Rice, J.R. 1983,Numerical Methods, Software, and Analysis(New York: McGraw-Hill),§9 .2. 16 .2 Adaptive ... organization about it:k1= hf(xn,yn)k 2 =hf(xn+h 2 ,yn+k1 2 )k3=hf(xn+h 2 ,yn+k 2 2)k4=hf(xn+ h, yn+ k3)yn+1= yn+k16+k 2 3+k33+k46+ O(h5)(16.1.3)The ... the whole interval. Figure 16.1 .2 illustrates theidea. In equations,k1= hf(xn,yn)k 2 =hfxn+1 2 h, yn+1 2 k1yn+1= yn+ k 2 + O(h3)(16.1 .2) As indicated in the error term,...