Formulas for numerical computation of derivatives f^(3) means third derivative f^(2)(x(1)) means second derivative computed at point x(1) f(x(2)) means function evaluated at point x(2) Points are x(0), x(1)=x(0)+h, x(2)=x0+2h, ... The error term gives the power of h and derivative order with z chosen for maximum value. The error terms, err=, are for test polynomial sum(x^pwr) with h = 0.125 or smaller, and order = 'pwr' computing order= 1, npoints= 2, at term= 1 f^( 1)(x( 1)= -1.00000000000000E+00 * f(x( 1)) + 1.00000000000000E+00 * f(x( 2) ) + O(h^ 1)f^( 2)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 0.0, x=1.0000, pwr= 0 err= 0.00000000000000E+00, h=0.1250, aderiv= 1.0, x=1.0000, pwr= 1 err=-1.25000000000000E-01, h=0.1250, aderiv= 2.0, x=1.0000, pwr= 2 err=-3.90625000000000E-01, h=0.1250, aderiv= 3.0, x=1.0000, pwr= 3 computing order= 1, npoints= 2, at term= 2 f^( 1)(x( 2)= -1.00000000000000E+00 * f(x( 1)) + 1.00000000000000E+00 * f(x( 2) ) + O(h^ 1)f^( 2)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 0.0, x=1.1250, pwr= 0 err= 0.00000000000000E+00, h=0.1250, aderiv= 1.0, x=1.1250, pwr= 1 err= 1.25000000000000E-01, h=0.1250, aderiv= 2.3, x=1.1250, pwr= 2 err= 4.06250000000000E-01, h=0.1250, aderiv= 3.8, x=1.1250, pwr= 3 computing order= 1, npoints= 3, at term= 1 f^( 1)(x( 1)= -1.50000000000000E+00 * f(x( 1)) + 2.00000000000000E+00 * f(x( 2) + -5.00000000000000E-01 * f(x( 3) ) + O(h^ 2)f^( 3)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 1.0, x=1.0000, pwr= 1 err= 0.00000000000000E+00, h=0.1250, aderiv= 2.0, x=1.0000, pwr= 2 err= 3.12500000000000E-02, h=0.1250, aderiv= 3.0, x=1.0000, pwr= 3 err= 1.36718750000000E-01, h=0.1250, aderiv= 4.0, x=1.0000, pwr= 4 computing order= 1, npoints= 3, at term= 2 f^( 1)(x( 2)= -5.00000000000000E-01 * f(x( 1)) + -0.00000000000000E+00 * f(x( 2) + 5.00000000000000E-01 * f(x( 3) ) + O(h^ 2)f^( 3)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 1.0, x=1.1250, pwr= 1 err= 0.00000000000000E+00, h=0.1250, aderiv= 2.3, x=1.1250, pwr= 2 err=-1.56250000000000E-02, h=0.1250, aderiv= 3.8, x=1.1250, pwr= 3 err=-7.03125000000000E-02, h=0.1250, aderiv= 5.7, x=1.1250, pwr= 4 computing order= 1, npoints= 3, at term= 3 f^( 1)(x( 3)= 5.00000000000000E-01 * f(x( 1)) + -2.00000000000000E+00 * f(x( 2) + 1.50000000000000E+00 * f(x( 3) ) + O(h^ 2)f^( 3)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 1.0, x=1.2500, pwr= 1 err= 0.00000000000000E+00, h=0.1250, aderiv= 2.5, x=1.2500, pwr= 2 err= 3.12500000000000E-02, h=0.1250, aderiv= 4.7, x=1.2500, pwr= 3 err= 1.44531250000000E-01, h=0.1250, aderiv= 7.8, x=1.2500, pwr= 4 computing order= 1, npoints= 4, at term= 1 f^( 1)(x( 1)= -1.83333333333333E+00 * f(x( 1)) + 3.00000000000000E+00 * f(x( 2) + -1.50000000000000E+00 * f(x( 3) + 3.33333333333333E-01 * f(x( 4) ) + O(h^ 3)f^( 4)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 2.0, x=1.0000, pwr= 2 err=-1.77635683940025E-15, h=0.1250, aderiv= 3.0, x=1.0000, pwr= 3 err=-1.17187500000018E-02, h=0.1250, aderiv= 4.0, x=1.0000, pwr= 4 err=-6.73828125000018E-02, h=0.1250, aderiv= 5.0, x=1.0000, pwr= 5 computing order= 1, npoints= 4, at term= 2 f^( 1)(x( 2)= -3.33333333333333E-01 * f(x( 1)) + -5.00000000000000E-01 * f(x( 2) + 1.00000000000000E+00 * f(x( 3) + -1.66666666666667E-01 * f(x( 4) ) + O(h^ 3)f^( 4)(z) err=-8.88178419700125E-16, h=0.1250, aderiv= 2.3, x=1.1250, pwr= 2 err=-8.88178419700125E-16, h=0.1250, aderiv= 3.8, x=1.1250, pwr= 3 err= 3.90624999999911E-03, h=0.1250, aderiv= 5.7, x=1.1250, pwr= 4 err= 2.29492187499991E-02, h=0.1250, aderiv= 8.0, x=1.1250, pwr= 5 computing order= 1, npoints= 4, at term= 3 f^( 1)(x( 3)= 1.66666666666667E-01 * f(x( 1)) + -1.00000000000000E+00 * f(x( 2) + 5.00000000000000E-01 * f(x( 3) + 3.33333333333333E-01 * f(x( 4) ) + O(h^ 3)f^( 4)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 2.5, x=1.2500, pwr= 2 err= 0.00000000000000E+00, h=0.1250, aderiv= 4.7, x=1.2500, pwr= 3 err=-3.90625000000000E-03, h=0.1250, aderiv= 7.8, x=1.2500, pwr= 4 err=-2.34375000000000E-02, h=0.1250, aderiv= 12.2, x=1.2500, pwr= 5 computing order= 1, npoints= 4, at term= 4 f^( 1)(x( 4)= -3.33333333333333E-01 * f(x( 1)) + 1.50000000000000E+00 * f(x( 2) + -3.00000000000000E+00 * f(x( 3) + 1.83333333333333E+00 * f(x( 4) ) + O(h^ 3)f^( 4)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 2.8, x=1.3750, pwr= 2 err= 0.00000000000000E+00, h=0.1250, aderiv= 5.7, x=1.3750, pwr= 3 err= 1.17187500000071E-02, h=0.1250, aderiv= 10.4, x=1.3750, pwr= 4 err= 7.17773437500142E-02, h=0.1250, aderiv= 17.9, x=1.3750, pwr= 5 computing order= 1, npoints= 5, at term= 1 f^( 1)(x( 1)= -2.08333333333333E+00 * f(x( 1)) + 4.00000000000000E+00 * f(x( 2) + -3.00000000000000E+00 * f(x( 3) + 1.33333333333333E+00 * f(x( 4) + -2.50000000000000E-01 * f(x( 5) ) + O(h^ 4)f^( 5)(z) err= 3.55271367880050E-15, h=0.1250, aderiv= 3.0, x=1.0000, pwr= 3 err= 7.10542735760100E-15, h=0.1250, aderiv= 4.0, x=1.0000, pwr= 4 err= 5.85937500000711E-03, h=0.1250, aderiv= 5.0, x=1.0000, pwr= 5 err= 4.24804687500142E-02, h=0.1250, aderiv= 6.0, x=1.0000, pwr= 6 computing order= 1, npoints= 5, at term= 2 f^( 1)(x( 2)= -2.50000000000000E-01 * f(x( 1)) + -8.33333333333333E-01 * f(x( 2) + 1.50000000000000E+00 * f(x( 3) + -5.00000000000000E-01 * f(x( 4) + 8.33333333333333E-02 * f(x( 5) ) + O(h^ 4)f^( 5)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 3.8, x=1.1250, pwr= 3 err= 0.00000000000000E+00, h=0.1250, aderiv= 5.7, x=1.1250, pwr= 4 err=-1.46484375000000E-03, h=0.1250, aderiv= 8.0, x=1.1250, pwr= 5 err=-1.08032226562500E-02, h=0.1250, aderiv= 10.8, x=1.1250, pwr= 6 computing order= 1, npoints= 5, at term= 3 f^( 1)(x( 3)= 8.33333333333333E-02 * f(x( 1)) + -6.66666666666667E-01 * f(x( 2) + 0.00000000000000E+00 * f(x( 3) + 6.66666666666667E-01 * f(x( 4) + -8.33333333333333E-02 * f(x( 5) ) + O(h^ 4)f^( 5)(z) err= 8.88178419700125E-16, h=0.1250, aderiv= 4.7, x=1.2500, pwr= 3 err= 0.00000000000000E+00, h=0.1250, aderiv= 7.8, x=1.2500, pwr= 4 err= 9.76562500000000E-04, h=0.1250, aderiv= 12.2, x=1.2500, pwr= 5 err= 7.32421875000711E-03, h=0.1250, aderiv= 18.3, x=1.2500, pwr= 6 computing order= 1, npoints= 5, at term= 4 f^( 1)(x( 4)= -8.33333333333333E-02 * f(x( 1)) + 5.00000000000000E-01 * f(x( 2) + -1.50000000000000E+00 * f(x( 3) + 8.33333333333333E-01 * f(x( 4) + 2.50000000000000E-01 * f(x( 5) ) + O(h^ 4)f^( 5)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 5.7, x=1.3750, pwr= 3 err= 0.00000000000000E+00, h=0.1250, aderiv= 10.4, x=1.3750, pwr= 4 err=-1.46484375000355E-03, h=0.1250, aderiv= 17.9, x=1.3750, pwr= 5 err=-1.11694335937571E-02, h=0.1250, aderiv= 29.5, x=1.3750, pwr= 6 computing order= 1, npoints= 5, at term= 5 f^( 1)(x( 5)= 2.50000000000000E-01 * f(x( 1)) + -1.33333333333333E+00 * f(x( 2) + 3.00000000000000E+00 * f(x( 3) + -4.00000000000000E+00 * f(x( 4) + 2.08333333333333E+00 * f(x( 5) ) + O(h^ 4)f^( 5)(z) err=-7.10542735760100E-15, h=0.1250, aderiv= 6.8, x=1.5000, pwr= 3 err= 0.00000000000000E+00, h=0.1250, aderiv= 13.5, x=1.5000, pwr= 4 err= 5.85937499998579E-03, h=0.1250, aderiv= 25.3, x=1.5000, pwr= 5 err= 4.54101562500000E-02, h=0.1250, aderiv= 45.6, x=1.5000, pwr= 6 computing order= 1, npoints= 6, at term= 1 f^( 1)(x( 1)= -2.28333333333333E+00 * f(x( 1)) + 5.00000000000000E+00 * f(x( 2) + -5.00000000000000E+00 * f(x( 3) + 3.33333333333333E+00 * f(x( 4) + -1.25000000000000E+00 * f(x( 5) + 2.00000000000000E-01 * f(x( 6) ) + O(h^ 5)f^( 6)(z) err=-7.10542735760100E-15, h=0.1250, aderiv= 4.0, x=1.0000, pwr= 4 err=-2.13162820728030E-14, h=0.1250, aderiv= 5.0, x=1.0000, pwr= 5 err=-3.66210937502132E-03, h=0.1250, aderiv= 6.0, x=1.0000, pwr= 6 err=-3.25012207031605E-02, h=0.1250, aderiv= 7.0, x=1.0000, pwr= 7 computing order= 1, npoints= 6, at term= 2 f^( 1)(x( 2)= -2.00000000000000E-01 * f(x( 1)) + -1.08333333333333E+00 * f(x( 2) + 2.00000000000000E+00 * f(x( 3) + -1.00000000000000E+00 * f(x( 4) + 3.33333333333333E-01 * f(x( 5) + -5.00000000000000E-02 * f(x( 6) ) + O(h^ 5)f^( 6)(z) err=-1.77635683940025E-15, h=0.1250, aderiv= 5.7, x=1.1250, pwr= 4 err=-1.77635683940025E-15, h=0.1250, aderiv= 8.0, x=1.1250, pwr= 5 err= 7.32421875001776E-04, h=0.1250, aderiv= 10.8, x=1.1250, pwr= 6 err= 6.59179687500178E-03, h=0.1250, aderiv= 14.2, x=1.1250, pwr= 7 computing order= 1, npoints= 6, at term= 3 f^( 1)(x( 3)= 5.00000000000000E-02 * f(x( 1)) + -5.00000000000000E-01 * f(x( 2) + -3.33333333333333E-01 * f(x( 3) + 1.00000000000000E+00 * f(x( 4) + -2.50000000000000E-01 * f(x( 5) + 3.33333333333333E-02 * f(x( 6) ) + O(h^ 5)f^( 6)(z) err=-8.88178419700125E-16, h=0.1250, aderiv= 7.8, x=1.2500, pwr= 4 err=-3.55271367880050E-15, h=0.1250, aderiv= 12.2, x=1.2500, pwr= 5 err=-3.66210937500000E-04, h=0.1250, aderiv= 18.3, x=1.2500, pwr= 6 err=-3.34167480469105E-03, h=0.1250, aderiv= 26.7, x=1.2500, pwr= 7 computing order= 1, npoints= 6, at term= 4 f^( 1)(x( 4)= -3.33333333333333E-02 * f(x( 1)) + 2.50000000000000E-01 * f(x( 2) + -1.00000000000000E+00 * f(x( 3) + 3.33333333333333E-01 * f(x( 4) + 5.00000000000000E-01 * f(x( 5) + -5.00000000000000E-02 * f(x( 6) ) + O(h^ 5)f^( 6)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 10.4, x=1.3750, pwr= 4 err= 0.00000000000000E+00, h=0.1250, aderiv= 17.9, x=1.3750, pwr= 5 err= 3.66210937500000E-04, h=0.1250, aderiv= 29.5, x=1.3750, pwr= 6 err= 3.38745117187500E-03, h=0.1250, aderiv= 47.3, x=1.3750, pwr= 7 computing order= 1, npoints= 6, at term= 5 f^( 1)(x( 5)= 5.00000000000000E-02 * f(x( 1)) + -3.33333333333333E-01 * f(x( 2) + 1.00000000000000E+00 * f(x( 3) + -2.00000000000000E+00 * f(x( 4) + 1.08333333333333E+00 * f(x( 5) + 2.00000000000000E-01 * f(x( 6) ) + O(h^ 5)f^( 6)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 13.5, x=1.5000, pwr= 4 err= 0.00000000000000E+00, h=0.1250, aderiv= 25.3, x=1.5000, pwr= 5 err=-7.32421875007105E-04, h=0.1250, aderiv= 45.6, x=1.5000, pwr= 6 err=-6.86645507812500E-03, h=0.1250, aderiv= 79.7, x=1.5000, pwr= 7 computing order= 1, npoints= 6, at term= 6 f^( 1)(x( 6)= -2.00000000000000E-01 * f(x( 1)) + 1.25000000000000E+00 * f(x( 2) + -3.33333333333333E+00 * f(x( 3) + 5.00000000000000E+00 * f(x( 4) + -5.00000000000000E+00 * f(x( 5) + 2.28333333333333E+00 * f(x( 6) ) + O(h^ 5)f^( 6)(z) err= 1.42108547152020E-14, h=0.1250, aderiv= 17.2, x=1.6250, pwr= 4 err= 2.84217094304040E-14, h=0.1250, aderiv= 34.9, x=1.6250, pwr= 5 err= 3.66210937500000E-03, h=0.1250, aderiv= 68.0, x=1.6250, pwr= 6 err= 3.47900390625568E-02, h=0.1250, aderiv= 128.9, x=1.6250, pwr= 7 computing order= 1, npoints= 7, at term= 1 f^( 1)(x( 1)= -2.45000000000000E+00 * f(x( 1)) + 6.00000000000000E+00 * f(x( 2) + -7.50000000000000E+00 * f(x( 3) + 6.66666666666667E+00 * f(x( 4) + -3.75000000000000E+00 * f(x( 5) + 1.20000000000000E+00 * f(x( 6) + -1.66666666666667E-01 * f(x( 7) ) + O(h^ 6)f^( 7)(z) err=-3.90798504668055E-14, h=0.1250, aderiv= 5.0, x=1.0000, pwr= 5 err=-4.26325641456060E-14, h=0.1250, aderiv= 6.0, x=1.0000, pwr= 6 err= 2.74658203126421E-03, h=0.1250, aderiv= 7.0, x=1.0000, pwr= 7 err= 2.91824340820455E-02, h=0.1250, aderiv= 8.0, x=1.0000, pwr= 8 computing order= 1, npoints= 7, at term= 2 f^( 1)(x( 2)= -1.66666666666667E-01 * f(x( 1)) + -1.28333333333333E+00 * f(x( 2) + 2.50000000000000E+00 * f(x( 3) + -1.66666666666667E+00 * f(x( 4) + 8.33333333333333E-01 * f(x( 5) + -2.50000000000000E-01 * f(x( 6) + 3.33333333333333E-02 * f(x( 7) ) + O(h^ 6)f^( 7)(z) err= 8.88178419700125E-15, h=0.1250, aderiv= 8.0, x=1.1250, pwr= 5 err= 1.06581410364015E-14, h=0.1250, aderiv= 10.8, x=1.1250, pwr= 6 err=-4.57763671867895E-04, h=0.1250, aderiv= 14.2, x=1.1250, pwr= 7 err=-4.92095947264204E-03, h=0.1250, aderiv= 18.2, x=1.1250, pwr= 8 computing order= 1, npoints= 7, at term= 3 f^( 1)(x( 3)= 3.33333333333333E-02 * f(x( 1)) + -4.00000000000000E-01 * f(x( 2) + -5.83333333333333E-01 * f(x( 3) + 1.33333333333333E+00 * f(x( 4) + -5.00000000000000E-01 * f(x( 5) + 1.33333333333333E-01 * f(x( 6) + -1.66666666666667E-02 * f(x( 7) ) + O(h^ 6)f^( 7)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 12.2, x=1.2500, pwr= 5 err= 1.42108547152020E-14, h=0.1250, aderiv= 18.3, x=1.2500, pwr= 6 err= 1.83105468760658E-04, h=0.1250, aderiv= 26.7, x=1.2500, pwr= 7 err= 1.99127197266336E-03, h=0.1250, aderiv= 38.1, x=1.2500, pwr= 8 computing order= 1, npoints= 7, at term= 4 f^( 1)(x( 4)= -1.66666666666667E-02 * f(x( 1)) + 1.50000000000000E-01 * f(x( 2) + -7.50000000000000E-01 * f(x( 3) + -0.00000000000000E+00 * f(x( 4) + 7.50000000000000E-01 * f(x( 5) + -1.50000000000000E-01 * f(x( 6) + 1.66666666666667E-02 * f(x( 7) ) + O(h^ 6)f^( 7)(z) err=-3.55271367880050E-15, h=0.1250, aderiv= 17.9, x=1.3750, pwr= 5 err=-7.10542735760100E-15, h=0.1250, aderiv= 29.5, x=1.3750, pwr= 6 err=-1.37329101562500E-04, h=0.1250, aderiv= 47.3, x=1.3750, pwr= 7 err=-1.51062011718750E-03, h=0.1250, aderiv= 74.3, x=1.3750, pwr= 8 computing order= 1, npoints= 7, at term= 5 f^( 1)(x( 5)= 1.66666666666667E-02 * f(x( 1)) + -1.33333333333333E-01 * f(x( 2) + 5.00000000000000E-01 * f(x( 3) + -1.33333333333333E+00 * f(x( 4) + 5.83333333333333E-01 * f(x( 5) + 4.00000000000000E-01 * f(x( 6) + -3.33333333333333E-02 * f(x( 7) ) + O(h^ 6)f^( 7)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 25.3, x=1.5000, pwr= 5 err=-1.42108547152020E-14, h=0.1250, aderiv= 45.6, x=1.5000, pwr= 6 err= 1.83105468750000E-04, h=0.1250, aderiv= 79.7, x=1.5000, pwr= 7 err= 2.03704833978691E-03, h=0.1250, aderiv= 136.7, x=1.5000, pwr= 8 computing order= 1, npoints= 7, at term= 6 f^( 1)(x( 6)= -3.33333333333333E-02 * f(x( 1)) + 2.50000000000000E-01 * f(x( 2) + -8.33333333333333E-01 * f(x( 3) + 1.66666666666667E+00 * f(x( 4) + -2.50000000000000E+00 * f(x( 5) + 1.28333333333333E+00 * f(x( 6) + 1.66666666666667E-01 * f(x( 7) ) + O(h^ 6)f^( 7)(z) err=-1.42108547152020E-14, h=0.1250, aderiv= 34.9, x=1.6250, pwr= 5 err=-1.42108547152020E-14, h=0.1250, aderiv= 68.0, x=1.6250, pwr= 6 err=-4.57763671875000E-04, h=0.1250, aderiv= 128.9, x=1.6250, pwr= 7 err=-5.14984130859375E-03, h=0.1250, aderiv= 239.4, x=1.6250, pwr= 8 computing order= 1, npoints= 7, at term= 7 f^( 1)(x( 7)= 1.66666666666667E-01 * f(x( 1)) + -1.20000000000000E+00 * f(x( 2) + 3.75000000000000E+00 * f(x( 3) + -6.66666666666667E+00 * f(x( 4) + 7.50000000000000E+00 * f(x( 5) + -6.00000000000000E+00 * f(x( 6) + 2.45000000000000E+00 * f(x( 7) ) + O(h^ 6)f^( 7)(z) err= 5.68434188608080E-14, h=0.1250, aderiv= 46.9, x=1.7500, pwr= 5 err=-5.68434188608080E-14, h=0.1250, aderiv= 98.5, x=1.7500, pwr= 6 err= 2.74658203125000E-03, h=0.1250, aderiv= 201.1, x=1.7500, pwr= 7 err= 3.12423706052414E-02, h=0.1250, aderiv= 402.1, x=1.7500, pwr= 8 computing order= 1, npoints= 8, at term= 1 f^( 1)(x( 1)= -2.59285714285714E+00 * f(x( 1)) + 7.00000000000000E+00 * f(x( 2) + -1.05000000000000E+01 * f(x( 3) + 1.16666666666667E+01 * f(x( 4) + -8.75000000000000E+00 * f(x( 5) + 4.20000000000000E+00 * f(x( 6) + -1.16666666666667E+00 * f(x( 7) + 1.42857142857143E-01 * f(x( 8) ) + O(h^ 7)f^( 8)(z) err= 8.52651282912120E-14, h=0.1250, aderiv= 6.0, x=1.0000, pwr= 6 err= 8.52651282912120E-14, h=0.1250, aderiv= 7.0, x=1.0000, pwr= 7 err=-2.40325927737217E-03, h=0.1250, aderiv= 8.0, x=1.0000, pwr= 8 err=-3.00407409665695E-02, h=0.1250, aderiv= 9.0, x=1.0000, pwr= 9 computing order= 1, npoints= 8, at term= 2 f^( 1)(x( 2)= -1.42857142857143E-01 * f(x( 1)) + -1.45000000000000E+00 * f(x( 2) + 3.00000000000000E+00 * f(x( 3) + -2.50000000000000E+00 * f(x( 4) + 1.66666666666667E+00 * f(x( 5) + -7.50000000000000E-01 * f(x( 6) + 2.00000000000000E-01 * f(x( 7) + -2.38095238095238E-02 * f(x( 8) ) + O(h^ 7)f^( 8)(z) err=-5.32907051820075E-15, h=0.1250, aderiv= 10.8, x=1.1250, pwr= 6 err= 8.88178419700125E-15, h=0.1250, aderiv= 14.2, x=1.1250, pwr= 7 err= 3.43322753899145E-04, h=0.1250, aderiv= 18.2, x=1.1250, pwr= 8 err= 4.33444976808772E-03, h=0.1250, aderiv= 23.1, x=1.1250, pwr= 9 computing order= 1, npoints= 8, at term= 3 f^( 1)(x( 3)= 2.38095238095238E-02 * f(x( 1)) + -3.33333333333333E-01 * f(x( 2) + -7.83333333333333E-01 * f(x( 3) + 1.66666666666667E+00 * f(x( 4) + -8.33333333333333E-01 * f(x( 5) + 3.33333333333333E-01 * f(x( 6) + -8.33333333333333E-02 * f(x( 7) + 9.52380952380952E-03 * f(x( 8) ) + O(h^ 7)f^( 8)(z) err=-3.55271367880050E-15, h=0.1250, aderiv= 18.3, x=1.2500, pwr= 6 err=-1.06581410364015E-14, h=0.1250, aderiv= 26.7, x=1.2500, pwr= 7 err=-1.14440917968750E-04, h=0.1250, aderiv= 38.1, x=1.2500, pwr= 8 err=-1.45912170410867E-03, h=0.1250, aderiv= 53.6, x=1.2500, pwr= 9 computing order= 1, npoints= 8, at term= 4 f^( 1)(x( 4)= -9.52380952380952E-03 * f(x( 1)) + 1.00000000000000E-01 * f(x( 2) + -6.00000000000000E-01 * f(x( 3) + -2.50000000000000E-01 * f(x( 4) + 1.00000000000000E+00 * f(x( 5) + -3.00000000000000E-01 * f(x( 6) + 6.66666666666667E-02 * f(x( 7) + -7.14285714285714E-03 * f(x( 8) ) + O(h^ 7)f^( 8)(z) err=-3.55271367880050E-15, h=0.1250, aderiv= 29.5, x=1.3750, pwr= 6 err=-7.10542735760100E-15, h=0.1250, aderiv= 47.3, x=1.3750, pwr= 7 err= 6.86645507954609E-05, h=0.1250, aderiv= 74.3, x=1.3750, pwr= 8 err= 8.84056091308594E-04, h=0.1250, aderiv= 115.0, x=1.3750, pwr= 9 computing order= 1, npoints= 8, at term= 5 f^( 1)(x( 5)= 7.14285714285714E-03 * f(x( 1)) + -6.66666666666667E-02 * f(x( 2) + 3.00000000000000E-01 * f(x( 3) + -1.00000000000000E+00 * f(x( 4) + 2.50000000000000E-01 * f(x( 5) + 6.00000000000000E-01 * f(x( 6) + -1.00000000000000E-01 * f(x( 7) + 9.52380952380952E-03 * f(x( 8) ) + O(h^ 7)f^( 8)(z) err= 7.10542735760100E-15, h=0.1250, aderiv= 45.6, x=1.5000, pwr= 6 err= 2.84217094304040E-14, h=0.1250, aderiv= 79.7, x=1.5000, pwr= 7 err=-6.86645507812500E-05, h=0.1250, aderiv= 136.7, x=1.5000, pwr= 8 err=-8.92639160156250E-04, h=0.1250, aderiv= 230.7, x=1.5000, pwr= 9 computing order= 1, npoints= 8, at term= 6 f^( 1)(x( 6)= -9.52380952380952E-03 * f(x( 1)) + 8.33333333333333E-02 * f(x( 2) + -3.33333333333333E-01 * f(x( 3) + 8.33333333333333E-01 * f(x( 4) + -1.66666666666667E+00 * f(x( 5) + 7.83333333333333E-01 * f(x( 6) + 3.33333333333333E-01 * f(x( 7) + -2.38095238095238E-02 * f(x( 8) ) + O(h^ 7)f^( 8)(z) err= 1.42108547152020E-14, h=0.1250, aderiv= 68.0, x=1.6250, pwr= 6 err= 0.00000000000000E+00, h=0.1250, aderiv= 128.9, x=1.6250, pwr= 7 err= 1.14440918054015E-04, h=0.1250, aderiv= 239.4, x=1.6250, pwr= 8 err= 1.50203704833984E-03, h=0.1250, aderiv= 437.6, x=1.6250, pwr= 9 computing order= 1, npoints= 8, at term= 7 f^( 1)(x( 7)= 2.38095238095238E-02 * f(x( 1)) + -2.00000000000000E-01 * f(x( 2) + 7.50000000000000E-01 * f(x( 3) + -1.66666666666667E+00 * f(x( 4) + 2.50000000000000E+00 * f(x( 5) + -3.00000000000000E+00 * f(x( 6) + 1.45000000000000E+00 * f(x( 7) + 1.42857142857143E-01 * f(x( 8) ) + O(h^ 7)f^( 8)(z) err= 2.84217094304040E-14, h=0.1250, aderiv= 98.5, x=1.7500, pwr= 6 err= 8.52651282912120E-14, h=0.1250, aderiv= 201.1, x=1.7500, pwr= 7 err=-3.43322753906250E-04, h=0.1250, aderiv= 402.1, x=1.7500, pwr= 8 err=-4.54902648925781E-03, h=0.1250, aderiv= 791.7, x=1.7500, pwr= 9 computing order= 1, npoints= 8, at term= 8 f^( 1)(x( 8)= -1.42857142857143E-01 * f(x( 1)) + 1.16666666666667E+00 * f(x( 2) + -4.20000000000000E+00 * f(x( 3) + 8.75000000000000E+00 * f(x( 4) + -1.16666666666667E+01 * f(x( 5) + 1.05000000000000E+01 * f(x( 6) + -7.00000000000000E+00 * f(x( 7) + 2.59285714285714E+00 * f(x( 8) ) + O(h^ 7)f^( 8)(z) err=-1.13686837721616E-13, h=0.1250, aderiv= 139.0, x=1.8750, pwr= 6 err= 0.00000000000000E+00, h=0.1250, aderiv= 304.2, x=1.8750, pwr= 7 err= 2.40325927688900E-03, h=0.1250, aderiv= 651.8, x=1.8750, pwr= 8 err= 3.21435928335632E-02, h=0.1250, aderiv=1374.8, x=1.8750, pwr= 9 computing order= 1, npoints= 9, at term= 1 f^( 1)(x( 1)= -2.71785714285714E+00 * f(x( 1)) + 8.00000000000000E+00 * f(x( 2) + -1.40000000000000E+01 * f(x( 3) + 1.86666666666667E+01 * f(x( 4) + -1.75000000000000E+01 * f(x( 5) + 1.12000000000000E+01 * f(x( 6) + -4.66666666666667E+00 * f(x( 7) + 1.14285714285714E+00 * f(x( 8) + -1.25000000000000E-01 * f(x( 9) ) + O(h^ 8)f^( 9)(z) err= 2.27373675443232E-13, h=0.1250, aderiv= 7.0, x=1.0000, pwr= 7 err= 6.82121026329696E-13, h=0.1250, aderiv= 8.0, x=1.0000, pwr= 8 err= 2.40325927870799E-03, h=0.1250, aderiv= 9.0, x=1.0000, pwr= 9 err= 3.48472595233034E-02, h=0.1250, aderiv= 10.0, x=1.0000, pwr= 10 computing order= 1, npoints= 9, at term= 2 f^( 1)(x( 2)= -1.25000000000000E-01 * f(x( 1)) + -1.59285714285714E+00 * f(x( 2) + 3.50000000000000E+00 * f(x( 3) + -3.50000000000000E+00 * f(x( 4) + 2.91666666666667E+00 * f(x( 5) + -1.75000000000000E+00 * f(x( 6) + 7.00000000000000E-01 * f(x( 7) + -1.66666666666667E-01 * f(x( 8) + 1.78571428571429E-02 * f(x( 9) ) + O(h^ 8)f^( 9)(z) err=-7.10542735760100E-15, h=0.1250, aderiv= 14.2, x=1.1250, pwr= 7 err= 4.26325641456060E-14, h=0.1250, aderiv= 18.2, x=1.1250, pwr= 8 err=-3.00407409582704E-04, h=0.1250, aderiv= 23.1, x=1.1250, pwr= 9 err=-4.39345836622351E-03, h=0.1250, aderiv= 28.9, x=1.1250, pwr= 10 computing order= 1, npoints= 9, at term= 3 f^( 1)(x( 3)= 1.78571428571429E-02 * f(x( 1)) + -2.85714285714286E-01 * f(x( 2) + -9.50000000000000E-01 * f(x( 3) + 2.00000000000000E+00 * f(x( 4) + -1.25000000000000E+00 * f(x( 5) + 6.66666666666667E-01 * f(x( 6) + -2.50000000000000E-01 * f(x( 7) + 5.71428571428571E-02 * f(x( 8) + -5.95238095238095E-03 * f(x( 9) ) + O(h^ 8)f^( 9)(z) err= 1.42108547152020E-14, h=0.1250, aderiv= 26.7, x=1.2500, pwr= 7 err= 1.42108547152020E-14, h=0.1250, aderiv= 38.1, x=1.2500, pwr= 8 err= 8.58306885049842E-05, h=0.1250, aderiv= 53.6, x=1.2500, pwr= 9 err= 1.26600265511456E-03, h=0.1250, aderiv= 74.5, x=1.2500, pwr= 10 computing order= 1, npoints= 9, at term= 4 f^( 1)(x( 4)= -5.95238095238095E-03 * f(x( 1)) + 7.14285714285714E-02 * f(x( 2) + -5.00000000000000E-01 * f(x( 3) + -4.50000000000000E-01 * f(x( 4) + 1.25000000000000E+00 * f(x( 5) + -5.00000000000000E-01 * f(x( 6) + 1.66666666666667E-01 * f(x( 7) + -3.57142857142857E-02 * f(x( 8) + 3.57142857142857E-03 * f(x( 9) ) + O(h^ 8)f^( 9)(z) err= 7.10542735760100E-15, h=0.1250, aderiv= 47.3, x=1.3750, pwr= 7 err= 2.84217094304040E-14, h=0.1250, aderiv= 74.3, x=1.3750, pwr= 8 err=-4.29153442382813E-05, h=0.1250, aderiv= 115.0, x=1.3750, pwr= 9 err=-6.38365745544434E-04, h=0.1250, aderiv= 175.7, x=1.3750, pwr= 10 computing order= 1, npoints= 9, at term= 5 f^( 1)(x( 5)= 3.57142857142857E-03 * f(x( 1)) + -3.80952380952381E-02 * f(x( 2) + 2.00000000000000E-01 * f(x( 3) + -8.00000000000000E-01 * f(x( 4) + 0.00000000000000E+00 * f(x( 5) + 8.00000000000000E-01 * f(x( 6) + -2.00000000000000E-01 * f(x( 7) + 3.80952380952381E-02 * f(x( 8) + -3.57142857142857E-03 * f(x( 9) ) + O(h^ 8)f^( 9)(z) err= 1.42108547152020E-14, h=0.1250, aderiv= 79.7, x=1.5000, pwr= 7 err=-2.84217094304040E-14, h=0.1250, aderiv= 136.7, x=1.5000, pwr= 8 err= 3.43322753622033E-05, h=0.1250, aderiv= 230.7, x=1.5000, pwr= 9 err= 5.14984130745688E-04, h=0.1250, aderiv= 384.4, x=1.5000, pwr= 10 computing order= 1, npoints= 9, at term= 6 f^( 1)(x( 6)= -3.57142857142857E-03 * f(x( 1)) + 3.57142857142857E-02 * f(x( 2) + -1.66666666666667E-01 * f(x( 3) + 5.00000000000000E-01 * f(x( 4) + -1.25000000000000E+00 * f(x( 5) + 4.50000000000000E-01 * f(x( 6) + 5.00000000000000E-01 * f(x( 7) + -7.14285714285714E-02 * f(x( 8) + 5.95238095238095E-03 * f(x( 9) ) + O(h^ 8)f^( 9)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 128.9, x=1.6250, pwr= 7 err=-5.68434188608080E-14, h=0.1250, aderiv= 239.4, x=1.6250, pwr= 8 err=-4.29153442382813E-05, h=0.1250, aderiv= 437.6, x=1.6250, pwr= 9 err=-6.49094581604004E-04, h=0.1250, aderiv= 790.1, x=1.6250, pwr= 10 computing order= 1, npoints= 9, at term= 7 f^( 1)(x( 7)= 5.95238095238095E-03 * f(x( 1)) + -5.71428571428571E-02 * f(x( 2) + 2.50000000000000E-01 * f(x( 3) + -6.66666666666667E-01 * f(x( 4) + 1.25000000000000E+00 * f(x( 5) + -2.00000000000000E+00 * f(x( 6) + 9.50000000000000E-01 * f(x( 7) + 2.85714285714286E-01 * f(x( 8) + -1.78571428571429E-02 * f(x( 9) ) + O(h^ 8)f^( 9)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 201.1, x=1.7500, pwr= 7 err=-5.68434188608080E-14, h=0.1250, aderiv= 402.1, x=1.7500, pwr= 8 err= 8.58306884765625E-05, h=0.1250, aderiv= 791.7, x=1.7500, pwr= 9 err= 1.30891799949495E-03, h=0.1250, aderiv=1539.4, x=1.7500, pwr= 10 computing order= 1, npoints= 9, at term= 8 f^( 1)(x( 8)= -1.78571428571429E-02 * f(x( 1)) + 1.66666666666667E-01 * f(x( 2) + -7.00000000000000E-01 * f(x( 3) + 1.75000000000000E+00 * f(x( 4) + -2.91666666666667E+00 * f(x( 5) + 3.50000000000000E+00 * f(x( 6) + -3.50000000000000E+00 * f(x( 7) + 1.59285714285714E+00 * f(x( 8) + 1.25000000000000E-01 * f(x( 9) ) + O(h^ 8)f^( 9)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 304.2, x=1.8750, pwr= 7 err= 0.00000000000000E+00, h=0.1250, aderiv= 651.8, x=1.8750, pwr= 8 err=-3.00407409667969E-04, h=0.1250, aderiv=1374.8, x=1.8750, pwr= 9 err=-4.61876392364502E-03, h=0.1250, aderiv=2864.3, x=1.8750, pwr= 10 computing order= 1, npoints= 9, at term= 9 f^( 1)(x( 9)= 1.25000000000000E-01 * f(x( 1)) + -1.14285714285714E+00 * f(x( 2) + 4.66666666666667E+00 * f(x( 3) + -1.12000000000000E+01 * f(x( 4) + 1.75000000000000E+01 * f(x( 5) + -1.86666666666667E+01 * f(x( 6) + 1.40000000000000E+01 * f(x( 7) + -8.00000000000000E+00 * f(x( 8) + 2.71785714285714E+00 * f(x( 9) ) + O(h^ 8)f^( 9)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 448.0, x=2.0000, pwr= 7 err= 0.00000000000000E+00, h=0.1250, aderiv=1024.0, x=2.0000, pwr= 8 err= 2.40325927734375E-03, h=0.1250, aderiv=2304.0, x=2.0000, pwr= 9 err= 3.72505187988281E-02, h=0.1250, aderiv=5120.0, x=2.0000, pwr= 10 computing order= 1, npoints= 10, at term= 1 f^( 1)(x( 1)= -2.82896825396825E+00 * f(x( 1)) + 9.00000000000000E+00 * f(x( 2) + -1.80000000000000E+01 * f(x( 3) + 2.80000000000000E+01 * f(x( 4) + -3.15000000000000E+01 * f(x( 5) + 2.52000000000000E+01 * f(x( 6) + -1.40000000000000E+01 * f(x( 7) + 5.14285714285714E+00 * f(x( 8) + -1.12500000000000E+00 * f(x( 9) + 1.11111111111111E-01 * f(x( 10) ) + O(h^ 9)f^( 10)(z) err=-1.08002495835535E-12, h=0.1250, aderiv= 8.0, x=1.0000, pwr= 8 err=-1.93267624126747E-12, h=0.1250, aderiv= 9.0, x=1.0000, pwr= 9 err=-2.70366668701172E-03, h=0.1250, aderiv= 10.0, x=1.0000, pwr= 10 err=-4.49484586770268E-02, h=0.1250, aderiv= 11.0, x=1.0000, pwr= 11 computing order= 1, npoints= 10, at term= 2 f^( 1)(x( 2)= -1.11111111111111E-01 * f(x( 1)) + -1.71785714285714E+00 * f(x( 2) + 4.00000000000000E+00 * f(x( 3) + -4.66666666666667E+00 * f(x( 4) + 4.66666666666667E+00 * f(x( 5) + -3.50000000000000E+00 * f(x( 6) + 1.86666666666667E+00 * f(x( 7) + -6.66666666666667E-01 * f(x( 8) + 1.42857142857143E-01 * f(x( 9) + -1.38888888888889E-02 * f(x( 10) ) + O(h^ 9)f^( 10)(z) err=-9.23705556488130E-14, h=0.1250, aderiv= 18.2, x=1.1250, pwr= 8 err= 2.41584530158434E-13, h=0.1250, aderiv= 23.1, x=1.1250, pwr= 9 err= 3.00407409440595E-04, h=0.1250, aderiv= 28.9, x=1.1250, pwr= 10 err= 5.03182411171110E-03, h=0.1250, aderiv= 35.7, x=1.1250, pwr= 11 computing order= 1, npoints= 10, at term= 3 f^( 1)(x( 3)= 1.38888888888889E-02 * f(x( 1)) + -2.50000000000000E-01 * f(x( 2) + -1.09285714285714E+00 * f(x( 3) + 2.33333333333333E+00 * f(x( 4) + -1.75000000000000E+00 * f(x( 5) + 1.16666666666667E+00 * f(x( 6) + -5.83333333333333E-01 * f(x( 7) + 2.00000000000000E-01 * f(x( 8) + -4.16666666666667E-02 * f(x( 9) + 3.96825396825397E-03 * f(x( 10) ) + O(h^ 9)f^( 10)(z) err=-2.84217094304040E-14, h=0.1250, aderiv= 38.1, x=1.2500, pwr= 8 err=-2.84217094304040E-14, h=0.1250, aderiv= 53.6, x=1.2500, pwr= 9 err=-7.51018524596248E-05, h=0.1250, aderiv= 74.5, x=1.2500, pwr= 10 err=-1.26734375955095E-03, h=0.1250, aderiv= 102.4, x=1.2500, pwr= 11 computing order= 1, npoints= 10, at term= 4 f^( 1)(x( 4)= -3.96825396825397E-03 * f(x( 1)) + 5.35714285714286E-02 * f(x( 2) + -4.28571428571429E-01 * f(x( 3) + -6.16666666666667E-01 * f(x( 4) + 1.50000000000000E+00 * f(x( 5) + -7.50000000000000E-01 * f(x( 6) + 3.33333333333333E-01 * f(x( 7) + -1.07142857142857E-01 * f(x( 8) + 2.14285714285714E-02 * f(x( 9) + -1.98412698412698E-03 * f(x( 10) ) + O(h^ 9)f^( 10)(z) err= 2.84217094304040E-14, h=0.1250, aderiv= 74.3, x=1.3750, pwr= 8 err= 0.00000000000000E+00, h=0.1250, aderiv= 115.0, x=1.3750, pwr= 9 err= 3.21865082355544E-05, h=0.1250, aderiv= 175.7, x=1.3750, pwr= 10 err= 5.47170639094929E-04, h=0.1250, aderiv= 265.7, x=1.3750, pwr= 11 computing order= 1, npoints= 10, at term= 5 f^( 1)(x( 5)= 1.98412698412698E-03 * f(x( 1)) + -2.38095238095238E-02 * f(x( 2) + 1.42857142857143E-01 * f(x( 3) + -6.66666666666667E-01 * f(x( 4) + -2.00000000000000E-01 * f(x( 5) + 1.00000000000000E+00 * f(x( 6) + -3.33333333333333E-01 * f(x( 7) + 9.52380952380952E-02 * f(x( 8) + -1.78571428571429E-02 * f(x( 9) + 1.58730158730159E-03 * f(x( 10) ) + O(h^ 9)f^( 10)(z) err= 2.84217094304040E-14, h=0.1250, aderiv= 136.7, x=1.5000, pwr= 8 err=-2.84217094304040E-14, h=0.1250, aderiv= 230.7, x=1.5000, pwr= 9 err=-2.14576721191406E-05, h=0.1250, aderiv= 384.4, x=1.5000, pwr= 10 err=-3.67462635040283E-04, h=0.1250, aderiv= 634.3, x=1.5000, pwr= 11 computing order= 1, npoints= 10, at term= 6 f^( 1)(x( 6)= -1.58730158730159E-03 * f(x( 1)) + 1.78571428571429E-02 * f(x( 2) + -9.52380952380952E-02 * f(x( 3) + 3.33333333333333E-01 * f(x( 4) + -1.00000000000000E+00 * f(x( 5) + 2.00000000000000E-01 * f(x( 6) + 6.66666666666667E-01 * f(x( 7) + -1.42857142857143E-01 * f(x( 8) + 2.38095238095238E-02 * f(x( 9) + -1.98412698412698E-03 * f(x( 10) ) + O(h^ 9)f^( 10)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 239.4, x=1.6250, pwr= 8 err= 5.68434188608080E-14, h=0.1250, aderiv= 437.6, x=1.6250, pwr= 9 err= 2.14576722328275E-05, h=0.1250, aderiv= 790.1, x=1.6250, pwr= 10 err= 3.70144844509923E-04, h=0.1250, aderiv=1412.3, x=1.6250, pwr= 11 computing order= 1, npoints= 10, at term= 7 f^( 1)(x( 7)= 1.98412698412698E-03 * f(x( 1)) + -2.14285714285714E-02 * f(x( 2) + 1.07142857142857E-01 * f(x( 3) + -3.33333333333333E-01 * f(x( 4) + 7.50000000000000E-01 * f(x( 5) + -1.50000000000000E+00 * f(x( 6) + 6.16666666666667E-01 * f(x( 7) + 4.28571428571429E-01 * f(x( 8) + -5.35714285714286E-02 * f(x( 9) + 3.96825396825397E-03 * f(x( 10) ) + O(h^ 9)f^( 10)(z) err=-5.68434188608080E-14, h=0.1250, aderiv= 402.1, x=1.7500, pwr= 8 err= 1.13686837721616E-13, h=0.1250, aderiv= 791.7, x=1.7500, pwr= 9 err=-3.21865079513373E-05, h=0.1250, aderiv=1539.4, x=1.7500, pwr= 10 err=-5.59240579150355E-04, h=0.1250, aderiv=2963.3, x=1.7500, pwr= 11 computing order= 1, npoints= 10, at term= 8 f^( 1)(x( 8)= -3.96825396825397E-03 * f(x( 1)) + 4.16666666666667E-02 * f(x( 2) + -2.00000000000000E-01 * f(x( 3) + 5.83333333333333E-01 * f(x( 4) + -1.16666666666667E+00 * f(x( 5) + 1.75000000000000E+00 * f(x( 6) + -2.33333333333333E+00 * f(x( 7) + 1.09285714285714E+00 * f(x( 8) + 2.50000000000000E-01 * f(x( 9) + -1.38888888888889E-02 * f(x( 10) ) + O(h^ 9)f^( 10)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 651.8, x=1.8750, pwr= 8 err= 6.82121026329696E-13, h=0.1250, aderiv=1374.8, x=1.8750, pwr= 9 err= 7.51018533264869E-05, h=0.1250, aderiv=2864.3, x=1.8750, pwr= 10 err= 1.31428241911635E-03, h=0.1250, aderiv=5907.5, x=1.8750, pwr= 11 computing order= 1, npoints= 10, at term= 9 f^( 1)(x( 9)= 1.38888888888889E-02 * f(x( 1)) + -1.42857142857143E-01 * f(x( 2) + 6.66666666666667E-01 * f(x( 3) + -1.86666666666667E+00 * f(x( 4) + 3.50000000000000E+00 * f(x( 5) + -4.66666666666667E+00 * f(x( 6) + 4.66666666666667E+00 * f(x( 7) + -4.00000000000000E+00 * f(x( 8) + 1.71785714285714E+00 * f(x( 9) + 1.11111111111111E-01 * f(x( 10) ) + O(h^ 9)f^( 10)(z) err= 2.27373675443232E-13, h=0.1250, aderiv=1024.0, x=2.0000, pwr= 8 err= 9.09494701772928E-13, h=0.1250, aderiv=2304.0, x=2.0000, pwr= 9 err=-3.00407409667969E-04, h=0.1250, aderiv=5120.0, x=2.0000, pwr= 10 err=-5.29468059721694E-03, h=0.1250, aderiv=11264.0, x=2.0000, pwr= 11 computing order= 1, npoints= 10, at term= 10 f^( 1)(x( 10)= -1.11111111111111E-01 * f(x( 1)) + 1.12500000000000E+00 * f(x( 2) + -5.14285714285714E+00 * f(x( 3) + 1.40000000000000E+01 * f(x( 4) + -2.52000000000000E+01 * f(x( 5) + 3.15000000000000E+01 * f(x( 6) + -2.80000000000000E+01 * f(x( 7) + 1.80000000000000E+01 * f(x( 8) + -9.00000000000000E+00 * f(x( 9) + 2.82896825396825E+00 * f(x( 10) ) + O(h^ 9)f^( 10)(z) err= 0.00000000000000E+00, h=0.1250, aderiv=1565.3, x=2.1250, pwr= 8 err=-3.63797880709171E-12, h=0.1250, aderiv=3742.1, x=2.1250, pwr= 9 err= 2.70366668701172E-03, h=0.1250, aderiv=8835.5, x=2.1250, pwr= 10 err= 4.79900836944580E-02, h=0.1250, aderiv=20652.9, x=2.1250, pwr= 11 computing order= 2, npoints= 3, at term= 1 f^( 2)(x( 1)= 1.00000000000000E+00 * f(x( 1)) + -2.00000000000000E+00 * f(x( 2) + 1.00000000000000E+00 * f(x( 3) ) + O(h^ 1)f^( 3)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 0.0, x=1.0000, pwr= 1 err= 0.00000000000000E+00, h=0.1250, aderiv= 2.0, x=1.0000, pwr= 2 err=-7.50000000000000E-01, h=0.1250, aderiv= 6.0, x=1.0000, pwr= 3 err=-3.21875000000000E+00, h=0.1250, aderiv= 12.0, x=1.0000, pwr= 4 computing order= 2, npoints= 3, at term= 2 f^( 2)(x( 2)= 1.00000000000000E+00 * f(x( 1)) + -2.00000000000000E+00 * f(x( 2) + 1.00000000000000E+00 * f(x( 3) ) + O(h^ 1)f^( 3)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 0.0, x=1.1250, pwr= 1 err= 0.00000000000000E+00, h=0.1250, aderiv= 2.0, x=1.1250, pwr= 2 err= 0.00000000000000E+00, h=0.1250, aderiv= 6.8, x=1.1250, pwr= 3 err=-3.12500000000000E-02, h=0.1250, aderiv= 15.2, x=1.1250, pwr= 4 computing order= 2, npoints= 3, at term= 3 f^( 2)(x( 3)= 1.00000000000000E+00 * f(x( 1)) + -2.00000000000000E+00 * f(x( 2) + 1.00000000000000E+00 * f(x( 3) ) + O(h^ 1)f^( 3)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 0.0, x=1.2500, pwr= 1 err= 0.00000000000000E+00, h=0.1250, aderiv= 2.0, x=1.2500, pwr= 2 err= 7.50000000000000E-01, h=0.1250, aderiv= 7.5, x=1.2500, pwr= 3 err= 3.53125000000000E+00, h=0.1250, aderiv= 18.8, x=1.2500, pwr= 4 computing order= 2, npoints= 4, at term= 1 f^( 2)(x( 1)= 2.00000000000000E+00 * f(x( 1)) + -5.00000000000000E+00 * f(x( 2) + 4.00000000000000E+00 * f(x( 3) + -1.00000000000000E+00 * f(x( 4) ) + O(h^ 2)f^( 4)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 2.0, x=1.0000, pwr= 2 err= 0.00000000000000E+00, h=0.1250, aderiv= 6.0, x=1.0000, pwr= 3 err= 3.43750000000000E-01, h=0.1250, aderiv= 12.0, x=1.0000, pwr= 4 err= 1.95312500000000E+00, h=0.1250, aderiv= 20.0, x=1.0000, pwr= 5 computing order= 2, npoints= 4, at term= 2 f^( 2)(x( 2)= 1.00000000000000E+00 * f(x( 1)) + -2.00000000000000E+00 * f(x( 2) + 1.00000000000000E+00 * f(x( 3) + 0.00000000000000E+00 * f(x( 4) ) + O(h^ 2)f^( 4)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 2.0, x=1.1250, pwr= 2 err= 0.00000000000000E+00, h=0.1250, aderiv= 6.8, x=1.1250, pwr= 3 err=-3.12500000000000E-02, h=0.1250, aderiv= 15.2, x=1.1250, pwr= 4 err=-1.75781250000000E-01, h=0.1250, aderiv= 28.5, x=1.1250, pwr= 5 computing order= 2, npoints= 4, at term= 3 f^( 2)(x( 3)= -0.00000000000000E+00 * f(x( 1)) + 1.00000000000000E+00 * f(x( 2) + -2.00000000000000E+00 * f(x( 3) + 1.00000000000000E+00 * f(x( 4) ) + O(h^ 2)f^( 4)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 2.0, x=1.2500, pwr= 2 err= 0.00000000000000E+00, h=0.1250, aderiv= 7.5, x=1.2500, pwr= 3 err=-3.12500000000000E-02, h=0.1250, aderiv= 18.8, x=1.2500, pwr= 4 err=-1.95312500000000E-01, h=0.1250, aderiv= 39.1, x=1.2500, pwr= 5 computing order= 2, npoints= 4, at term= 4 f^( 2)(x( 4)= -1.00000000000000E+00 * f(x( 1)) + 4.00000000000000E+00 * f(x( 2) + -5.00000000000000E+00 * f(x( 3) + 2.00000000000000E+00 * f(x( 4) ) + O(h^ 2)f^( 4)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 2.0, x=1.3750, pwr= 2 err= 0.00000000000000E+00, h=0.1250, aderiv= 8.3, x=1.3750, pwr= 3 err= 3.43750000000000E-01, h=0.1250, aderiv= 22.7, x=1.3750, pwr= 4 err= 2.12890625000000E+00, h=0.1250, aderiv= 52.0, x=1.3750, pwr= 5 computing order= 2, npoints= 5, at term= 1 f^( 2)(x( 1)= 2.91666666666667E+00 * f(x( 1)) + -8.66666666666667E+00 * f(x( 2) + 9.50000000000000E+00 * f(x( 3) + -4.66666666666667E+00 * f(x( 4) + 9.16666666666667E-01 * f(x( 5) ) + O(h^ 3)f^( 5)(z) err= 1.13686837721616E-13, h=0.1250, aderiv= 6.0, x=1.0000, pwr= 3 err= 0.00000000000000E+00, h=0.1250, aderiv= 12.0, x=1.0000, pwr= 4 err=-1.95312500000000E-01, h=0.1250, aderiv= 20.0, x=1.0000, pwr= 5 err=-1.40429687499977E+00, h=0.1250, aderiv= 30.0, x=1.0000, pwr= 6 computing order= 2, npoints= 5, at term= 2 f^( 2)(x( 2)= 9.16666666666667E-01 * f(x( 1)) + -1.66666666666667E+00 * f(x( 2) + 5.00000000000000E-01 * f(x( 3) + 3.33333333333333E-01 * f(x( 4) + -8.33333333333333E-02 * f(x( 5) ) + O(h^ 3)f^( 5)(z) err= 1.42108547152020E-14, h=0.1250, aderiv= 6.8, x=1.1250, pwr= 3 err= 1.42108547152020E-14, h=0.1250, aderiv= 15.2, x=1.1250, pwr= 4 err= 1.95312500000142E-02, h=0.1250, aderiv= 28.5, x=1.1250, pwr= 5 err= 1.41113281250028E-01, h=0.1250, aderiv= 48.1, x=1.1250, pwr= 6 computing order= 2, npoints= 5, at term= 3 f^( 2)(x( 3)= -8.33333333333333E-02 * f(x( 1)) + 1.33333333333333E+00 * f(x( 2) + -2.50000000000000E+00 * f(x( 3) + 1.33333333333333E+00 * f(x( 4) + -8.33333333333333E-02 * f(x( 5) ) + O(h^ 3)f^( 5)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 7.5, x=1.2500, pwr= 3 err= 5.68434188608080E-14, h=0.1250, aderiv= 18.8, x=1.2500, pwr= 4 err= 5.68434188608080E-14, h=0.1250, aderiv= 39.1, x=1.2500, pwr= 5 err= 1.95312500011369E-03, h=0.1250, aderiv= 73.2, x=1.2500, pwr= 6 computing order= 2, npoints= 5, at term= 4 f^( 2)(x( 4)= -8.33333333333333E-02 * f(x( 1)) + 3.33333333333333E-01 * f(x( 2) + 5.00000000000000E-01 * f(x( 3) + -1.66666666666667E+00 * f(x( 4) + 9.16666666666667E-01 * f(x( 5) ) + O(h^ 3)f^( 5)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 8.3, x=1.3750, pwr= 3 err= 0.00000000000000E+00, h=0.1250, aderiv= 22.7, x=1.3750, pwr= 4 err=-1.95312499998863E-02, h=0.1250, aderiv= 52.0, x=1.3750, pwr= 5 err=-1.51855468749886E-01, h=0.1250, aderiv= 107.2, x=1.3750, pwr= 6 computing order= 2, npoints= 5, at term= 5 f^( 2)(x( 5)= 9.16666666666667E-01 * f(x( 1)) + -4.66666666666667E+00 * f(x( 2) + 9.50000000000000E+00 * f(x( 3) + -8.66666666666667E+00 * f(x( 4) + 2.91666666666667E+00 * f(x( 5) ) + O(h^ 3)f^( 5)(z) err=-2.27373675443232E-13, h=0.1250, aderiv= 9.0, x=1.5000, pwr= 3 err= 0.00000000000000E+00, h=0.1250, aderiv= 27.0, x=1.5000, pwr= 4 err= 1.95312499999773E-01, h=0.1250, aderiv= 67.5, x=1.5000, pwr= 5 err= 1.52539062500000E+00, h=0.1250, aderiv= 151.9, x=1.5000, pwr= 6 computing order= 2, npoints= 6, at term= 1 f^( 2)(x( 1)= 3.75000000000000E+00 * f(x( 1)) + -1.28333333333333E+01 * f(x( 2) + 1.78333333333333E+01 * f(x( 3) + -1.30000000000000E+01 * f(x( 4) + 5.08333333333333E+00 * f(x( 5) + -8.33333333333333E-01 * f(x( 6) ) + O(h^ 4)f^( 6)(z) err= 3.41060513164848E-13, h=0.1250, aderiv= 12.0, x=1.0000, pwr= 4 err= 2.27373675443232E-13, h=0.1250, aderiv= 20.0, x=1.0000, pwr= 5 err= 1.33789062500341E-01, h=0.1250, aderiv= 30.0, x=1.0000, pwr= 6 err= 1.18005371093818E+00, h=0.1250, aderiv= 42.0, x=1.0000, pwr= 7 computing order= 2, npoints= 6, at term= 2 f^( 2)(x( 2)= 8.33333333333333E-01 * f(x( 1)) + -1.25000000000000E+00 * f(x( 2) + -3.33333333333333E-01 * f(x( 3) + 1.16666666666667E+00 * f(x( 4) + -5.00000000000000E-01 * f(x( 5) + 8.33333333333333E-02 * f(x( 6) ) + O(h^ 4)f^( 6)(z) err=-1.42108547152020E-14, h=0.1250, aderiv= 15.2, x=1.1250, pwr= 4 err=-3.55271367880050E-14, h=0.1250, aderiv= 28.5, x=1.1250, pwr= 5 err=-1.26953125000142E-02, h=0.1250, aderiv= 48.1, x=1.1250, pwr= 6 err=-1.12792968750085E-01, h=0.1250, aderiv= 75.7, x=1.1250, pwr= 7 computing order= 2, npoints= 6, at term= 3 f^( 2)(x( 3)= -8.33333333333333E-02 * f(x( 1)) + 1.33333333333333E+00 * f(x( 2) + -2.50000000000000E+00 * f(x( 3) + 1.33333333333333E+00 * f(x( 4) + -8.33333333333333E-02 * f(x( 5) + 0.00000000000000E+00 * f(x( 6) ) + O(h^ 4)f^( 6)(z) err= 5.68434188608080E-14, h=0.1250, aderiv= 18.8, x=1.2500, pwr= 4 err= 5.68434188608080E-14, h=0.1250, aderiv= 39.1, x=1.2500, pwr= 5 err= 1.95312500011369E-03, h=0.1250, aderiv= 73.2, x=1.2500, pwr= 6 err= 1.70898437501137E-02, h=0.1250, aderiv= 128.2, x=1.2500, pwr= 7 computing order= 2, npoints= 6, at term= 4 f^( 2)(x( 4)= -0.00000000000000E+00 * f(x( 1)) + -8.33333333333333E-02 * f(x( 2) + 1.33333333333333E+00 * f(x( 3) + -2.50000000000000E+00 * f(x( 4) + 1.33333333333333E+00 * f(x( 5) + -8.33333333333333E-02 * f(x( 6) ) + O(h^ 4)f^( 6)(z) err= 1.42108547152020E-14, h=0.1250, aderiv= 22.7, x=1.3750, pwr= 4 err= 3.55271367880050E-14, h=0.1250, aderiv= 52.0, x=1.3750, pwr= 5 err= 1.95312500007105E-03, h=0.1250, aderiv= 107.2, x=1.3750, pwr= 6 err= 1.87988281251421E-02, h=0.1250, aderiv= 206.4, x=1.3750, pwr= 7 computing order= 2, npoints= 6, at term= 5 f^( 2)(x( 5)= 8.33333333333333E-02 * f(x( 1)) + -5.00000000000000E-01 * f(x( 2) + 1.16666666666667E+00 * f(x( 3) + -3.33333333333333E-01 * f(x( 4) + -1.25000000000000E+00 * f(x( 5) + 8.33333333333333E-01 * f(x( 6) ) + O(h^ 4)f^( 6)(z) err=-5.68434188608080E-14, h=0.1250, aderiv= 27.0, x=1.5000, pwr= 4 err=-1.13686837721616E-13, h=0.1250, aderiv= 67.5, x=1.5000, pwr= 5 err=-1.26953125001137E-02, h=0.1250, aderiv= 151.9, x=1.5000, pwr= 6 err=-1.20483398437727E-01, h=0.1250, aderiv= 318.9, x=1.5000, pwr= 7 computing order= 2, npoints= 6, at term= 6 f^( 2)(x( 6)= -8.33333333333333E-01 * f(x( 1)) + 5.08333333333333E+00 * f(x( 2) + -1.30000000000000E+01 * f(x( 3) + 1.78333333333333E+01 * f(x( 4) + -1.28333333333333E+01 * f(x( 5) + 3.75000000000000E+00 * f(x( 6) ) + O(h^ 4)f^( 6)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 31.7, x=1.6250, pwr= 4 err= 4.54747350886464E-13, h=0.1250, aderiv= 85.8, x=1.6250, pwr= 5 err= 1.33789062500000E-01, h=0.1250, aderiv= 209.2, x=1.6250, pwr= 6 err= 1.27832031250182E+00, h=0.1250, aderiv= 475.9, x=1.6250, pwr= 7 computing order= 2, npoints= 7, at term= 1 f^( 2)(x( 1)= 4.51111111111111E+00 * f(x( 1)) + -1.74000000000000E+01 * f(x( 2) + 2.92500000000000E+01 * f(x( 3) + -2.82222222222222E+01 * f(x( 4) + 1.65000000000000E+01 * f(x( 5) + -5.40000000000000E+00 * f(x( 6) + 7.61111111111111E-01 * f(x( 7) ) + O(h^ 5)f^( 7)(z) err=-2.27373675443232E-13, h=0.1250, aderiv= 20.0, x=1.0000, pwr= 5 err= 0.00000000000000E+00, h=0.1250, aderiv= 30.0, x=1.0000, pwr= 6 err=-1.07666015625000E-01, h=0.1250, aderiv= 42.0, x=1.0000, pwr= 7 err=-1.13845825195313E+00, h=0.1250, aderiv= 56.0, x=1.0000, pwr= 8 computing order= 2, npoints= 7, at term= 2 f^( 2)(x( 2)= 7.61111111111111E-01 * f(x( 1)) + -8.16666666666667E-01 * f(x( 2) + -1.41666666666667E+00 * f(x( 3) + 2.61111111111111E+00 * f(x( 4) + -1.58333333333333E+00 * f(x( 5) + 5.16666666666667E-01 * f(x( 6) + -7.22222222222222E-02 * f(x( 7) ) + O(h^ 5)f^( 7)(z) err= 1.42108547152020E-14, h=0.1250, aderiv= 28.5, x=1.1250, pwr= 5 err= 2.84217094304040E-14, h=0.1250, aderiv= 48.1, x=1.1250, pwr= 6 err= 9.39941406238631E-03, h=0.1250, aderiv= 75.7, x=1.1250, pwr= 7 err= 1.00128173828125E-01, h=0.1250, aderiv= 113.5, x=1.1250, pwr= 8 computing order= 2, npoints= 7, at term= 3 f^( 2)(x( 3)= -7.22222222222222E-02 * f(x( 1)) + 1.26666666666667E+00 * f(x( 2) + -2.33333333333333E+00 * f(x( 3) + 1.11111111111111E+00 * f(x( 4) + 8.33333333333333E-02 * f(x( 5) + -6.66666666666667E-02 * f(x( 6) + 1.11111111111111E-02 * f(x( 7) ) + O(h^ 5)f^( 7)(z) err= 5.68434188608080E-14, h=0.1250, aderiv= 39.1, x=1.2500, pwr= 5 err= 2.84217094304040E-14, h=0.1250, aderiv= 73.2, x=1.2500, pwr= 6 err=-1.70898437488631E-03, h=0.1250, aderiv= 128.2, x=1.2500, pwr= 7 err=-1.82189941405682E-02, h=0.1250, aderiv= 213.6, x=1.2500, pwr= 8 computing order= 2, npoints= 7, at term= 4 f^( 2)(x( 4)= 1.11111111111111E-02 * f(x( 1)) + -1.50000000000000E-01 * f(x( 2) + 1.50000000000000E+00 * f(x( 3) + -2.72222222222222E+00 * f(x( 4) + 1.50000000000000E+00 * f(x( 5) + -1.50000000000000E-01 * f(x( 6) + 1.11111111111111E-02 * f(x( 7) ) + O(h^ 5)f^( 7)(z) err= 9.23705556488130E-14, h=0.1250, aderiv= 52.0, x=1.3750, pwr= 5 err=-8.52651282912120E-14, h=0.1250, aderiv= 107.2, x=1.3750, pwr= 6 err= 0.00000000000000E+00, h=0.1250, aderiv= 206.4, x=1.3750, pwr= 7 err=-2.74658202897626E-04, h=0.1250, aderiv= 378.4, x=1.3750, pwr= 8 computing order= 2, npoints= 7, at term= 5 f^( 2)(x( 5)= 1.11111111111111E-02 * f(x( 1)) + -6.66666666666667E-02 * f(x( 2) + 8.33333333333333E-02 * f(x( 3) + 1.11111111111111E+00 * f(x( 4) + -2.33333333333333E+00 * f(x( 5) + 1.26666666666667E+00 * f(x( 6) + -7.22222222222222E-02 * f(x( 7) ) + O(h^ 5)f^( 7)(z) err= 7.10542735760100E-14, h=0.1250, aderiv= 67.5, x=1.5000, pwr= 5 err= 1.42108547152020E-13, h=0.1250, aderiv= 151.9, x=1.5000, pwr= 6 err= 1.70898437500000E-03, h=0.1250, aderiv= 318.9, x=1.5000, pwr= 7 err= 1.93786621096024E-02, h=0.1250, aderiv= 637.9, x=1.5000, pwr= 8 computing order= 2, npoints= 7, at term= 6 f^( 2)(x( 6)= -7.22222222222222E-02 * f(x( 1)) + 5.16666666666667E-01 * f(x( 2) + -1.58333333333333E+00 * f(x( 3) + 2.61111111111111E+00 * f(x( 4) + -1.41666666666667E+00 * f(x( 5) + -8.16666666666667E-01 * f(x( 6) + 7.61111111111111E-01 * f(x( 7) ) + O(h^ 5)f^( 7)(z) err=-1.13686837721616E-13, h=0.1250, aderiv= 85.8, x=1.6250, pwr= 5 err= 0.00000000000000E+00, h=0.1250, aderiv= 209.2, x=1.6250, pwr= 6 err=-9.39941406227263E-03, h=0.1250, aderiv= 475.9, x=1.6250, pwr= 7 err=-1.06658935546875E-01, h=0.1250, aderiv=1031.1, x=1.6250, pwr= 8 computing order= 2, npoints= 7, at term= 7 f^( 2)(x( 7)= 7.61111111111111E-01 * f(x( 1)) + -5.40000000000000E+00 * f(x( 2) + 1.65000000000000E+01 * f(x( 3) + -2.82222222222222E+01 * f(x( 4) + 2.92500000000000E+01 * f(x( 5) + -1.74000000000000E+01 * f(x( 6) + 4.51111111111111E+00 * f(x( 7) ) + O(h^ 5)f^( 7)(z) err=-9.09494701772928E-13, h=0.1250, aderiv= 107.2, x=1.7500, pwr= 5 err= 0.00000000000000E+00, h=0.1250, aderiv= 281.4, x=1.7500, pwr= 6 err= 1.07666015623181E-01, h=0.1250, aderiv= 689.3, x=1.7500, pwr= 7 err= 1.23019409179324E+00, h=0.1250, aderiv=1608.5, x=1.7500, pwr= 8 computing order= 2, npoints= 8, at term= 1 f^( 2)(x( 1)= 5.21111111111111E+00 * f(x( 1)) + -2.23000000000000E+01 * f(x( 2) + 4.39500000000000E+01 * f(x( 3) + -5.27222222222222E+01 * f(x( 4) + 4.10000000000000E+01 * f(x( 5) + -2.01000000000000E+01 * f(x( 6) + 5.66111111111111E+00 * f(x( 7) + -7.00000000000000E-01 * f(x( 8) ) + O(h^ 6)f^( 8)(z) err=-2.27373675443232E-13, h=0.1250, aderiv= 30.0, x=1.0000, pwr= 6 err=-4.54747350886464E-13, h=0.1250, aderiv= 42.0, x=1.0000, pwr= 7 err= 9.97009277416510E-02, h=0.1250, aderiv= 56.0, x=1.0000, pwr= 8 err= 1.24145507812500E+00, h=0.1250, aderiv= 72.0, x=1.0000, pwr= 9 computing order= 2, npoints= 8, at term= 2 f^( 2)(x( 2)= 7.00000000000000E-01 * f(x( 1)) + -3.88888888888889E-01 * f(x( 2) + -2.70000000000000E+00 * f(x( 3) + 4.75000000000000E+00 * f(x( 4) + -3.72222222222222E+00 * f(x( 5) + 1.80000000000000E+00 * f(x( 6) + -5.00000000000000E-01 * f(x( 7) + 6.11111111111111E-02 * f(x( 8) ) + O(h^ 6)f^( 8)(z) err=-2.27373675443232E-13, h=0.1250, aderiv= 48.1, x=1.1250, pwr= 6 err=-4.54747350886464E-13, h=0.1250, aderiv= 75.7, x=1.1250, pwr= 7 err=-7.96508789062500E-03, h=0.1250, aderiv= 113.5, x=1.1250, pwr= 8 err=-9.98725891122376E-02, h=0.1250, aderiv= 164.2, x=1.1250, pwr= 9 computing order= 2, npoints= 8, at term= 3 f^( 2)(x( 3)= -6.11111111111111E-02 * f(x( 1)) + 1.18888888888889E+00 * f(x( 2) + -2.10000000000000E+00 * f(x( 3) + 7.22222222222222E-01 * f(x( 4) + 4.72222222222222E-01 * f(x( 5) + -3.00000000000000E-01 * f(x( 6) + 8.88888888888889E-02 * f(x( 7) + -1.11111111111111E-02 * f(x( 8) ) + O(h^ 6)f^( 8)(z) err=-5.68434188608080E-14, h=0.1250, aderiv= 73.2, x=1.2500, pwr= 6 err=-1.13686837721616E-13, h=0.1250, aderiv= 128.2, x=1.2500, pwr= 7 err= 1.43432617176131E-03, h=0.1250, aderiv= 213.6, x=1.2500, pwr= 8 err= 1.80587768555824E-02, h=0.1250, aderiv= 343.3, x=1.2500, pwr= 9 computing order= 2, npoints= 8, at term= 4 f^( 2)(x( 4)= 1.11111111111111E-02 * f(x( 1)) + -1.50000000000000E-01 * f(x( 2) + 1.50000000000000E+00 * f(x( 3) + -2.72222222222222E+00 * f(x( 4) + 1.50000000000000E+00 * f(x( 5) + -1.50000000000000E-01 * f(x( 6) + 1.11111111111111E-02 * f(x( 7) + 0.00000000000000E+00 * f(x( 8) ) + O(h^ 6)f^( 8)(z) err=-8.52651282912120E-14, h=0.1250, aderiv= 107.2, x=1.3750, pwr= 6 err= 0.00000000000000E+00, h=0.1250, aderiv= 206.4, x=1.3750, pwr= 7 err=-2.74658202897626E-04, h=0.1250, aderiv= 378.4, x=1.3750, pwr= 8 err=-3.39889526367188E-03, h=0.1250, aderiv= 669.0, x=1.3750, pwr= 9 computing order= 2, npoints= 8, at term= 5 f^( 2)(x( 5)= -0.00000000000000E+00 * f(x( 1)) + 1.11111111111111E-02 * f(x( 2) + -1.50000000000000E-01 * f(x( 3) + 1.50000000000000E+00 * f(x( 4) + -2.72222222222222E+00 * f(x( 5) + 1.50000000000000E+00 * f(x( 6) + -1.50000000000000E-01 * f(x( 7) + 1.11111111111111E-02 * f(x( 8) ) + O(h^ 6)f^( 8)(z) err=-1.13686837721616E-13, h=0.1250, aderiv= 151.9, x=1.5000, pwr= 6 err= 1.70530256582424E-13, h=0.1250, aderiv= 318.9, x=1.5000, pwr= 7 err=-2.74658203011313E-04, h=0.1250, aderiv= 637.9, x=1.5000, pwr= 8 err=-3.70788574241487E-03, h=0.1250, aderiv=1230.2, x=1.5000, pwr= 9 computing order= 2, npoints= 8, at term= 6 f^( 2)(x( 6)= -1.11111111111111E-02 * f(x( 1)) + 8.88888888888889E-02 * f(x( 2) + -3.00000000000000E-01 * f(x( 3) + 4.72222222222222E-01 * f(x( 4) + 7.22222222222222E-01 * f(x( 5) + -2.10000000000000E+00 * f(x( 6) + 1.18888888888889E+00 * f(x( 7) + -6.11111111111111E-02 * f(x( 8) ) + O(h^ 6)f^( 8)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 209.2, x=1.6250, pwr= 6 err= 0.00000000000000E+00, h=0.1250, aderiv= 475.9, x=1.6250, pwr= 7 err= 1.43432617278449E-03, h=0.1250, aderiv=1031.1, x=1.6250, pwr= 8 err= 1.90544128417969E-02, h=0.1250, aderiv=2154.3, x=1.6250, pwr= 9 computing order= 2, npoints= 8, at term= 7 f^( 2)(x( 7)= 6.11111111111111E-02 * f(x( 1)) + -5.00000000000000E-01 * f(x( 2) + 1.80000000000000E+00 * f(x( 3) + -3.72222222222222E+00 * f(x( 4) + 4.75000000000000E+00 * f(x( 5) + -2.70000000000000E+00 * f(x( 6) + -3.88888888888889E-01 * f(x( 7) + 7.00000000000000E-01 * f(x( 8) ) + O(h^ 6)f^( 8)(z) err= 4.54747350886464E-13, h=0.1250, aderiv= 281.4, x=1.7500, pwr= 6 err= 9.09494701772928E-13, h=0.1250, aderiv= 689.3, x=1.7500, pwr= 7 err=-7.96508788880601E-03, h=0.1250, aderiv=1608.5, x=1.7500, pwr= 8 err=-1.06224060056775E-01, h=0.1250, aderiv=3619.1, x=1.7500, pwr= 9 computing order= 2, npoints= 8, at term= 8 f^( 2)(x( 8)= -7.00000000000000E-01 * f(x( 1)) + 5.66111111111111E+00 * f(x( 2) + -2.01000000000000E+01 * f(x( 3) + 4.10000000000000E+01 * f(x( 4) + -5.27222222222222E+01 * f(x( 5) + 4.39500000000000E+01 * f(x( 6) + -2.23000000000000E+01 * f(x( 7) + 5.21111111111111E+00 * f(x( 8) ) + O(h^ 6)f^( 8)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 370.8, x=1.8750, pwr= 6 err= 7.27595761418343E-12, h=0.1250, aderiv= 973.3, x=1.8750, pwr= 7 err= 9.97009277198231E-02, h=0.1250, aderiv=2433.3, x=1.8750, pwr= 8 err= 1.33830642698740E+00, h=0.1250, aderiv=5866.0, x=1.8750, pwr= 9 computing order= 3, npoints= 4, at term= 1 f^( 3)(x( 1)= -1.00000000000000E+00 * f(x( 1)) + 3.00000000000000E+00 * f(x( 2) + -3.00000000000000E+00 * f(x( 3) + 1.00000000000000E+00 * f(x( 4) ) + O(h^ 1)f^( 4)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 0.0, x=1.0000, pwr= 2 err= 0.00000000000000E+00, h=0.1250, aderiv= 6.0, x=1.0000, pwr= 3 err=-4.50000000000000E+00, h=0.1250, aderiv= 24.0, x=1.0000, pwr= 4 err=-2.48437500000000E+01, h=0.1250, aderiv= 60.0, x=1.0000, pwr= 5 computing order= 3, npoints= 4, at term= 2 f^( 3)(x( 2)= -1.00000000000000E+00 * f(x( 1)) + 3.00000000000000E+00 * f(x( 2) + -3.00000000000000E+00 * f(x( 3) + 1.00000000000000E+00 * f(x( 4) ) + O(h^ 1)f^( 4)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 0.0, x=1.1250, pwr= 2 err= 0.00000000000000E+00, h=0.1250, aderiv= 6.0, x=1.1250, pwr= 3 err=-1.50000000000000E+00, h=0.1250, aderiv= 27.0, x=1.1250, pwr= 4 err=-8.90625000000000E+00, h=0.1250, aderiv= 75.9, x=1.1250, pwr= 5 computing order= 3, npoints= 4, at term= 3 f^( 3)(x( 3)= -1.00000000000000E+00 * f(x( 1)) + 3.00000000000000E+00 * f(x( 2) + -3.00000000000000E+00 * f(x( 3) + 1.00000000000000E+00 * f(x( 4) ) + O(h^ 1)f^( 4)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 0.0, x=1.2500, pwr= 2 err= 0.00000000000000E+00, h=0.1250, aderiv= 6.0, x=1.2500, pwr= 3 err= 1.50000000000000E+00, h=0.1250, aderiv= 30.0, x=1.2500, pwr= 4 err= 8.90625000000000E+00, h=0.1250, aderiv= 93.8, x=1.2500, pwr= 5 computing order= 3, npoints= 4, at term= 4 f^( 3)(x( 4)= -1.00000000000000E+00 * f(x( 1)) + 3.00000000000000E+00 * f(x( 2) + -3.00000000000000E+00 * f(x( 3) + 1.00000000000000E+00 * f(x( 4) ) + O(h^ 1)f^( 4)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 0.0, x=1.3750, pwr= 2 err= 0.00000000000000E+00, h=0.1250, aderiv= 6.0, x=1.3750, pwr= 3 err= 4.50000000000000E+00, h=0.1250, aderiv= 33.0, x=1.3750, pwr= 4 err= 2.85937500000000E+01, h=0.1250, aderiv= 113.4, x=1.3750, pwr= 5 computing order= 3, npoints= 5, at term= 1 f^( 3)(x( 1)= -2.50000000000000E+00 * f(x( 1)) + 9.00000000000000E+00 * f(x( 2) + -1.20000000000000E+01 * f(x( 3) + 7.00000000000000E+00 * f(x( 4) + -1.50000000000000E+00 * f(x( 5) ) + O(h^ 2)f^( 5)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 6.0, x=1.0000, pwr= 3 err= 0.00000000000000E+00, h=0.1250, aderiv= 24.0, x=1.0000, pwr= 4 err= 3.28125000000000E+00, h=0.1250, aderiv= 60.0, x=1.0000, pwr= 5 err= 2.32031250000000E+01, h=0.1250, aderiv= 120.0, x=1.0000, pwr= 6 computing order= 3, npoints= 5, at term= 2 f^( 3)(x( 2)= -1.50000000000000E+00 * f(x( 1)) + 5.00000000000000E+00 * f(x( 2) + -6.00000000000000E+00 * f(x( 3) + 3.00000000000000E+00 * f(x( 4) + -5.00000000000000E-01 * f(x( 5) ) + O(h^ 2)f^( 5)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 6.0, x=1.1250, pwr= 3 err= 0.00000000000000E+00, h=0.1250, aderiv= 27.0, x=1.1250, pwr= 4 err= 4.68750000000000E-01, h=0.1250, aderiv= 75.9, x=1.1250, pwr= 5 err= 3.51562500000000E+00, h=0.1250, aderiv= 170.9, x=1.1250, pwr= 6 computing order= 3, npoints= 5, at term= 3 f^( 3)(x( 3)= -5.00000000000000E-01 * f(x( 1)) + 1.00000000000000E+00 * f(x( 2) + 0.00000000000000E+00 * f(x( 3) + -1.00000000000000E+00 * f(x( 4) + 5.00000000000000E-01 * f(x( 5) ) + O(h^ 2)f^( 5)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 6.0, x=1.2500, pwr= 3 err= 0.00000000000000E+00, h=0.1250, aderiv= 30.0, x=1.2500, pwr= 4 err=-4.68750000000000E-01, h=0.1250, aderiv= 93.8, x=1.2500, pwr= 5 err=-3.51562500000000E+00, h=0.1250, aderiv= 234.4, x=1.2500, pwr= 6 computing order= 3, npoints= 5, at term= 4 f^( 3)(x( 4)= 5.00000000000000E-01 * f(x( 1)) + -3.00000000000000E+00 * f(x( 2) + 6.00000000000000E+00 * f(x( 3) + -5.00000000000000E+00 * f(x( 4) + 1.50000000000000E+00 * f(x( 5) ) + O(h^ 2)f^( 5)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 6.0, x=1.3750, pwr= 3 err= 0.00000000000000E+00, h=0.1250, aderiv= 33.0, x=1.3750, pwr= 4 err= 4.68750000000000E-01, h=0.1250, aderiv= 113.4, x=1.3750, pwr= 5 err= 3.51562500000000E+00, h=0.1250, aderiv= 312.0, x=1.3750, pwr= 6 computing order= 3, npoints= 5, at term= 5 f^( 3)(x( 5)= 1.50000000000000E+00 * f(x( 1)) + -7.00000000000000E+00 * f(x( 2) + 1.20000000000000E+01 * f(x( 3) + -9.00000000000000E+00 * f(x( 4) + 2.50000000000000E+00 * f(x( 5) ) + O(h^ 2)f^( 5)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 6.0, x=1.5000, pwr= 3 err= 0.00000000000000E+00, h=0.1250, aderiv= 36.0, x=1.5000, pwr= 4 err= 3.28125000000000E+00, h=0.1250, aderiv= 135.0, x=1.5000, pwr= 5 err= 2.60156250000000E+01, h=0.1250, aderiv= 405.0, x=1.5000, pwr= 6 computing order= 3, npoints= 6, at term= 1 f^( 3)(x( 1)= -4.25000000000000E+00 * f(x( 1)) + 1.77500000000000E+01 * f(x( 2) + -2.95000000000000E+01 * f(x( 3) + 2.45000000000000E+01 * f(x( 4) + -1.02500000000000E+01 * f(x( 5) + 1.75000000000000E+00 * f(x( 6) ) + O(h^ 3)f^( 6)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 24.0, x=1.0000, pwr= 4 err= 0.00000000000000E+00, h=0.1250, aderiv= 60.0, x=1.0000, pwr= 5 err=-2.63671875000000E+00, h=0.1250, aderiv= 120.0, x=1.0000, pwr= 6 err=-2.29995117187500E+01, h=0.1250, aderiv= 210.0, x=1.0000, pwr= 7 computing order= 3, npoints= 6, at term= 2 f^( 3)(x( 2)= -1.75000000000000E+00 * f(x( 1)) + 6.25000000000000E+00 * f(x( 2) + -8.50000000000000E+00 * f(x( 3) + 5.50000000000000E+00 * f(x( 4) + -1.75000000000000E+00 * f(x( 5) + 2.50000000000000E-01 * f(x( 6) ) + O(h^ 3)f^( 6)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 27.0, x=1.1250, pwr= 4 err= 0.00000000000000E+00, h=0.1250, aderiv= 75.9, x=1.1250, pwr= 5 err=-1.75781250000000E-01, h=0.1250, aderiv= 170.9, x=1.1250, pwr= 6 err=-1.62011718750000E+00, h=0.1250, aderiv= 336.4, x=1.1250, pwr= 7 computing order= 3, npoints= 6, at term= 3 f^( 3)(x( 3)= -2.50000000000000E-01 * f(x( 1)) + -2.50000000000000E-01 * f(x( 2) + 2.50000000000000E+00 * f(x( 3) + -3.50000000000000E+00 * f(x( 4) + 1.75000000000000E+00 * f(x( 5) + -2.50000000000000E-01 * f(x( 6) ) + O(h^ 3)f^( 6)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 30.0, x=1.2500, pwr= 4 err= 0.00000000000000E+00, h=0.1250, aderiv= 93.8, x=1.2500, pwr= 5 err= 1.75781250000000E-01, h=0.1250, aderiv= 234.4, x=1.2500, pwr= 6 err= 1.60986328125000E+00, h=0.1250, aderiv= 512.7, x=1.2500, pwr= 7 computing order= 3, npoints= 6, at term= 4 f^( 3)(x( 4)= 2.50000000000000E-01 * f(x( 1)) + -1.75000000000000E+00 * f(x( 2) + 3.50000000000000E+00 * f(x( 3) + -2.50000000000000E+00 * f(x( 4) + 2.50000000000000E-01 * f(x( 5) + 2.50000000000000E-01 * f(x( 6) ) + O(h^ 3)f^( 6)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 33.0, x=1.3750, pwr= 4 err= 0.00000000000000E+00, h=0.1250, aderiv= 113.4, x=1.3750, pwr= 5 err=-1.75781250000000E-01, h=0.1250, aderiv= 312.0, x=1.3750, pwr= 6 err=-1.62011718750000E+00, h=0.1250, aderiv= 750.6, x=1.3750, pwr= 7 computing order= 3, npoints= 6, at term= 5 f^( 3)(x( 5)= -2.50000000000000E-01 * f(x( 1)) + 1.75000000000000E+00 * f(x( 2) + -5.50000000000000E+00 * f(x( 3) + 8.50000000000000E+00 * f(x( 4) + -6.25000000000000E+00 * f(x( 5) + 1.75000000000000E+00 * f(x( 6) ) + O(h^ 3)f^( 6)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 36.0, x=1.5000, pwr= 4 err= 0.00000000000000E+00, h=0.1250, aderiv= 135.0, x=1.5000, pwr= 5 err= 1.75781250000000E-01, h=0.1250, aderiv= 405.0, x=1.5000, pwr= 6 err= 1.60986328125000E+00, h=0.1250, aderiv=1063.1, x=1.5000, pwr= 7 computing order= 3, npoints= 6, at term= 6 f^( 3)(x( 6)= -1.75000000000000E+00 * f(x( 1)) + 1.02500000000000E+01 * f(x( 2) + -2.45000000000000E+01 * f(x( 3) + 2.95000000000000E+01 * f(x( 4) + -1.77500000000000E+01 * f(x( 5) + 4.25000000000000E+00 * f(x( 6) ) + O(h^ 3)f^( 6)(z) err= 0.00000000000000E+00, h=0.1250, aderiv= 39.0, x=1.6250, pwr= 4 err= 0.00000000000000E+00, h=0.1250, aderiv= 158.4, x=1.6250, pwr= 5 err= 2.63671875000000E+00, h=0.1250, aderiv= 514.9, x=1.6250, pwr= 6 err= 2.54501953125000E+01, h=0.1250, aderiv=1464.3, x=1.6250, pwr= 7 computing order= 4, npoints= 5, at term= 1 f^( 4)(x( 1)= 1.00000000000000E+00 * f(x( 1)) + -4.00000000000000E+00 * f(x( 2) + 6.00000000000000E+00 * f(x( 3) + -4.00000000000000E+00 * f(x( 4) + 1.00000000000000E+00 * f(x( 5) ) + O(h^ 1)f^( 5)(z) err= 0.00000000000000E+00, h=0.0625, aderiv= 0.0, x=1.0000, pwr= 3 err= 0.00000000000000E+00, h=0.0625, aderiv= 24.0, x=1.0000, pwr= 4 err=-1.50000000000000E+01, h=0.0625, aderiv= 120.0, x=1.0000, pwr= 5 err=-9.60937500000000E+01, h=0.0625, aderiv= 360.0, x=1.0000, pwr= 6 computing order= 4, npoints= 5, at term= 2 f^( 4)(x( 2)= 1.00000000000000E+00 * f(x( 1)) + -4.00000000000000E+00 * f(x( 2) + 6.00000000000000E+00 * f(x( 3) + -4.00000000000000E+00 * f(x( 4) + 1.00000000000000E+00 * f(x( 5) ) + O(h^ 1)f^( 5)(z) err= 0.00000000000000E+00, h=0.0625, aderiv= 0.0, x=1.0625, pwr= 3 err= 0.00000000000000E+00, h=0.0625, aderiv= 24.0, x=1.0625, pwr= 4 err=-7.50000000000000E+00, h=0.0625, aderiv= 127.5, x=1.0625, pwr= 5 err=-4.96875000000000E+01, h=0.0625, aderiv= 406.4, x=1.0625, pwr= 6 computing order= 4, npoints= 5, at term= 3 f^( 4)(x( 3)= 1.00000000000000E+00 * f(x( 1)) + -4.00000000000000E+00 * f(x( 2) + 6.00000000000000E+00 * f(x( 3) + -4.00000000000000E+00 * f(x( 4) + 1.00000000000000E+00 * f(x( 5) ) + O(h^ 1)f^( 5)(z) err= 0.00000000000000E+00, h=0.0625, aderiv= 0.0, x=1.1250, pwr= 3 err= 0.00000000000000E+00, h=0.0625, aderiv= 24.0, x=1.1250, pwr= 4 err= 0.00000000000000E+00, h=0.0625, aderiv= 135.0, x=1.1250, pwr= 5 err=-4.68750000000000E-01, h=0.0625, aderiv= 455.6, x=1.1250, pwr= 6 computing order= 4, npoints= 5, at term= 4 f^( 4)(x( 4)= 1.00000000000000E+00 * f(x( 1)) + -4.00000000000000E+00 * f(x( 2) + 6.00000000000000E+00 * f(x( 3) + -4.00000000000000E+00 * f(x( 4) + 1.00000000000000E+00 * f(x( 5) ) + O(h^ 1)f^( 5)(z) err= 0.00000000000000E+00, h=0.0625, aderiv= 0.0, x=1.1875, pwr= 3 err= 0.00000000000000E+00, h=0.0625, aderiv= 24.0, x=1.1875, pwr= 4 err= 7.50000000000000E+00, h=0.0625, aderiv= 142.5, x=1.1875, pwr= 5 err= 5.15625000000000E+01, h=0.0625, aderiv= 507.7, x=1.1875, pwr= 6 computing order= 4, npoints= 5, at term= 5 f^( 4)(x( 5)= 1.00000000000000E+00 * f(x( 1)) + -4.00000000000000E+00 * f(x( 2) + 6.00000000000000E+00 * f(x( 3) + -4.00000000000000E+00 * f(x( 4) + 1.00000000000000E+00 * f(x( 5) ) + O(h^ 1)f^( 5)(z) err= 0.00000000000000E+00, h=0.0625, aderiv= 0.0, x=1.2500, pwr= 3 err= 0.00000000000000E+00, h=0.0625, aderiv= 24.0, x=1.2500, pwr= 4 err= 1.50000000000000E+01, h=0.0625, aderiv= 150.0, x=1.2500, pwr= 5 err= 1.06406250000000E+02, h=0.0625, aderiv= 562.5, x=1.2500, pwr= 6 computing order= 4, npoints= 6, at term= 1 f^( 4)(x( 1)= 3.00000000000000E+00 * f(x( 1)) + -1.40000000000000E+01 * f(x( 2) + 2.60000000000000E+01 * f(x( 3) + -2.40000000000000E+01 * f(x( 4) + 1.10000000000000E+01 * f(x( 5) + -2.00000000000000E+00 * f(x( 6) ) + O(h^ 2)f^( 6)(z) err= 0.00000000000000E+00, h=0.0625, aderiv= 24.0, x=1.0000, pwr= 4 err= 0.00000000000000E+00, h=0.0625, aderiv= 120.0, x=1.0000, pwr= 5 err= 7.96875000000000E+00, h=0.0625, aderiv= 360.0, x=1.0000, pwr= 6 err= 6.19335937500000E+01, h=0.0625, aderiv= 840.0, x=1.0000, pwr= 7 computing order= 4, npoints= 6, at term= 2 f^( 4)(x( 2)= 2.00000000000000E+00 * f(x( 1)) + -9.00000000000000E+00 * f(x( 2) + 1.60000000000000E+01 * f(x( 3) + -1.40000000000000E+01 * f(x( 4) + 6.00000000000000E+00 * f(x( 5) + -1.00000000000000E+00 * f(x( 6) ) + O(h^ 2)f^( 6)(z) err= 0.00000000000000E+00, h=0.0625, aderiv= 24.0, x=1.0625, pwr= 4 err= 0.00000000000000E+00, h=0.0625, aderiv= 127.5, x=1.0625, pwr= 5 err= 2.34375000000000E+00, h=0.0625, aderiv= 406.4, x=1.0625, pwr= 6 err= 1.86621093750000E+01, h=0.0625, aderiv=1007.5, x=1.0625, pwr= 7 computing order= 4, npoints= 6, at term= 3 f^( 4)(x( 3)= 1.00000000000000E+00 * f(x( 1)) + -4.00000000000000E+00 * f(x( 2) + 6.00000000000000E+00 * f(x( 3) + -4.00000000000000E+00 * f(x( 4) + 1.00000000000000E+00 * f(x( 5) + 0.00000000000000E+00 * f(x( 6) ) + O(h^ 2)f^( 6)(z) err= 0.00000000000000E+00, h=0.0625, aderiv= 24.0, x=1.1250, pwr= 4 err= 0.00000000000000E+00, h=0.0625, aderiv= 135.0, x=1.1250, pwr= 5 err=-4.68750000000000E-01, h=0.0625, aderiv= 455.6, x=1.1250, pwr= 6 err=-3.69140625000000E+00, h=0.0625, aderiv=1196.0, x=1.1250, pwr= 7 computing order= 4, npoints= 6, at term= 4 f^( 4)(x( 4)= -0.00000000000000E+00 * f(x( 1)) + 1.00000000000000E+00 * f(x( 2) + -4.00000000000000E+00 * f(x( 3) + 6.00000000000000E+00 * f(x( 4) + -4.00000000000000E+00 * f(x( 5) + 1.00000000000000E+00 * f(x( 6) ) + O(h^ 2)f^( 6)(z) err= 0.00000000000000E+00, h=0.0625, aderiv= 24.0, x=1.1875, pwr= 4 err= 0.00000000000000E+00, h=0.0625, aderiv= 142.5, x=1.1875, pwr= 5 err=-4.68750000000000E-01, h=0.0625, aderiv= 507.7, x=1.1875, pwr= 6 err=-3.89648437500000E+00, h=0.0625, aderiv=1406.6, x=1.1875, pwr= 7 computing order= 4, npoints= 6, at term= 5 f^( 4)(x( 5)= -1.00000000000000E+00 * f(x( 1)) + 6.00000000000000E+00 * f(x( 2) + -1.40000000000000E+01 * f(x( 3) + 1.60000000000000E+01 * f(x( 4) + -9.00000000000000E+00 * f(x( 5) + 2.00000000000000E+00 * f(x( 6) ) + O(h^ 2)f^( 6)(z) err= 0.00000000000000E+00, h=0.0625, aderiv= 24.0, x=1.2500, pwr= 4 err= 0.00000000000000E+00, h=0.0625, aderiv= 150.0, x=1.2500, pwr= 5 err= 2.34375000000000E+00, h=0.0625, aderiv= 562.5, x=1.2500, pwr= 6 err= 1.92773437500000E+01, h=0.0625, aderiv=1640.6, x=1.2500, pwr= 7 computing order= 4, npoints= 6, at term= 6 f^( 4)(x( 6)= -2.00000000000000E+00 * f(x( 1)) + 1.10000000000000E+01 * f(x( 2) + -2.40000000000000E+01 * f(x( 3) + 2.60000000000000E+01 * f(x( 4) + -1.40000000000000E+01 * f(x( 5) + 3.00000000000000E+00 * f(x( 6) ) + O(h^ 2)f^( 6)(z) err= 0.00000000000000E+00, h=0.0625, aderiv= 24.0, x=1.3125, pwr= 4 err= 0.00000000000000E+00, h=0.0625, aderiv= 157.5, x=1.3125, pwr= 5 err= 7.96875000000000E+00, h=0.0625, aderiv= 620.2, x=1.3125, pwr= 6 err= 6.70605468750000E+01, h=0.0625, aderiv=1899.2, x=1.3125, pwr= 7 computing order= 5, npoints= 6, at term= 1 f^( 5)(x( 1)= -1.00000000000000E+00 * f(x( 1)) + 5.00000000000000E+00 * f(x( 2) + -1.00000000000000E+01 * f(x( 3) + 1.00000000000000E+01 * f(x( 4) + -5.00000000000000E+00 * f(x( 5) + 1.00000000000000E+00 * f(x( 6) ) + O(h^ 1)f^( 6)(z) err= 0.00000000000000E+00, h=0.0313, aderiv= 0.0, x=1.0000, pwr= 4 err= 0.00000000000000E+00, h=0.0313, aderiv= 120.0, x=1.0000, pwr= 5 err=-5.62500000000000E+01, h=0.0313, aderiv= 720.0, x=1.0000, pwr= 6 err=-4.10156250000000E+02, h=0.0313, aderiv=2520.0, x=1.0000, pwr= 7 computing order= 5, npoints= 6, at term= 2 f^( 5)(x( 2)= -1.00000000000000E+00 * f(x( 1)) + 5.00000000000000E+00 * f(x( 2) + -1.00000000000000E+01 * f(x( 3) + 1.00000000000000E+01 * f(x( 4) + -5.00000000000000E+00 * f(x( 5) + 1.00000000000000E+00 * f(x( 6) ) + O(h^ 1)f^( 6)(z) err= 0.00000000000000E+00, h=0.0313, aderiv= 0.0, x=1.0313, pwr= 4 err= 0.00000000000000E+00, h=0.0313, aderiv= 120.0, x=1.0313, pwr= 5 err=-3.37500000000000E+01, h=0.0313, aderiv= 742.5, x=1.0313, pwr= 6 err=-2.50195312500000E+02, h=0.0313, aderiv=2680.0, x=1.0313, pwr= 7 computing order= 5, npoints= 6, at term= 3 f^( 5)(x( 3)= -1.00000000000000E+00 * f(x( 1)) + 5.00000000000000E+00 * f(x( 2) + -1.00000000000000E+01 * f(x( 3) + 1.00000000000000E+01 * f(x( 4) + -5.00000000000000E+00 * f(x( 5) + 1.00000000000000E+00 * f(x( 6) ) + O(h^ 1)f^( 6)(z) err= 0.00000000000000E+00, h=0.0313, aderiv= 0.0, x=1.0625, pwr= 4 err= 0.00000000000000E+00, h=0.0313, aderiv= 120.0, x=1.0625, pwr= 5 err=-1.12500000000000E+01, h=0.0313, aderiv= 765.0, x=1.0625, pwr= 6 err=-8.53125000000000E+01, h=0.0313, aderiv=2844.8, x=1.0625, pwr= 7 computing order= 5, npoints= 6, at term= 4 f^( 5)(x( 4)= -1.00000000000000E+00 * f(x( 1)) + 5.00000000000000E+00 * f(x( 2) + -1.00000000000000E+01 * f(x( 3) + 1.00000000000000E+01 * f(x( 4) + -5.00000000000000E+00 * f(x( 5) + 1.00000000000000E+00 * f(x( 6) ) + O(h^ 1)f^( 6)(z) err= 0.00000000000000E+00, h=0.0313, aderiv= 0.0, x=1.0938, pwr= 4 err= 0.00000000000000E+00, h=0.0313, aderiv= 120.0, x=1.0938, pwr= 5 err= 1.12500000000000E+01, h=0.0313, aderiv= 787.5, x=1.0938, pwr= 6 err= 8.44921875000000E+01, h=0.0313, aderiv=3014.6, x=1.0938, pwr= 7 computing order= 5, npoints= 6, at term= 5 f^( 5)(x( 5)= -1.00000000000000E+00 * f(x( 1)) + 5.00000000000000E+00 * f(x( 2) + -1.00000000000000E+01 * f(x( 3) + 1.00000000000000E+01 * f(x( 4) + -5.00000000000000E+00 * f(x( 5) + 1.00000000000000E+00 * f(x( 6) ) + O(h^ 1)f^( 6)(z) err= 0.00000000000000E+00, h=0.0313, aderiv= 0.0, x=1.1250, pwr= 4 err= 0.00000000000000E+00, h=0.0313, aderiv= 120.0, x=1.1250, pwr= 5 err= 3.37500000000000E+01, h=0.0313, aderiv= 810.0, x=1.1250, pwr= 6 err= 2.59218750000000E+02, h=0.0313, aderiv=3189.4, x=1.1250, pwr= 7 computing order= 5, npoints= 6, at term= 6 f^( 5)(x( 6)= -1.00000000000000E+00 * f(x( 1)) + 5.00000000000000E+00 * f(x( 2) + -1.00000000000000E+01 * f(x( 3) + 1.00000000000000E+01 * f(x( 4) + -5.00000000000000E+00 * f(x( 5) + 1.00000000000000E+00 * f(x( 6) ) + O(h^ 1)f^( 6)(z) err= 0.00000000000000E+00, h=0.0313, aderiv= 0.0, x=1.1563, pwr= 4 err= 0.00000000000000E+00, h=0.0313, aderiv= 120.0, x=1.1563, pwr= 5 err= 5.62500000000000E+01, h=0.0313, aderiv= 832.5, x=1.1563, pwr= 6 err= 4.38867187500000E+02, h=0.0313, aderiv=3369.0, x=1.1563, pwr= 7 computing order= 5, npoints= 7, at term= 1 f^( 5)(x( 1)= -3.50000000000000E+00 * f(x( 1)) + 2.00000000000000E+01 * f(x( 2) + -4.75000000000000E+01 * f(x( 3) + 6.00000000000000E+01 * f(x( 4) + -4.25000000000000E+01 * f(x( 5) + 1.60000000000000E+01 * f(x( 6) + -2.50000000000000E+00 * f(x( 7) ) + O(h^ 2)f^( 7)(z) err= 0.00000000000000E+00, h=0.0313, aderiv= 120.0, x=1.0000, pwr= 5 err= 0.00000000000000E+00, h=0.0313, aderiv= 720.0, x=1.0000, pwr= 6 err= 2.05078125000000E+01, h=0.0313, aderiv=2520.0, x=1.0000, pwr= 7 err= 1.74829101562500E+02, h=0.0313, aderiv=6720.0, x=1.0000, pwr= 8 computing order= 5, npoints= 7, at term= 2 f^( 5)(x( 2)= -2.50000000000000E+00 * f(x( 1)) + 1.40000000000000E+01 * f(x( 2) + -3.25000000000000E+01 * f(x( 3) + 4.00000000000000E+01 * f(x( 4) + -2.75000000000000E+01 * f(x( 5) + 1.00000000000000E+01 * f(x( 6) + -1.50000000000000E+00 * f(x( 7) ) + O(h^ 2)f^( 7)(z) err= 0.00000000000000E+00, h=0.0313, aderiv= 120.0, x=1.0313, pwr= 5 err= 0.00000000000000E+00, h=0.0313, aderiv= 742.5, x=1.0313, pwr= 6 err= 8.20312500000000E+00, h=0.0313, aderiv=2680.0, x=1.0313, pwr= 7 err= 7.07519531250000E+01, h=0.0313, aderiv=7369.9, x=1.0313, pwr= 8 computing order= 5, npoints= 7, at term= 3 f^( 5)(x( 3)= -1.50000000000000E+00 * f(x( 1)) + 8.00000000000000E+00 * f(x( 2) + -1.75000000000000E+01 * f(x( 3) + 2.00000000000000E+01 * f(x( 4) + -1.25000000000000E+01 * f(x( 5) + 4.00000000000000E+00 * f(x( 6) + -5.00000000000000E-01 * f(x( 7) ) + O(h^ 2)f^( 7)(z) err= 0.00000000000000E+00, h=0.0313, aderiv= 120.0, x=1.0625, pwr= 5 err= 0.00000000000000E+00, h=0.0313, aderiv= 765.0, x=1.0625, pwr= 6 err= 8.20312500000000E-01, h=0.0313, aderiv=2844.8, x=1.0625, pwr= 7 err= 7.28027343750000E+00, h=0.0313, aderiv=8060.4, x=1.0625, pwr= 8 computing order= 5, npoints= 7, at term= 4 f^( 5)(x( 4)= -5.00000000000000E-01 * f(x( 1)) + 2.00000000000000E+00 * f(x( 2) + -2.50000000000000E+00 * f(x( 3) + -0.00000000000000E+00 * f(x( 4) + 2.50000000000000E+00 * f(x( 5) + -2.00000000000000E+00 * f(x( 6) + 5.00000000000000E-01 * f(x( 7) ) + O(h^ 2)f^( 7)(z) err= 0.00000000000000E+00, h=0.0313, aderiv= 120.0, x=1.0938, pwr= 5 err= 0.00000000000000E+00, h=0.0313, aderiv= 787.5, x=1.0938, pwr= 6 err=-1.64062500000000E+00, h=0.0313, aderiv=3014.6, x=1.0938, pwr= 7 err=-1.43554687500000E+01, h=0.0313, aderiv=8792.7, x=1.0938, pwr= 8 computing order= 5, npoints= 7, at term= 5 f^( 5)(x( 5)= 5.00000000000000E-01 * f(x( 1)) + -4.00000000000000E+00 * f(x( 2) + 1.25000000000000E+01 * f(x( 3) + -2.00000000000000E+01 * f(x( 4) + 1.75000000000000E+01 * f(x( 5) + -8.00000000000000E+00 * f(x( 6) + 1.50000000000000E+00 * f(x( 7) ) + O(h^ 2)f^( 7)(z) err= 0.00000000000000E+00, h=0.0313, aderiv= 120.0, x=1.1250, pwr= 5 err= 0.00000000000000E+00, h=0.0313, aderiv= 810.0, x=1.1250, pwr= 6 err= 8.20312500000000E-01, h=0.0313, aderiv=3189.4, x=1.1250, pwr= 7 err= 7.07519531250000E+00, h=0.0313, aderiv=9568.1, x=1.1250, pwr= 8 computing order= 5, npoints= 7, at term= 6 f^( 5)(x( 6)= 1.50000000000000E+00 * f(x( 1)) + -1.00000000000000E+01 * f(x( 2) + 2.75000000000000E+01 * f(x( 3) + -4.00000000000000E+01 * f(x( 4) + 3.25000000000000E+01 * f(x( 5) + -1.40000000000000E+01 * f(x( 6) + 2.50000000000000E+00 * f(x( 7) ) + O(h^ 2)f^( 7)(z) err= 0.00000000000000E+00, h=0.0313, aderiv= 120.0, x=1.1563, pwr= 5 err= 0.00000000000000E+00, h=0.0313, aderiv= 832.5, x=1.1563, pwr= 6 err= 8.20312500000000E+00, h=0.0313, aderiv=3369.0, x=1.1563, pwr= 7 err= 7.28027343750000E+01, h=0.0313, aderiv=10387.8, x=1.1563, pwr= 8 computing order= 5, npoints= 7, at term= 7 f^( 5)(x( 7)= 2.50000000000000E+00 * f(x( 1)) + -1.60000000000000E+01 * f(x( 2) + 4.25000000000000E+01 * f(x( 3) + -6.00000000000000E+01 * f(x( 4) + 4.75000000000000E+01 * f(x( 5) + -2.00000000000000E+01 * f(x( 6) + 3.50000000000000E+00 * f(x( 7) ) + O(h^ 2)f^( 7)(z) err= 0.00000000000000E+00, h=0.0313, aderiv= 120.0, x=1.1875, pwr= 5 err= 0.00000000000000E+00, h=0.0313, aderiv= 855.0, x=1.1875, pwr= 6 err= 2.05078125000000E+01, h=0.0313, aderiv=3553.6, x=1.1875, pwr= 7 err= 1.84057617187500E+02, h=0.0313, aderiv=11253.0, x=1.1875, pwr= 8 computing order= 6, npoints= 7, at term= 1 f^( 6)(x( 1)= 1.00000000000000E+00 * f(x( 1)) + -6.00000000000000E+00 * f(x( 2) + 1.50000000000000E+01 * f(x( 3) + -2.00000000000000E+01 * f(x( 4) + 1.50000000000000E+01 * f(x( 5) + -6.00000000000000E+00 * f(x( 6) + 1.00000000000000E+00 * f(x( 7) ) + O(h^ 1)f^( 7)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 0.0, x=1.0000, pwr= 5 err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0000, pwr= 6 err=-2.36250000000000E+02, h=0.0156, aderiv=5040.0, x=1.0000, pwr= 7 err=-1.93675781250000E+03, h=0.0156, aderiv=20160.0, x=1.0000, pwr= 8 computing order= 6, npoints= 7, at term= 2 f^( 6)(x( 2)= 1.00000000000000E+00 * f(x( 1)) + -6.00000000000000E+00 * f(x( 2) + 1.50000000000000E+01 * f(x( 3) + -2.00000000000000E+01 * f(x( 4) + 1.50000000000000E+01 * f(x( 5) + -6.00000000000000E+00 * f(x( 6) + 1.00000000000000E+00 * f(x( 7) ) + O(h^ 1)f^( 7)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 0.0, x=1.0156, pwr= 5 err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0156, pwr= 6 err=-1.57500000000000E+02, h=0.0156, aderiv=5118.8, x=1.0156, pwr= 7 err=-1.30183593750000E+03, h=0.0156, aderiv=20794.9, x=1.0156, pwr= 8 computing order= 6, npoints= 7, at term= 3 f^( 6)(x( 3)= 1.00000000000000E+00 * f(x( 1)) + -6.00000000000000E+00 * f(x( 2) + 1.50000000000000E+01 * f(x( 3) + -2.00000000000000E+01 * f(x( 4) + 1.50000000000000E+01 * f(x( 5) + -6.00000000000000E+00 * f(x( 6) + 1.00000000000000E+00 * f(x( 7) ) + O(h^ 1)f^( 7)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 0.0, x=1.0313, pwr= 5 err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0313, pwr= 6 err=-7.87500000000000E+01, h=0.0156, aderiv=5197.5, x=1.0313, pwr= 7 err=-6.57070312500000E+02, h=0.0156, aderiv=21439.7, x=1.0313, pwr= 8 computing order= 6, npoints= 7, at term= 4 f^( 6)(x( 4)= 1.00000000000000E+00 * f(x( 1)) + -6.00000000000000E+00 * f(x( 2) + 1.50000000000000E+01 * f(x( 3) + -2.00000000000000E+01 * f(x( 4) + 1.50000000000000E+01 * f(x( 5) + -6.00000000000000E+00 * f(x( 6) + 1.00000000000000E+00 * f(x( 7) ) + O(h^ 1)f^( 7)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 0.0, x=1.0469, pwr= 5 err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0469, pwr= 6 err= 0.00000000000000E+00, h=0.0156, aderiv=5276.3, x=1.0469, pwr= 7 err=-2.46093750000000E+00, h=0.0156, aderiv=22094.3, x=1.0469, pwr= 8 computing order= 6, npoints= 7, at term= 5 f^( 6)(x( 5)= 1.00000000000000E+00 * f(x( 1)) + -6.00000000000000E+00 * f(x( 2) + 1.50000000000000E+01 * f(x( 3) + -2.00000000000000E+01 * f(x( 4) + 1.50000000000000E+01 * f(x( 5) + -6.00000000000000E+00 * f(x( 6) + 1.00000000000000E+00 * f(x( 7) ) + O(h^ 1)f^( 7)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 0.0, x=1.0625, pwr= 5 err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0625, pwr= 6 err= 7.87500000000000E+01, h=0.0156, aderiv=5355.0, x=1.0625, pwr= 7 err= 6.61992187500000E+02, h=0.0156, aderiv=22758.8, x=1.0625, pwr= 8 computing order= 6, npoints= 7, at term= 6 f^( 6)(x( 6)= 1.00000000000000E+00 * f(x( 1)) + -6.00000000000000E+00 * f(x( 2) + 1.50000000000000E+01 * f(x( 3) + -2.00000000000000E+01 * f(x( 4) + 1.50000000000000E+01 * f(x( 5) + -6.00000000000000E+00 * f(x( 6) + 1.00000000000000E+00 * f(x( 7) ) + O(h^ 1)f^( 7)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 0.0, x=1.0781, pwr= 5 err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0781, pwr= 6 err= 1.57500000000000E+02, h=0.0156, aderiv=5433.8, x=1.0781, pwr= 7 err= 1.33628906250000E+03, h=0.0156, aderiv=23433.0, x=1.0781, pwr= 8 computing order= 6, npoints= 7, at term= 7 f^( 6)(x( 7)= 1.00000000000000E+00 * f(x( 1)) + -6.00000000000000E+00 * f(x( 2) + 1.50000000000000E+01 * f(x( 3) + -2.00000000000000E+01 * f(x( 4) + 1.50000000000000E+01 * f(x( 5) + -6.00000000000000E+00 * f(x( 6) + 1.00000000000000E+00 * f(x( 7) ) + O(h^ 1)f^( 7)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 0.0, x=1.0938, pwr= 5 err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0938, pwr= 6 err= 2.36250000000000E+02, h=0.0156, aderiv=5512.5, x=1.0938, pwr= 7 err= 2.02042968750000E+03, h=0.0156, aderiv=24117.2, x=1.0938, pwr= 8 computing order= 6, npoints= 8, at term= 1 f^( 6)(x( 1)= 4.00000000000000E+00 * f(x( 1)) + -2.70000000000000E+01 * f(x( 2) + 7.80000000000000E+01 * f(x( 3) + -1.25000000000000E+02 * f(x( 4) + 1.20000000000000E+02 * f(x( 5) + -6.90000000000000E+01 * f(x( 6) + 2.20000000000000E+01 * f(x( 7) + -3.00000000000000E+00 * f(x( 8) ) + O(h^ 2)f^( 8)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0000, pwr= 6 err= 0.00000000000000E+00, h=0.0156, aderiv=5040.0, x=1.0000, pwr= 7 err= 5.66022949218750E+01, h=0.0156, aderiv=20160.0, x=1.0000, pwr= 8 err= 5.28794555664063E+02, h=0.0156, aderiv=60480.0, x=1.0000, pwr= 9 computing order= 6, npoints= 8, at term= 2 f^( 6)(x( 2)= 3.00000000000000E+00 * f(x( 1)) + -2.00000000000000E+01 * f(x( 2) + 5.70000000000000E+01 * f(x( 3) + -9.00000000000000E+01 * f(x( 4) + 8.50000000000000E+01 * f(x( 5) + -4.80000000000000E+01 * f(x( 6) + 1.50000000000000E+01 * f(x( 7) + -2.00000000000000E+00 * f(x( 8) ) + O(h^ 2)f^( 8)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0156, pwr= 6 err= 0.00000000000000E+00, h=0.0156, aderiv=5118.8, x=1.0156, pwr= 7 err= 2.70698242187500E+01, h=0.0156, aderiv=20794.9, x=1.0156, pwr= 8 err= 2.54360656738281E+02, h=0.0156, aderiv=63359.5, x=1.0156, pwr= 9 computing order= 6, npoints= 8, at term= 3 f^( 6)(x( 3)= 2.00000000000000E+00 * f(x( 1)) + -1.30000000000000E+01 * f(x( 2) + 3.60000000000000E+01 * f(x( 3) + -5.50000000000000E+01 * f(x( 4) + 5.00000000000000E+01 * f(x( 5) + -2.70000000000000E+01 * f(x( 6) + 8.00000000000000E+00 * f(x( 7) + -1.00000000000000E+00 * f(x( 8) ) + O(h^ 2)f^( 8)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0313, pwr= 6 err= 0.00000000000000E+00, h=0.0156, aderiv=5197.5, x=1.0313, pwr= 7 err= 7.38305664062500E+00, h=0.0156, aderiv=21439.7, x=1.0313, pwr= 8 err= 6.99062805175781E+01, h=0.0156, aderiv=66329.0, x=1.0313, pwr= 9 computing order= 6, npoints= 8, at term= 4 f^( 6)(x( 4)= 1.00000000000000E+00 * f(x( 1)) + -6.00000000000000E+00 * f(x( 2) + 1.50000000000000E+01 * f(x( 3) + -2.00000000000000E+01 * f(x( 4) + 1.50000000000000E+01 * f(x( 5) + -6.00000000000000E+00 * f(x( 6) + 1.00000000000000E+00 * f(x( 7) + 0.00000000000000E+00 * f(x( 8) ) + O(h^ 2)f^( 8)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0469, pwr= 6 err= 0.00000000000000E+00, h=0.0156, aderiv=5276.3, x=1.0469, pwr= 7 err=-2.46093750000000E+00, h=0.0156, aderiv=22094.3, x=1.0469, pwr= 8 err=-2.31866455078125E+01, h=0.0156, aderiv=69389.9, x=1.0469, pwr= 9 computing order= 6, npoints= 8, at term= 5 f^( 6)(x( 5)= -0.00000000000000E+00 * f(x( 1)) + 1.00000000000000E+00 * f(x( 2) + -6.00000000000000E+00 * f(x( 3) + 1.50000000000000E+01 * f(x( 4) + -2.00000000000000E+01 * f(x( 5) + 1.50000000000000E+01 * f(x( 6) + -6.00000000000000E+00 * f(x( 7) + 1.00000000000000E+00 * f(x( 8) ) + O(h^ 2)f^( 8)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0625, pwr= 6 err= 0.00000000000000E+00, h=0.0156, aderiv=5355.0, x=1.0625, pwr= 7 err=-2.46093750000000E+00, h=0.0156, aderiv=22758.8, x=1.0625, pwr= 8 err=-2.35325012207031E+01, h=0.0156, aderiv=72543.5, x=1.0625, pwr= 9 computing order= 6, npoints= 8, at term= 6 f^( 6)(x( 6)= -1.00000000000000E+00 * f(x( 1)) + 8.00000000000000E+00 * f(x( 2) + -2.70000000000000E+01 * f(x( 3) + 5.00000000000000E+01 * f(x( 4) + -5.50000000000000E+01 * f(x( 5) + 3.60000000000000E+01 * f(x( 6) + -1.30000000000000E+01 * f(x( 7) + 2.00000000000000E+00 * f(x( 8) ) + O(h^ 2)f^( 8)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0781, pwr= 6 err= 0.00000000000000E+00, h=0.0156, aderiv=5433.8, x=1.0781, pwr= 7 err= 7.38330078125000E+00, h=0.0156, aderiv=23433.0, x=1.0781, pwr= 8 err= 7.02521362304688E+01, h=0.0156, aderiv=75791.3, x=1.0781, pwr= 9 computing order= 6, npoints= 8, at term= 7 f^( 6)(x( 7)= -2.00000000000000E+00 * f(x( 1)) + 1.50000000000000E+01 * f(x( 2) + -4.80000000000000E+01 * f(x( 3) + 8.50000000000000E+01 * f(x( 4) + -9.00000000000000E+01 * f(x( 5) + 5.70000000000000E+01 * f(x( 6) + -2.00000000000000E+01 * f(x( 7) + 3.00000000000000E+00 * f(x( 8) ) + O(h^ 2)f^( 8)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.0938, pwr= 6 err= 0.00000000000000E+00, h=0.0156, aderiv=5512.5, x=1.0938, pwr= 7 err= 2.70705566406250E+01, h=0.0156, aderiv=24117.2, x=1.0938, pwr= 8 err= 2.59552124023438E+02, h=0.0156, aderiv=79134.5, x=1.0938, pwr= 9 computing order= 6, npoints= 8, at term= 8 f^( 6)(x( 8)= -3.00000000000000E+00 * f(x( 1)) + 2.20000000000000E+01 * f(x( 2) + -6.90000000000000E+01 * f(x( 3) + 1.20000000000000E+02 * f(x( 4) + -1.25000000000000E+02 * f(x( 5) + 7.80000000000000E+01 * f(x( 6) + -2.70000000000000E+01 * f(x( 7) + 4.00000000000000E+00 * f(x( 8) ) + O(h^ 2)f^( 8)(z) err= 0.00000000000000E+00, h=0.0156, aderiv= 720.0, x=1.1094, pwr= 6 err= 0.00000000000000E+00, h=0.0156, aderiv=5591.3, x=1.1094, pwr= 7 err= 5.66005859375000E+01, h=0.0156, aderiv=24811.2, x=1.1094, pwr= 8 err= 5.45751586914063E+02, h=0.0156, aderiv=82574.7, x=1.1094, pwr= 9 numerical accuracy fails at 4,? multiple precision required order=4, points=9, term=1 err=-2.70083546638489E-08, h=0.0625, aderiv= 840.0, x=1.0000, pwr= 7 err=-4.74974513053894E-08, h=0.0625, aderiv=1680.0, x=1.0000, pwr= 8 err=-1.54000859521329E+00, h=0.0625, aderiv=3024.0, x=1.0000, pwr= 9 err=-1.86961269564927E+01, h=0.0625, aderiv=5040.0, x=1.0000, pwr= 10 order=4, points=9, term=2 err= 6.98491930961609E-10, h=0.0625, aderiv=1007.5, x=1.0625, pwr= 7 err= 4.42378222942352E-09, h=0.0625, aderiv=2141.0, x=1.0625, pwr= 8 err=-1.45637511042878E-01, h=0.0625, aderiv=4094.7, x=1.0625, pwr= 9 err=-1.79325342108496E+00, h=0.0625, aderiv=7251.1, x=1.0625, pwr= 10 order=4, points=10, term=1 err=-5.96046447753906E-08, h=0.0625, aderiv=1680.0, x=1.0000, pwr= 8 err=-1.49011611938477E-08, h=0.0625, aderiv=3024.0, x=1.0000, pwr= 9 err= 1.03523258864880E+00, h=0.0625, aderiv=5040.0, x=1.0000, pwr= 10 err= 1.41943017393351E+01, h=0.0625, aderiv=7920.0, x=1.0000, pwr= 11 order=4, points=10, term=2 err= 3.49245965480804E-10, h=0.0625, aderiv=2141.0, x=1.0625, pwr= 8 err= 6.75208866596222E-09, h=0.0625, aderiv=4094.7, x=1.0625, pwr= 9 err= 7.27272122167051E-02, h=0.0625, aderiv=7251.1, x=1.0625, pwr= 10 err= 1.00985349155962E+00, h=0.0625, aderiv=12106.8, x=1.0625, pwr= 11 order=4, points=11, term=1 err=-1.34110450744629E-07, h=0.0625, aderiv=3024.0, x=1.0000, pwr= 9 err= 9.68575477600098E-08, h=0.0625, aderiv=5040.0, x=1.0000, pwr= 10 err=-7.51868031919003E-01, h=0.0625, aderiv=7920.0, x=1.0000, pwr= 11 err=-1.15356941223145E+01, h=0.0625, aderiv=11880.0, x=1.0000, pwr= 12 order=4, points=11, term=2 err= 2.86381691694260E-08, h=0.0625, aderiv=4094.7, x=1.0625, pwr= 9 err= 3.91155481338501E-08, h=0.0625, aderiv=7251.1, x=1.0625, pwr= 10 err=-4.01455238461494E-02, h=0.0625, aderiv=12106.8, x=1.0625, pwr= 11 err=-6.23072045855224E-01, h=0.0625, aderiv=19295.1, x=1.0625, pwr= 12 order=4, points=12, term=1 err= 8.34465026855469E-07, h=0.0625, aderiv=5040.0, x=1.0000, pwr= 10 err=-6.40749931335449E-07, h=0.0625, aderiv=7920.0, x=1.0000, pwr= 11 err= 5.88177889585495E-01, h=0.0625, aderiv=11880.0, x=1.0000, pwr= 12 err= 1.00197954177856E+01, h=0.0625, aderiv=17160.0, x=1.0000, pwr= 13 order=4, points=12, term=2 err= 4.09781932830811E-08, h=0.0625, aderiv=7251.1, x=1.0625, pwr= 10 err=-7.45058059692383E-08, h=0.0625, aderiv=12106.8, x=1.0625, pwr= 11 err= 2.42743995040655E-02, h=0.0625, aderiv=19295.1, x=1.0625, pwr= 12 err= 4.17958608013578E-01, h=0.0625, aderiv=29612.7, x=1.0625, pwr= 13 order=4, points=13, term=1 err=-7.40069038454278E+08, h=0.0625, aderiv=7920.0, x=1.0000, pwr= 11 err=-8.84583515597460E+08, h=0.0625, aderiv=11880.0, x=1.0000, pwr= 12 err=-1.05708671187897E+09, h=0.0625, aderiv=17160.0, x=1.0000, pwr= 13 err=-1.26295135013829E+09, h=0.0625, aderiv=24024.0, x=1.0000, pwr= 14 ERROR sum a(i) /= 0.0, is 1.56913729891508E+03 order=4, points=13, term=2 err= 5.87628468695397E+05, h=0.0625, aderiv=12106.8, x=1.0625, pwr= 11 err= 5.87628468695415E+05, h=0.0625, aderiv=19295.1, x=1.0625, pwr= 12 err= 5.87628452748975E+05, h=0.0625, aderiv=29612.7, x=1.0625, pwr= 13 err= 5.87628166044905E+05, h=0.0625, aderiv=44048.9, x=1.0625, pwr= 14 ERROR sum a(i) /= 0.0, is -8.96649885094335E+00 order=4, points=14, term=1 err= 6.13917240195935E+08, h=0.0625, aderiv=11880.0, x=1.0000, pwr= 12 err= 9.18069147647843E+08, h=0.0625, aderiv=17160.0, x=1.0000, pwr= 13 err= 1.35047502268347E+09, h=0.0625, aderiv=24024.0, x=1.0000, pwr= 14 err= 1.95334770576523E+09, h=0.0625, aderiv=32760.0, x=1.0000, pwr= 15 ERROR sum a(i) /= 0.0, is 6.51281884516447E+01 order=4, points=14, term=2 err=-3.69317456353910E+09, h=0.0625, aderiv=19295.1, x=1.0625, pwr= 12 err=-5.37297194164015E+09, h=0.0625, aderiv=29612.7, x=1.0625, pwr= 13 err=-7.81503560448436E+09, h=0.0625, aderiv=44048.9, x=1.0625, pwr= 14 err=-1.13613147020667E+10, h=0.0625, aderiv=63820.8, x=1.0625, pwr= 15 ERROR sum a(i) /= 0.0, is 7.11259711334740E+02 many formulas checked against Abramowitz and Stegun, Handbook of Mathematical Functions Table 25.2