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=0 f^(1)(x[0])=(1/1h^1)( -1 * f(x[0]) + 1 * f(x[1] ) + O(h^1)f^(2)(z) err=0, h=0.125, deriv=0, x=1, pwr=0 err=0, h=0.125, deriv=1, x=1, pwr=1 err=-0.125, h=0.125, deriv=2, x=1, pwr=2 err=-0.0625, h=0.0625, deriv=2, x=1, pwr=2 err=-0.390625, h=0.125, deriv=3, x=1, pwr=3 err=-0.191406, h=0.0625, deriv=3, x=1, pwr=3 computing order=1, npoints=2, at term=1 f^(1)(x[1])=(1/1h^1)( -1 * f(x[0]) + 1 * f(x[1] ) + O(h^1)f^(2)(z) err=0, h=0.125, deriv=0, x=1.125, pwr=0 err=0, h=0.125, deriv=1, x=1.125, pwr=1 err=0.125, h=0.125, deriv=2.25, x=1.125, pwr=2 err=0.0625, h=0.0625, deriv=2.125, x=1.0625, pwr=2 err=0.40625, h=0.125, deriv=3.79688, x=1.125, pwr=3 err=0.195312, h=0.0625, deriv=3.38672, x=1.0625, pwr=3 computing order=1, npoints=3, at term=0 f^(1)(x[0])=(1/2h^1)( -3 * f(x[0]) + 4 * f(x[1]) + -1 * f(x[2] ) + O(h^2)f^(3)(z) err=0, h=0.125, deriv=1, x=1, pwr=1 err=0, h=0.125, deriv=2, x=1, pwr=2 err=0.03125, h=0.125, deriv=3, x=1, pwr=3 err=0.0078125, h=0.0625, deriv=3, x=1, pwr=3 err=0.136719, h=0.125, deriv=4, x=1, pwr=4 err=0.0327148, h=0.0625, deriv=4, x=1, pwr=4 computing order=1, npoints=3, at term=1 f^(1)(x[1])=(1/2h^1)( -1 * f(x[0]) + 0 * f(x[1]) + 1 * f(x[2] ) + O(h^2)f^(3)(z) err=0, h=0.125, deriv=1, x=1.125, pwr=1 err=0, h=0.125, deriv=2.25, x=1.125, pwr=2 err=-0.015625, h=0.125, deriv=3.79688, x=1.125, pwr=3 err=-0.00390625, h=0.0625, deriv=3.38672, x=1.0625, pwr=3 err=-0.0703125, h=0.125, deriv=5.69531, x=1.125, pwr=4 err=-0.0166016, h=0.0625, deriv=4.79785, x=1.0625, pwr=4 computing order=1, npoints=3, at term=2 f^(1)(x[2])=(1/2h^1)( 1 * f(x[0]) + -4 * f(x[1]) + 3 * f(x[2] ) + O(h^2)f^(3)(z) err=0, h=0.125, deriv=1, x=1.25, pwr=1 err=0, h=0.125, deriv=2.5, x=1.25, pwr=2 err=0.03125, h=0.125, deriv=4.6875, x=1.25, pwr=3 err=0.0078125, h=0.0625, deriv=3.79688, x=1.125, pwr=3 err=0.144531, h=0.125, deriv=7.8125, x=1.25, pwr=4 err=0.0336914, h=0.0625, deriv=5.69531, x=1.125, pwr=4 computing order=1, npoints=4, at term=0 f^(1)(x[0])=(1/6h^1)( -11 * f(x[0]) + 18 * f(x[1]) + -9 * f(x[2]) + 2 * f(x[3] ) + O(h^3)f^(4)(z) err=0, h=0.125, deriv=2, x=1, pwr=2 err=0, h=0.125, deriv=3, x=1, pwr=3 err=-0.0117188, h=0.125, deriv=4, x=1, pwr=4 err=-0.00146484, h=0.0625, deriv=4, x=1, pwr=4 err=-0.0673828, h=0.125, deriv=5, x=1, pwr=5 err=-0.00787354, h=0.0625, deriv=5, x=1, pwr=5 computing order=1, npoints=4, at term=1 f^(1)(x[1])=(1/6h^1)( -2 * f(x[0]) + -3 * f(x[1]) + 6 * f(x[2]) + -1 * f(x[3] ) + O(h^3)f^(4)(z) err=0, h=0.125, deriv=2.25, x=1.125, pwr=2 err=0, h=0.125, deriv=3.79688, x=1.125, pwr=3 err=0.00390625, h=0.125, deriv=5.69531, x=1.125, pwr=4 err=0.000488281, h=0.0625, deriv=4.79785, x=1.0625, pwr=4 err=0.0229492, h=0.125, deriv=8.00903, x=1.125, pwr=5 err=0.00265503, h=0.0625, deriv=6.37215, x=1.0625, pwr=5 computing order=1, npoints=4, at term=2 f^(1)(x[2])=(1/6h^1)( 1 * f(x[0]) + -6 * f(x[1]) + 3 * f(x[2]) + 2 * f(x[3] ) + O(h^3)f^(4)(z) err=0, h=0.125, deriv=2.5, x=1.25, pwr=2 err=0, h=0.125, deriv=4.6875, x=1.25, pwr=3 err=-0.00390625, h=0.125, deriv=7.8125, x=1.25, pwr=4 err=-0.000488281, h=0.0625, deriv=5.69531, x=1.125, pwr=4 err=-0.0234375, h=0.125, deriv=12.207, x=1.25, pwr=5 err=-0.00268555, h=0.0625, deriv=8.00903, x=1.125, pwr=5 computing order=1, npoints=4, at term=3 f^(1)(x[3])=(1/6h^1)( -2 * f(x[0]) + 9 * f(x[1]) + -18 * f(x[2]) + 11 * f(x[3] ) + O(h^3)f^(4)(z) err=0, h=0.125, deriv=2.75, x=1.375, pwr=2 err=0, h=0.125, deriv=5.67188, x=1.375, pwr=3 err=0.0117188, h=0.125, deriv=10.3984, x=1.375, pwr=4 err=0.00146484, h=0.0625, deriv=6.69824, x=1.1875, pwr=4 err=0.0717773, h=0.125, deriv=17.8723, x=1.375, pwr=5 err=0.00814819, h=0.0625, deriv=9.9427, x=1.1875, pwr=5 computing order=1, npoints=5, at term=0 f^(1)(x[0])=(1/12h^1)( -25 * f(x[0]) + 48 * f(x[1]) + -36 * f(x[2]) + 16 * f(x[3]) + -3 * f(x[4] ) + O(h^4)f^(5)(z) err=0, h=0.125, deriv=3, x=1, pwr=3 err=0, h=0.125, deriv=4, x=1, pwr=4 err=0.00585938, h=0.125, deriv=5, x=1, pwr=5 err=0.000366211, h=0.0625, deriv=5, x=1, pwr=5 err=0.0424805, h=0.125, deriv=6, x=1, pwr=6 err=0.00242615, h=0.0625, deriv=6, x=1, pwr=6 computing order=1, npoints=5, at term=1 f^(1)(x[1])=(1/12h^1)( -3 * f(x[0]) + -10 * f(x[1]) + 18 * f(x[2]) + -6 * f(x[3]) + 1 * f(x[4] ) + O(h^4)f^(5)(z) err=0, h=0.125, deriv=3.79688, x=1.125, pwr=3 err=0, h=0.125, deriv=5.69531, x=1.125, pwr=4 err=-0.00146484, h=0.125, deriv=8.00903, x=1.125, pwr=5 err=-9.15527e-05, h=0.0625, deriv=6.37215, x=1.0625, pwr=5 err=-0.0108032, h=0.125, deriv=10.8122, x=1.125, pwr=6 err=-0.000612259, h=0.0625, deriv=8.12449, x=1.0625, pwr=6 computing order=1, npoints=5, at term=2 f^(1)(x[2])=(1/12h^1)( 1 * f(x[0]) + -8 * f(x[1]) + 0 * f(x[2]) + 8 * f(x[3]) + -1 * f(x[4] ) + O(h^4)f^(5)(z) err=0, h=0.125, deriv=4.6875, x=1.25, pwr=3 err=0, h=0.125, deriv=7.8125, x=1.25, pwr=4 err=0.000976562, h=0.125, deriv=12.207, x=1.25, pwr=5 err=6.10352e-05, h=0.0625, deriv=8.00903, x=1.125, pwr=5 err=0.00732422, h=0.125, deriv=18.3105, x=1.25, pwr=6 err=0.000411987, h=0.0625, deriv=10.8122, x=1.125, pwr=6 computing order=1, npoints=5, at term=3 f^(1)(x[3])=(1/12h^1)( -1 * f(x[0]) + 6 * f(x[1]) + -18 * f(x[2]) + 10 * f(x[3]) + 3 * f(x[4] ) + O(h^4)f^(5)(z) err=0, h=0.125, deriv=5.67188, x=1.375, pwr=3 err=0, h=0.125, deriv=10.3984, x=1.375, pwr=4 err=-0.00146484, h=0.125, deriv=17.8723, x=1.375, pwr=5 err=-9.15527e-05, h=0.0625, deriv=9.9427, x=1.1875, pwr=5 err=-0.0111694, h=0.125, deriv=29.4893, x=1.375, pwr=6 err=-0.000623703, h=0.0625, deriv=14.1684, x=1.1875, pwr=6 computing order=1, npoints=5, at term=4 f^(1)(x[4])=(1/12h^1)( 3 * f(x[0]) + -16 * f(x[1]) + 36 * f(x[2]) + -48 * f(x[3]) + 25 * f(x[4] ) + O(h^4)f^(5)(z) err=0, h=0.125, deriv=6.75, x=1.5, pwr=3 err=0, h=0.125, deriv=13.5, x=1.5, pwr=4 err=0.00585938, h=0.125, deriv=25.3125, x=1.5, pwr=5 err=0.000366211, h=0.0625, deriv=12.207, x=1.25, pwr=5 err=0.0454102, h=0.125, deriv=45.5625, x=1.5, pwr=6 err=0.0025177, h=0.0625, deriv=18.3105, x=1.25, pwr=6 computing order=1, npoints=6, at term=0 f^(1)(x[0])=(1/60h^1)( -137 * f(x[0]) + 300 * f(x[1]) + -300 * f(x[2]) + 200 * f(x[3]) + -75 * f(x[4]) + 12 * f(x[5] ) + O(h^5)f^(6)(z) err=0, h=0.125, deriv=4, x=1, pwr=4 err=0, h=0.125, deriv=5, x=1, pwr=5 err=-0.00366211, h=0.125, deriv=6, x=1, pwr=6 err=-0.000114441, h=0.0625, deriv=6, x=1, pwr=6 err=-0.0325012, h=0.125, deriv=7, x=1, pwr=7 err=-0.000908375, h=0.0625, deriv=7, x=1, pwr=7 computing order=1, npoints=6, at term=1 f^(1)(x[1])=(1/60h^1)( -12 * f(x[0]) + -65 * f(x[1]) + 120 * f(x[2]) + -60 * f(x[3]) + 20 * f(x[4]) + -3 * f(x[5] ) + O(h^5)f^(6)(z) err=0, h=0.125, deriv=5.69531, x=1.125, pwr=4 err=0, h=0.125, deriv=8.00903, x=1.125, pwr=5 err=0.000732422, h=0.125, deriv=10.8122, x=1.125, pwr=6 err=2.28882e-05, h=0.0625, deriv=8.12449, x=1.0625, pwr=6 err=0.0065918, h=0.125, deriv=14.191, x=1.125, pwr=7 err=0.000183105, h=0.0625, deriv=10.071, x=1.0625, pwr=7 computing order=1, npoints=6, at term=2 f^(1)(x[2])=(1/60h^1)( 3 * f(x[0]) + -30 * f(x[1]) + -20 * f(x[2]) + 60 * f(x[3]) + -15 * f(x[4]) + 2 * f(x[5] ) + O(h^5)f^(6)(z) err=0, h=0.125, deriv=7.8125, x=1.25, pwr=4 err=0, h=0.125, deriv=12.207, x=1.25, pwr=5 err=-0.000366211, h=0.125, deriv=18.3105, x=1.25, pwr=6 err=-1.14441e-05, h=0.0625, deriv=10.8122, x=1.125, pwr=6 err=-0.00334167, h=0.125, deriv=26.7029, x=1.25, pwr=7 err=-9.2268e-05, h=0.0625, deriv=14.191, x=1.125, pwr=7 computing order=1, npoints=6, at term=3 f^(1)(x[3])=(1/60h^1)( -2 * f(x[0]) + 15 * f(x[1]) + -60 * f(x[2]) + 20 * f(x[3]) + 30 * f(x[4]) + -3 * f(x[5] ) + O(h^5)f^(6)(z) err=0, h=0.125, deriv=10.3984, x=1.375, pwr=4 err=0, h=0.125, deriv=17.8723, x=1.375, pwr=5 err=0.000366211, h=0.125, deriv=29.4893, x=1.375, pwr=6 err=1.14441e-05, h=0.0625, deriv=14.1684, x=1.1875, pwr=6 err=0.00338745, h=0.125, deriv=47.3058, x=1.375, pwr=7 err=9.29832e-05, h=0.0625, deriv=19.6291, x=1.1875, pwr=7 computing order=1, npoints=6, at term=4 f^(1)(x[4])=(1/60h^1)( 3 * f(x[0]) + -20 * f(x[1]) + 60 * f(x[2]) + -120 * f(x[3]) + 65 * f(x[4]) + 12 * f(x[5] ) + O(h^5)f^(6)(z) err=0, h=0.125, deriv=13.5, x=1.5, pwr=4 err=0, h=0.125, deriv=25.3125, x=1.5, pwr=5 err=-0.000732422, h=0.125, deriv=45.5625, x=1.5, pwr=6 err=-2.28882e-05, h=0.0625, deriv=18.3105, x=1.25, pwr=6 err=-0.00686646, h=0.125, deriv=79.7344, x=1.5, pwr=7 err=-0.000187397, h=0.0625, deriv=26.7029, x=1.25, pwr=7 computing order=1, npoints=6, at term=5 f^(1)(x[5])=(1/60h^1)( -12 * f(x[0]) + 75 * f(x[1]) + -200 * f(x[2]) + 300 * f(x[3]) + -300 * f(x[4]) + 137 * f(x[5] ) + O(h^5)f^(6)(z) err=0, h=0.125, deriv=17.1641, x=1.625, pwr=4 err=0, h=0.125, deriv=34.8645, x=1.625, pwr=5 err=0.00366211, h=0.125, deriv=67.9858, x=1.625, pwr=6 err=0.000114441, h=0.0625, deriv=23.3694, x=1.3125, pwr=6 err=0.03479, h=0.125, deriv=128.89, x=1.625, pwr=7 err=0.000944138, h=0.0625, deriv=35.7844, x=1.3125, pwr=7 computing order=1, npoints=7, at term=0 f^(1)(x[0])=(1/60h^1)( -147 * f(x[0]) + 360 * f(x[1]) + -450 * f(x[2]) + 400 * f(x[3]) + -225 * f(x[4]) + 72 * f(x[5]) + -10 * f(x[6] ) + O(h^6)f^(7)(z) err=0, h=0.125, deriv=5, x=1, pwr=5 err=0, h=0.125, deriv=6, x=1, pwr=6 err=0.00274658, h=0.125, deriv=7, x=1, pwr=7 err=4.29153e-05, h=0.0625, deriv=7, x=1, pwr=7 err=0.0291824, h=0.125, deriv=8, x=1, pwr=8 err=0.000399649, h=0.0625, deriv=8, x=1, pwr=8 computing order=1, npoints=7, at term=1 f^(1)(x[1])=(1/60h^1)( -10 * f(x[0]) + -77 * f(x[1]) + 150 * f(x[2]) + -100 * f(x[3]) + 50 * f(x[4]) + -15 * f(x[5]) + 2 * f(x[6] ) + O(h^6)f^(7)(z) err=0, h=0.125, deriv=8.00903, x=1.125, pwr=5 err=0, h=0.125, deriv=10.8122, x=1.125, pwr=6 err=-0.000457764, h=0.125, deriv=14.191, x=1.125, pwr=7 err=-7.15256e-06, h=0.0625, deriv=10.071, x=1.0625, pwr=7 err=-0.00492096, h=0.125, deriv=18.2456, x=1.125, pwr=8 err=-6.70552e-05, h=0.0625, deriv=12.229, x=1.0625, pwr=8 computing order=1, npoints=7, at term=2 f^(1)(x[2])=(1/60h^1)( 2 * f(x[0]) + -24 * f(x[1]) + -35 * f(x[2]) + 80 * f(x[3]) + -30 * f(x[4]) + 8 * f(x[5]) + -1 * f(x[6] ) + O(h^6)f^(7)(z) err=0, h=0.125, deriv=12.207, x=1.25, pwr=5 err=0, h=0.125, deriv=18.3105, x=1.25, pwr=6 err=0.000183105, h=0.125, deriv=26.7029, x=1.25, pwr=7 err=2.86102e-06, h=0.0625, deriv=14.191, x=1.125, pwr=7 err=0.00199127, h=0.125, deriv=38.147, x=1.25, pwr=8 err=2.70009e-05, h=0.0625, deriv=18.2456, x=1.125, pwr=8 computing order=1, npoints=7, at term=3 f^(1)(x[3])=(1/60h^1)( -1 * f(x[0]) + 9 * f(x[1]) + -45 * f(x[2]) + 0 * f(x[3]) + 45 * f(x[4]) + -9 * f(x[5]) + 1 * f(x[6] ) + O(h^6)f^(7)(z) err=0, h=0.125, deriv=17.8723, x=1.375, pwr=5 err=0, h=0.125, deriv=29.4893, x=1.375, pwr=6 err=-0.000137329, h=0.125, deriv=47.3058, x=1.375, pwr=7 err=-2.14577e-06, h=0.0625, deriv=19.6291, x=1.1875, pwr=7 err=-0.00151062, h=0.125, deriv=74.3377, x=1.375, pwr=8 err=-2.03848e-05, h=0.0625, deriv=26.6395, x=1.1875, pwr=8 computing order=1, npoints=7, at term=4 f^(1)(x[4])=(1/60h^1)( 1 * f(x[0]) + -8 * f(x[1]) + 30 * f(x[2]) + -80 * f(x[3]) + 35 * f(x[4]) + 24 * f(x[5]) + -2 * f(x[6] ) + O(h^6)f^(7)(z) err=0, h=0.125, deriv=25.3125, x=1.5, pwr=5 err=0, h=0.125, deriv=45.5625, x=1.5, pwr=6 err=0.000183105, h=0.125, deriv=79.7344, x=1.5, pwr=7 err=2.86102e-06, h=0.0625, deriv=26.7029, x=1.25, pwr=7 err=0.00203705, h=0.125, deriv=136.688, x=1.5, pwr=8 err=2.73585e-05, h=0.0625, deriv=38.147, x=1.25, pwr=8 computing order=1, npoints=7, at term=5 f^(1)(x[5])=(1/60h^1)( -2 * f(x[0]) + 15 * f(x[1]) + -50 * f(x[2]) + 100 * f(x[3]) + -150 * f(x[4]) + 77 * f(x[5]) + 10 * f(x[6] ) + O(h^6)f^(7)(z) err=0, h=0.125, deriv=34.8645, x=1.625, pwr=5 err=0, h=0.125, deriv=67.9858, x=1.625, pwr=6 err=-0.000457764, h=0.125, deriv=128.89, x=1.625, pwr=7 err=-7.15256e-06, h=0.0625, deriv=35.7844, x=1.3125, pwr=7 err=-0.00514984, h=0.125, deriv=239.367, x=1.625, pwr=8 err=-6.88434e-05, h=0.0625, deriv=53.6766, x=1.3125, pwr=8 computing order=1, npoints=7, at term=6 f^(1)(x[6])=(1/60h^1)( 10 * f(x[0]) + -72 * f(x[1]) + 225 * f(x[2]) + -400 * f(x[3]) + 450 * f(x[4]) + -360 * f(x[5]) + 147 * f(x[6] ) + O(h^6)f^(7)(z) err=0, h=0.125, deriv=46.8945, x=1.75, pwr=5 err=0, h=0.125, deriv=98.4785, x=1.75, pwr=6 err=0.00274658, h=0.125, deriv=201.06, x=1.75, pwr=7 err=4.29153e-05, h=0.0625, deriv=47.3058, x=1.375, pwr=7 err=0.0312424, h=0.125, deriv=402.121, x=1.75, pwr=8 err=0.000415742, h=0.0625, deriv=74.3377, x=1.375, pwr=8 computing order=1, npoints=8, at term=0 f^(1)(x[0])=(1/420h^1)( -1089 * f(x[0]) + 2940 * f(x[1]) + -4410 * f(x[2]) + 4900 * f(x[3]) + -3675 * f(x[4]) + 1764 * f(x[5]) + -490 * f(x[6]) + 60 * f(x[7] ) + O(h^7)f^(8)(z) err=0, h=0.125, deriv=6, x=1, pwr=6 err=0, h=0.125, deriv=7, x=1, pwr=7 err=-0.00240326, h=0.125, deriv=8, x=1, pwr=8 err=-1.87755e-05, h=0.0625, deriv=8, x=1, pwr=8 err=-0.0300407, h=0.125, deriv=9, x=1, pwr=9 err=-0.000201836, h=0.0625, deriv=9, x=1, pwr=9 computing order=1, npoints=8, at term=1 f^(1)(x[1])=(1/420h^1)( -60 * f(x[0]) + -609 * f(x[1]) + 1260 * f(x[2]) + -1050 * f(x[3]) + 700 * f(x[4]) + -315 * f(x[5]) + 84 * f(x[6]) + -10 * f(x[7] ) + O(h^7)f^(8)(z) err=0, h=0.125, deriv=10.8122, x=1.125, pwr=6 err=0, h=0.125, deriv=14.191, x=1.125, pwr=7 err=0.000343323, h=0.125, deriv=18.2456, x=1.125, pwr=8 err=2.68221e-06, h=0.0625, deriv=12.229, x=1.0625, pwr=8 err=0.00433445, h=0.125, deriv=23.0921, x=1.125, pwr=9 err=2.90014e-05, h=0.0625, deriv=14.6175, x=1.0625, pwr=9 computing order=1, npoints=8, at term=2 f^(1)(x[2])=(1/420h^1)( 10 * f(x[0]) + -140 * f(x[1]) + -329 * f(x[2]) + 700 * f(x[3]) + -350 * f(x[4]) + 140 * f(x[5]) + -35 * f(x[6]) + 4 * f(x[7] ) + O(h^7)f^(8)(z) err=0, h=0.125, deriv=18.3105, x=1.25, pwr=6 err=0, h=0.125, deriv=26.7029, x=1.25, pwr=7 err=-0.000114441, h=0.125, deriv=38.147, x=1.25, pwr=8 err=-8.9407e-07, h=0.0625, deriv=18.2456, x=1.125, pwr=8 err=-0.00145912, h=0.125, deriv=53.6442, x=1.25, pwr=9 err=-9.72301e-06, h=0.0625, deriv=23.0921, x=1.125, pwr=9 computing order=1, npoints=8, at term=3 f^(1)(x[3])=(1/420h^1)( -4 * f(x[0]) + 42 * f(x[1]) + -252 * f(x[2]) + -105 * f(x[3]) + 420 * f(x[4]) + -126 * f(x[5]) + 28 * f(x[6]) + -3 * f(x[7] ) + O(h^7)f^(8)(z) err=0, h=0.125, deriv=29.4893, x=1.375, pwr=6 err=0, h=0.125, deriv=47.3058, x=1.375, pwr=7 err=6.86646e-05, h=0.125, deriv=74.3377, x=1.375, pwr=8 err=5.36442e-07, h=0.0625, deriv=26.6395, x=1.1875, pwr=8 err=0.000884056, h=0.125, deriv=114.991, x=1.375, pwr=9 err=5.86733e-06, h=0.0625, deriv=35.5886, x=1.1875, pwr=9 computing order=1, npoints=8, at term=4 f^(1)(x[4])=(1/420h^1)( 3 * f(x[0]) + -28 * f(x[1]) + 126 * f(x[2]) + -420 * f(x[3]) + 105 * f(x[4]) + 252 * f(x[5]) + -42 * f(x[6]) + 4 * f(x[7] ) + O(h^7)f^(8)(z) err=0, h=0.125, deriv=45.5625, x=1.5, pwr=6 err=0, h=0.125, deriv=79.7344, x=1.5, pwr=7 err=-6.86646e-05, h=0.125, deriv=136.688, x=1.5, pwr=8 err=-5.36442e-07, h=0.0625, deriv=38.147, x=1.25, pwr=8 err=-0.000892639, h=0.125, deriv=230.66, x=1.5, pwr=9 err=-5.90086e-06, h=0.0625, deriv=53.6442, x=1.25, pwr=9 computing order=1, npoints=8, at term=5 f^(1)(x[5])=(1/420h^1)( -4 * f(x[0]) + 35 * f(x[1]) + -140 * f(x[2]) + 350 * f(x[3]) + -700 * f(x[4]) + 329 * f(x[5]) + 140 * f(x[6]) + -10 * f(x[7] ) + O(h^7)f^(8)(z) err=0, h=0.125, deriv=67.9858, x=1.625, pwr=6 err=0, h=0.125, deriv=128.89, x=1.625, pwr=7 err=0.000114441, h=0.125, deriv=239.367, x=1.625, pwr=8 err=8.9407e-07, h=0.0625, deriv=53.6766, x=1.3125, pwr=8 err=0.00150204, h=0.125, deriv=437.592, x=1.625, pwr=9 err=9.89065e-06, h=0.0625, deriv=79.2569, x=1.3125, pwr=9 computing order=1, npoints=8, at term=6 f^(1)(x[6])=(1/420h^1)( 10 * f(x[0]) + -84 * f(x[1]) + 315 * f(x[2]) + -700 * f(x[3]) + 1050 * f(x[4]) + -1260 * f(x[5]) + 609 * f(x[6]) + 60 * f(x[7] ) + O(h^7)f^(8)(z) err=0, h=0.125, deriv=98.4785, x=1.75, pwr=6 err=0, h=0.125, deriv=201.06, x=1.75, pwr=7 err=-0.000343323, h=0.125, deriv=402.121, x=1.75, pwr=8 err=-2.68221e-06, h=0.0625, deriv=74.3377, x=1.375, pwr=8 err=-0.00454903, h=0.125, deriv=791.675, x=1.75, pwr=9 err=-2.98396e-05, h=0.0625, deriv=114.991, x=1.375, pwr=9 computing order=1, npoints=8, at term=7 f^(1)(x[7])=(1/420h^1)( -60 * f(x[0]) + 490 * f(x[1]) + -1764 * f(x[2]) + 3675 * f(x[3]) + -4900 * f(x[4]) + 4410 * f(x[5]) + -2940 * f(x[6]) + 1089 * f(x[7] ) + O(h^7)f^(8)(z) err=0, h=0.125, deriv=139.046, x=1.875, pwr=6 err=0, h=0.125, deriv=304.163, x=1.875, pwr=7 err=0.00240326, h=0.125, deriv=651.777, x=1.875, pwr=8 err=1.87755e-05, h=0.0625, deriv=101.472, x=1.4375, pwr=8 err=0.0321436, h=0.125, deriv=1374.84, x=1.875, pwr=9 err=0.00021005, h=0.0625, deriv=164.099, x=1.4375, pwr=9 computing order=1, npoints=9, at term=0 f^(1)(x[0])=(1/840h^1)( -2283 * f(x[0]) + 6720 * f(x[1]) + -11760 * f(x[2]) + 15680 * f(x[3]) + -14700 * f(x[4]) + 9408 * f(x[5]) + -3920 * f(x[6]) + 960 * f(x[7]) + -105 * f(x[8] ) + O(h^8)f^(9)(z) err=0, h=0.125, deriv=7, x=1, pwr=7 err=0, h=0.125, deriv=8, x=1, pwr=8 err=0.00240326, h=0.125, deriv=9, x=1, pwr=9 err=9.38773e-06, h=0.0625, deriv=9, x=1, pwr=9 err=0.0348473, h=0.125, deriv=10, x=1, pwr=10 err=0.000115, h=0.0625, deriv=10, x=1, pwr=10 computing order=1, npoints=9, at term=1 f^(1)(x[1])=(1/840h^1)( -105 * f(x[0]) + -1338 * f(x[1]) + 2940 * f(x[2]) + -2940 * f(x[3]) + 2450 * f(x[4]) + -1470 * f(x[5]) + 588 * f(x[6]) + -140 * f(x[7]) + 15 * f(x[8] ) + O(h^8)f^(9)(z) err=0, h=0.125, deriv=14.191, x=1.125, pwr=7 err=0, h=0.125, deriv=18.2456, x=1.125, pwr=8 err=-0.000300407, h=0.125, deriv=23.0921, x=1.125, pwr=9 err=-1.17347e-06, h=0.0625, deriv=14.6175, x=1.0625, pwr=9 err=-0.00439346, h=0.125, deriv=28.8651, x=1.125, pwr=10 err=-1.44483e-05, h=0.0625, deriv=17.2568, x=1.0625, pwr=10 computing order=1, npoints=9, at term=2 f^(1)(x[2])=(1/840h^1)( 15 * f(x[0]) + -240 * f(x[1]) + -798 * f(x[2]) + 1680 * f(x[3]) + -1050 * f(x[4]) + 560 * f(x[5]) + -210 * f(x[6]) + 48 * f(x[7]) + -5 * f(x[8] ) + O(h^8)f^(9)(z) err=0, h=0.125, deriv=26.7029, x=1.25, pwr=7 err=0, h=0.125, deriv=38.147, x=1.25, pwr=8 err=8.58307e-05, h=0.125, deriv=53.6442, x=1.25, pwr=9 err=3.35276e-07, h=0.0625, deriv=23.0921, x=1.125, pwr=9 err=0.001266, h=0.125, deriv=74.5058, x=1.25, pwr=10 err=4.14904e-06, h=0.0625, deriv=28.8651, x=1.125, pwr=10 computing order=1, npoints=9, at term=3 f^(1)(x[3])=(1/840h^1)( -5 * f(x[0]) + 60 * f(x[1]) + -420 * f(x[2]) + -378 * f(x[3]) + 1050 * f(x[4]) + -420 * f(x[5]) + 140 * f(x[6]) + -30 * f(x[7]) + 3 * f(x[8] ) + O(h^8)f^(9)(z) err=0, h=0.125, deriv=47.3058, x=1.375, pwr=7 err=0, h=0.125, deriv=74.3377, x=1.375, pwr=8 err=-4.29153e-05, h=0.125, deriv=114.991, x=1.375, pwr=9 err=-1.67638e-07, h=0.0625, deriv=35.5886, x=1.1875, pwr=9 err=-0.000638366, h=0.125, deriv=175.681, x=1.375, pwr=10 err=-2.085e-06, h=0.0625, deriv=46.9572, x=1.1875, pwr=10 computing order=1, npoints=9, at term=4 f^(1)(x[4])=(1/840h^1)( 3 * f(x[0]) + -32 * f(x[1]) + 168 * f(x[2]) + -672 * f(x[3]) + 0 * f(x[4]) + 672 * f(x[5]) + -168 * f(x[6]) + 32 * f(x[7]) + -3 * f(x[8] ) + O(h^8)f^(9)(z) err=0, h=0.125, deriv=79.7344, x=1.5, pwr=7 err=0, h=0.125, deriv=136.688, x=1.5, pwr=8 err=3.43323e-05, h=0.125, deriv=230.66, x=1.5, pwr=9 err=1.3411e-07, h=0.0625, deriv=53.6442, x=1.25, pwr=9 err=0.000514984, h=0.125, deriv=384.434, x=1.5, pwr=10 err=1.67638e-06, h=0.0625, deriv=74.5058, x=1.25, pwr=10 computing order=1, npoints=9, at term=5 f^(1)(x[5])=(1/840h^1)( -3 * f(x[0]) + 30 * f(x[1]) + -140 * f(x[2]) + 420 * f(x[3]) + -1050 * f(x[4]) + 378 * f(x[5]) + 420 * f(x[6]) + -60 * f(x[7]) + 5 * f(x[8] ) + O(h^8)f^(9)(z) err=0, h=0.125, deriv=128.89, x=1.625, pwr=7 err=0, h=0.125, deriv=239.367, x=1.625, pwr=8 err=-4.29153e-05, h=0.125, deriv=437.592, x=1.625, pwr=9 err=-1.67638e-07, h=0.0625, deriv=79.2569, x=1.3125, pwr=9 err=-0.000649095, h=0.125, deriv=790.097, x=1.625, pwr=10 err=-2.10595e-06, h=0.0625, deriv=115.583, x=1.3125, pwr=10 computing order=1, npoints=9, at term=6 f^(1)(x[6])=(1/840h^1)( 5 * f(x[0]) + -48 * f(x[1]) + 210 * f(x[2]) + -560 * f(x[3]) + 1050 * f(x[4]) + -1680 * f(x[5]) + 798 * f(x[6]) + 240 * f(x[7]) + -15 * f(x[8] ) + O(h^8)f^(9)(z) err=0, h=0.125, deriv=201.06, x=1.75, pwr=7 err=0, h=0.125, deriv=402.121, x=1.75, pwr=8 err=8.58307e-05, h=0.125, deriv=791.675, x=1.75, pwr=9 err=3.35276e-07, h=0.0625, deriv=114.991, x=1.375, pwr=9 err=0.00130892, h=0.125, deriv=1539.37, x=1.75, pwr=10 err=4.23286e-06, h=0.0625, deriv=175.681, x=1.375, pwr=10 computing order=1, npoints=9, at term=7 f^(1)(x[7])=(1/840h^1)( -15 * f(x[0]) + 140 * f(x[1]) + -588 * f(x[2]) + 1470 * f(x[3]) + -2450 * f(x[4]) + 2940 * f(x[5]) + -2940 * f(x[6]) + 1338 * f(x[7]) + 105 * f(x[8] ) + O(h^8)f^(9)(z) err=0, h=0.125, deriv=304.163, x=1.875, pwr=7 err=0, h=0.125, deriv=651.777, x=1.875, pwr=8 err=-0.000300407, h=0.125, deriv=1374.84, x=1.875, pwr=9 err=-1.17347e-06, h=0.0625, deriv=164.099, x=1.4375, pwr=9 err=-0.00461876, h=0.125, deriv=2864.25, x=1.875, pwr=10 err=-1.48884e-05, h=0.0625, deriv=262.102, x=1.4375, pwr=10 computing order=1, npoints=9, at term=8 f^(1)(x[8])=(1/840h^1)( 105 * f(x[0]) + -960 * f(x[1]) + 3920 * f(x[2]) + -9408 * f(x[3]) + 14700 * f(x[4]) + -15680 * f(x[5]) + 11760 * f(x[6]) + -6720 * f(x[7]) + 2283 * f(x[8] ) + O(h^8)f^(9)(z) err=0, h=0.125, deriv=448, x=2, pwr=7 err=0, h=0.125, deriv=1024, x=2, pwr=8 err=0.00240326, h=0.125, deriv=2304, x=2, pwr=9 err=9.38773e-06, h=0.0625, deriv=230.66, x=1.5, pwr=9 err=0.0372505, h=0.125, deriv=5120, x=2, pwr=10 err=0.000119694, h=0.0625, deriv=384.434, x=1.5, pwr=10 computing order=1, npoints=10, at term=0 f^(1)(x[0])=(1/2520h^1)( -7129 * f(x[0]) + 22680 * f(x[1]) + -45360 * f(x[2]) + 70560 * f(x[3]) + -79380 * f(x[4]) + 63504 * f(x[5]) + -35280 * f(x[6]) + 12960 * f(x[7]) + -2835 * f(x[8]) + 280 * f(x[9] ) + O(h^9)f^(10)(z) err=0, h=0.125, deriv=8, x=1, pwr=8 err=0, h=0.125, deriv=9, x=1, pwr=9 err=-0.00270367, h=0.125, deriv=10, x=1, pwr=10 err=-5.2806e-06, h=0.0625, deriv=10, x=1, pwr=10 err=-0.0449485, h=0.125, deriv=11, x=1, pwr=11 err=-7.29383e-05, h=0.0625, deriv=11, x=1, pwr=11 computing order=1, npoints=10, at term=1 f^(1)(x[1])=(1/2520h^1)( -280 * f(x[0]) + -4329 * f(x[1]) + 10080 * f(x[2]) + -11760 * f(x[3]) + 11760 * f(x[4]) + -8820 * f(x[5]) + 4704 * f(x[6]) + -1680 * f(x[7]) + 360 * f(x[8]) + -35 * f(x[9] ) + O(h^9)f^(10)(z) err=0, h=0.125, deriv=18.2456, x=1.125, pwr=8 err=0, h=0.125, deriv=23.0921, x=1.125, pwr=9 err=0.000300407, h=0.125, deriv=28.8651, x=1.125, pwr=10 err=5.86733e-07, h=0.0625, deriv=17.2568, x=1.0625, pwr=10 err=0.00503182, h=0.125, deriv=35.7205, x=1.125, pwr=11 err=8.14092e-06, h=0.0625, deriv=20.1689, x=1.0625, pwr=11 computing order=1, npoints=10, at term=2 f^(1)(x[2])=(1/2520h^1)( 35 * f(x[0]) + -630 * f(x[1]) + -2754 * f(x[2]) + 5880 * f(x[3]) + -4410 * f(x[4]) + 2940 * f(x[5]) + -1470 * f(x[6]) + 504 * f(x[7]) + -105 * f(x[8]) + 10 * f(x[9] ) + O(h^9)f^(10)(z) err=0, h=0.125, deriv=38.147, x=1.25, pwr=8 err=0, h=0.125, deriv=53.6442, x=1.25, pwr=9 err=-7.51019e-05, h=0.125, deriv=74.5058, x=1.25, pwr=10 err=-1.46683e-07, h=0.0625, deriv=28.8651, x=1.125, pwr=10 err=-0.00126734, h=0.125, deriv=102.445, x=1.25, pwr=11 err=-2.0444e-06, h=0.0625, deriv=35.7205, x=1.125, pwr=11 computing order=1, npoints=10, at term=3 f^(1)(x[3])=(1/2520h^1)( -10 * f(x[0]) + 135 * f(x[1]) + -1080 * f(x[2]) + -1554 * f(x[3]) + 3780 * f(x[4]) + -1890 * f(x[5]) + 840 * f(x[6]) + -270 * f(x[7]) + 54 * f(x[8]) + -5 * f(x[9] ) + O(h^9)f^(10)(z) err=0, h=0.125, deriv=74.3377, x=1.375, pwr=8 err=0, h=0.125, deriv=114.991, x=1.375, pwr=9 err=3.21865e-05, h=0.125, deriv=175.681, x=1.375, pwr=10 err=6.28643e-08, h=0.0625, deriv=46.9572, x=1.1875, pwr=10 err=0.000547171, h=0.125, deriv=265.717, x=1.375, pwr=11 err=8.801e-07, h=0.0625, deriv=61.3379, x=1.1875, pwr=11 computing order=1, npoints=10, at term=4 f^(1)(x[4])=(1/2520h^1)( 5 * f(x[0]) + -60 * f(x[1]) + 360 * f(x[2]) + -1680 * f(x[3]) + -504 * f(x[4]) + 2520 * f(x[5]) + -840 * f(x[6]) + 240 * f(x[7]) + -45 * f(x[8]) + 4 * f(x[9] ) + O(h^9)f^(10)(z) err=0, h=0.125, deriv=136.688, x=1.5, pwr=8 err=0, h=0.125, deriv=230.66, x=1.5, pwr=9 err=-2.14577e-05, h=0.125, deriv=384.434, x=1.5, pwr=10 err=-4.19095e-08, h=0.0625, deriv=74.5058, x=1.25, pwr=10 err=-0.000367463, h=0.125, deriv=634.315, x=1.5, pwr=11 err=-5.89353e-07, h=0.0625, deriv=102.445, x=1.25, pwr=11 computing order=1, npoints=10, at term=5 f^(1)(x[5])=(1/2520h^1)( -4 * f(x[0]) + 45 * f(x[1]) + -240 * f(x[2]) + 840 * f(x[3]) + -2520 * f(x[4]) + 504 * f(x[5]) + 1680 * f(x[6]) + -360 * f(x[7]) + 60 * f(x[8]) + -5 * f(x[9] ) + O(h^9)f^(10)(z) err=0, h=0.125, deriv=239.367, x=1.625, pwr=8 err=0, h=0.125, deriv=437.592, x=1.625, pwr=9 err=2.14577e-05, h=0.125, deriv=790.097, x=1.625, pwr=10 err=4.19095e-08, h=0.0625, deriv=115.583, x=1.3125, pwr=10 err=0.000370145, h=0.125, deriv=1412.3, x=1.625, pwr=11 err=5.91972e-07, h=0.0625, deriv=166.873, x=1.3125, pwr=11 computing order=1, npoints=10, at term=6 f^(1)(x[6])=(1/2520h^1)( 5 * f(x[0]) + -54 * f(x[1]) + 270 * f(x[2]) + -840 * f(x[3]) + 1890 * f(x[4]) + -3780 * f(x[5]) + 1554 * f(x[6]) + 1080 * f(x[7]) + -135 * f(x[8]) + 10 * f(x[9] ) + O(h^9)f^(10)(z) err=0, h=0.125, deriv=402.121, x=1.75, pwr=8 err=0, h=0.125, deriv=791.675, x=1.75, pwr=9 err=-3.21865e-05, h=0.125, deriv=1539.37, x=1.75, pwr=10 err=-6.28643e-08, h=0.0625, deriv=175.681, x=1.375, pwr=10 err=-0.000559241, h=0.125, deriv=2963.28, x=1.75, pwr=11 err=-8.91887e-07, h=0.0625, deriv=265.717, x=1.375, pwr=11 computing order=1, npoints=10, at term=7 f^(1)(x[7])=(1/2520h^1)( -10 * f(x[0]) + 105 * f(x[1]) + -504 * f(x[2]) + 1470 * f(x[3]) + -2940 * f(x[4]) + 4410 * f(x[5]) + -5880 * f(x[6]) + 2754 * f(x[7]) + 630 * f(x[8]) + -35 * f(x[9] ) + O(h^9)f^(10)(z) err=0, h=0.125, deriv=651.777, x=1.875, pwr=8 err=0, h=0.125, deriv=1374.84, x=1.875, pwr=9 err=7.51019e-05, h=0.125, deriv=2864.25, x=1.875, pwr=10 err=1.46683e-07, h=0.0625, deriv=262.102, x=1.4375, pwr=10 err=0.00131428, h=0.125, deriv=5907.52, x=1.875, pwr=11 err=2.09024e-06, h=0.0625, deriv=414.449, x=1.4375, pwr=11 computing order=1, npoints=10, at term=8 f^(1)(x[8])=(1/2520h^1)( 35 * f(x[0]) + -360 * f(x[1]) + 1680 * f(x[2]) + -4704 * f(x[3]) + 8820 * f(x[4]) + -11760 * f(x[5]) + 11760 * f(x[6]) + -10080 * f(x[7]) + 4329 * f(x[8]) + 280 * f(x[9] ) + O(h^9)f^(10)(z) err=0, h=0.125, deriv=1024, x=2, pwr=8 err=0, h=0.125, deriv=2304, x=2, pwr=9 err=-0.000300407, h=0.125, deriv=5120, x=2, pwr=10 err=-5.86733e-07, h=0.0625, deriv=384.434, x=1.5, pwr=10 err=-0.00529468, h=0.125, deriv=11264, x=2, pwr=11 err=-8.39762e-06, h=0.0625, deriv=634.315, x=1.5, pwr=11 computing order=1, npoints=10, at term=9 f^(1)(x[9])=(1/2520h^1)( -280 * f(x[0]) + 2835 * f(x[1]) + -12960 * f(x[2]) + 35280 * f(x[3]) + -63504 * f(x[4]) + 79380 * f(x[5]) + -70560 * f(x[6]) + 45360 * f(x[7]) + -22680 * f(x[8]) + 7129 * f(x[9] ) + O(h^9)f^(10)(z) err=0, h=0.125, deriv=1565.32, x=2.125, pwr=8 err=0, h=0.125, deriv=3742.09, x=2.125, pwr=9 err=0.00270367, h=0.125, deriv=8835.49, x=2.125, pwr=10 err=5.2806e-06, h=0.0625, deriv=555.112, x=1.5625, pwr=10 err=0.0479901, h=0.125, deriv=20652.9, x=2.125, pwr=11 err=7.59086e-05, h=0.0625, deriv=954.098, x=1.5625, pwr=11 computing order=2, npoints=3, at term=0 f^(2)(x[0])=(1/1h^2)( 1 * f(x[0]) + -2 * f(x[1]) + 1 * f(x[2] ) + O(h^1)f^(3)(z) err=0, h=0.125, deriv=0, x=1, pwr=1 err=0, h=0.125, deriv=2, x=1, pwr=2 err=-0.75, h=0.125, deriv=6, x=1, pwr=3 err=-0.375, h=0.0625, deriv=6, x=1, pwr=3 err=-3.21875, h=0.125, deriv=12, x=1, pwr=4 err=-1.55469, h=0.0625, deriv=12, x=1, pwr=4 computing order=2, npoints=3, at term=1 f^(2)(x[1])=(1/1h^2)( 1 * f(x[0]) + -2 * f(x[1]) + 1 * f(x[2] ) + O(h^1)f^(3)(z) err=0, h=0.125, deriv=0, x=1.125, pwr=1 err=0, h=0.125, deriv=2, x=1.125, pwr=2 err=0, h=0.125, deriv=6.75, x=1.125, pwr=3 err=0, h=0.0625, deriv=6.375, x=1.0625, pwr=3 err=-0.03125, h=0.125, deriv=15.1875, x=1.125, pwr=4 err=-0.0078125, h=0.0625, deriv=13.5469, x=1.0625, pwr=4 computing order=2, npoints=3, at term=2 f^(2)(x[2])=(1/1h^2)( 1 * f(x[0]) + -2 * f(x[1]) + 1 * f(x[2] ) + O(h^1)f^(3)(z) err=0, h=0.125, deriv=0, x=1.25, pwr=1 err=0, h=0.125, deriv=2, x=1.25, pwr=2 err=0.75, h=0.125, deriv=7.5, x=1.25, pwr=3 err=0.375, h=0.0625, deriv=6.75, x=1.125, pwr=3 err=3.53125, h=0.125, deriv=18.75, x=1.25, pwr=4 err=1.63281, h=0.0625, deriv=15.1875, x=1.125, pwr=4 computing order=2, npoints=4, at term=0 f^(2)(x[0])=(1/1h^2)( 2 * f(x[0]) + -5 * f(x[1]) + 4 * f(x[2]) + -1 * f(x[3] ) + O(h^2)f^(4)(z) err=0, h=0.125, deriv=2, x=1, pwr=2 err=0, h=0.125, deriv=6, x=1, pwr=3 err=0.34375, h=0.125, deriv=12, x=1, pwr=4 err=0.0859375, h=0.0625, deriv=12, x=1, pwr=4 err=1.95312, h=0.125, deriv=20, x=1, pwr=5 err=0.458984, h=0.0625, deriv=20, x=1, pwr=5 computing order=2, npoints=4, at term=1 f^(2)(x[1])=(1/1h^2)( 1 * f(x[0]) + -2 * f(x[1]) + 1 * f(x[2]) + 0 * f(x[3] ) + O(h^2)f^(4)(z) err=0, h=0.125, deriv=2, x=1.125, pwr=2 err=0, h=0.125, deriv=6.75, x=1.125, pwr=3 err=-0.03125, h=0.125, deriv=15.1875, x=1.125, pwr=4 err=-0.0078125, h=0.0625, deriv=13.5469, x=1.0625, pwr=4 err=-0.175781, h=0.125, deriv=28.4766, x=1.125, pwr=5 err=-0.0415039, h=0.0625, deriv=23.9893, x=1.0625, pwr=5 computing order=2, npoints=4, at term=2 f^(2)(x[2])=(1/1h^2)( 0 * f(x[0]) + 1 * f(x[1]) + -2 * f(x[2]) + 1 * f(x[3] ) + O(h^2)f^(4)(z) err=0, h=0.125, deriv=2, x=1.25, pwr=2 err=0, h=0.125, deriv=7.5, x=1.25, pwr=3 err=-0.03125, h=0.125, deriv=18.75, x=1.25, pwr=4 err=-0.0078125, h=0.0625, deriv=15.1875, x=1.125, pwr=4 err=-0.195312, h=0.125, deriv=39.0625, x=1.25, pwr=5 err=-0.0439453, h=0.0625, deriv=28.4766, x=1.125, pwr=5 computing order=2, npoints=4, at term=3 f^(2)(x[3])=(1/1h^2)( -1 * f(x[0]) + 4 * f(x[1]) + -5 * f(x[2]) + 2 * f(x[3] ) + O(h^2)f^(4)(z) err=0, h=0.125, deriv=2, x=1.375, pwr=2 err=0, h=0.125, deriv=8.25, x=1.375, pwr=3 err=0.34375, h=0.125, deriv=22.6875, x=1.375, pwr=4 err=0.0859375, h=0.0625, deriv=16.9219, x=1.1875, pwr=4 err=2.12891, h=0.125, deriv=51.9922, x=1.375, pwr=5 err=0.480957, h=0.0625, deriv=33.4912, x=1.1875, pwr=5 computing order=2, npoints=5, at term=0 f^(2)(x[0])=(1/12h^2)( 35 * f(x[0]) + -104 * f(x[1]) + 114 * f(x[2]) + -56 * f(x[3]) + 11 * f(x[4] ) + O(h^3)f^(5)(z) err=0, h=0.125, deriv=6, x=1, pwr=3 err=0, h=0.125, deriv=12, x=1, pwr=4 err=-0.195312, h=0.125, deriv=20, x=1, pwr=5 err=-0.0244141, h=0.0625, deriv=20, x=1, pwr=5 err=-1.4043, h=0.125, deriv=30, x=1, pwr=6 err=-0.161011, h=0.0625, deriv=30, x=1, pwr=6 computing order=2, npoints=5, at term=1 f^(2)(x[1])=(1/12h^2)( 11 * f(x[0]) + -20 * f(x[1]) + 6 * f(x[2]) + 4 * f(x[3]) + -1 * f(x[4] ) + O(h^3)f^(5)(z) err=0, h=0.125, deriv=6.75, x=1.125, pwr=3 err=0, h=0.125, deriv=15.1875, x=1.125, pwr=4 err=0.0195312, h=0.125, deriv=28.4766, x=1.125, pwr=5 err=0.00244141, h=0.0625, deriv=23.9893, x=1.0625, pwr=5 err=0.141113, h=0.125, deriv=48.0542, x=1.125, pwr=6 err=0.0161438, h=0.0625, deriv=38.2329, x=1.0625, pwr=6 computing order=2, npoints=5, at term=2 f^(2)(x[2])=(1/12h^2)( -1 * f(x[0]) + 16 * f(x[1]) + -30 * f(x[2]) + 16 * f(x[3]) + -1 * f(x[4] ) + O(h^3)f^(5)(z) err=0, h=0.125, deriv=7.5, x=1.25, pwr=3 err=0, h=0.125, deriv=18.75, x=1.25, pwr=4 err=0, h=0.125, deriv=39.0625, x=1.25, pwr=5 err=0, h=0.0625, deriv=28.4766, x=1.125, pwr=5 err=0.00195312, h=0.125, deriv=73.2422, x=1.25, pwr=6 err=0.00012207, h=0.0625, deriv=48.0542, x=1.125, pwr=6 computing order=2, npoints=5, at term=3 f^(2)(x[3])=(1/12h^2)( -1 * f(x[0]) + 4 * f(x[1]) + 6 * f(x[2]) + -20 * f(x[3]) + 11 * f(x[4] ) + O(h^3)f^(5)(z) err=0, h=0.125, deriv=8.25, x=1.375, pwr=3 err=0, h=0.125, deriv=22.6875, x=1.375, pwr=4 err=-0.0195312, h=0.125, deriv=51.9922, x=1.375, pwr=5 err=-0.00244141, h=0.0625, deriv=33.4912, x=1.1875, pwr=5 err=-0.151855, h=0.125, deriv=107.234, x=1.375, pwr=6 err=-0.0168152, h=0.0625, deriv=59.6562, x=1.1875, pwr=6 computing order=2, npoints=5, at term=4 f^(2)(x[4])=(1/12h^2)( 11 * f(x[0]) + -56 * f(x[1]) + 114 * f(x[2]) + -104 * f(x[3]) + 35 * f(x[4] ) + O(h^3)f^(5)(z) err=0, h=0.125, deriv=9, x=1.5, pwr=3 err=0, h=0.125, deriv=27, x=1.5, pwr=4 err=0.195312, h=0.125, deriv=67.5, x=1.5, pwr=5 err=0.0244141, h=0.0625, deriv=39.0625, x=1.25, pwr=5 err=1.52539, h=0.125, deriv=151.875, x=1.5, pwr=6 err=0.168579, h=0.0625, deriv=73.2422, x=1.25, pwr=6 computing order=2, npoints=6, at term=0 f^(2)(x[0])=(1/12h^2)( 45 * f(x[0]) + -154 * f(x[1]) + 214 * f(x[2]) + -156 * f(x[3]) + 61 * f(x[4]) + -10 * f(x[5] ) + O(h^4)f^(6)(z) err=0, h=0.125, deriv=12, x=1, pwr=4 err=0, h=0.125, deriv=20, x=1, pwr=5 err=0.133789, h=0.125, deriv=30, x=1, pwr=6 err=0.00836182, h=0.0625, deriv=30, x=1, pwr=6 err=1.18005, h=0.125, deriv=42, x=1, pwr=7 err=0.066143, h=0.0625, deriv=42, x=1, pwr=7 computing order=2, npoints=6, at term=1 f^(2)(x[1])=(1/12h^2)( 10 * f(x[0]) + -15 * f(x[1]) + -4 * f(x[2]) + 14 * f(x[3]) + -6 * f(x[4]) + 1 * f(x[5] ) + O(h^4)f^(6)(z) err=0, h=0.125, deriv=15.1875, x=1.125, pwr=4 err=0, h=0.125, deriv=28.4766, x=1.125, pwr=5 err=-0.0126953, h=0.125, deriv=48.0542, x=1.125, pwr=6 err=-0.000793457, h=0.0625, deriv=38.2329, x=1.0625, pwr=6 err=-0.112793, h=0.125, deriv=75.6854, x=1.125, pwr=7 err=-0.00630188, h=0.0625, deriv=56.8714, x=1.0625, pwr=7 computing order=2, npoints=6, at term=2 f^(2)(x[2])=(1/12h^2)( -1 * f(x[0]) + 16 * f(x[1]) + -30 * f(x[2]) + 16 * f(x[3]) + -1 * f(x[4]) + 0 * f(x[5] ) + O(h^4)f^(6)(z) err=0, h=0.125, deriv=18.75, x=1.25, pwr=4 err=0, h=0.125, deriv=39.0625, x=1.25, pwr=5 err=0.00195312, h=0.125, deriv=73.2422, x=1.25, pwr=6 err=0.00012207, h=0.0625, deriv=48.0542, x=1.125, pwr=6 err=0.0170898, h=0.125, deriv=128.174, x=1.25, pwr=7 err=0.000961304, h=0.0625, deriv=75.6854, x=1.125, pwr=7 computing order=2, npoints=6, at term=3 f^(2)(x[3])=(1/12h^2)( 0 * f(x[0]) + -1 * f(x[1]) + 16 * f(x[2]) + -30 * f(x[3]) + 16 * f(x[4]) + -1 * f(x[5] ) + O(h^4)f^(6)(z) err=0, h=0.125, deriv=22.6875, x=1.375, pwr=4 err=0, h=0.125, deriv=51.9922, x=1.375, pwr=5 err=0.00195312, h=0.125, deriv=107.234, x=1.375, pwr=6 err=0.00012207, h=0.0625, deriv=59.6562, x=1.1875, pwr=6 err=0.0187988, h=0.125, deriv=206.425, x=1.375, pwr=7 err=0.00101471, h=0.0625, deriv=99.1785, x=1.1875, pwr=7 computing order=2, npoints=6, at term=4 f^(2)(x[4])=(1/12h^2)( 1 * f(x[0]) + -6 * f(x[1]) + 14 * f(x[2]) + -4 * f(x[3]) + -15 * f(x[4]) + 10 * f(x[5] ) + O(h^4)f^(6)(z) err=0, h=0.125, deriv=27, x=1.5, pwr=4 err=0, h=0.125, deriv=67.5, x=1.5, pwr=5 err=-0.0126953, h=0.125, deriv=151.875, x=1.5, pwr=6 err=-0.000793457, h=0.0625, deriv=73.2422, x=1.25, pwr=6 err=-0.120483, h=0.125, deriv=318.938, x=1.5, pwr=7 err=-0.00654221, h=0.0625, deriv=128.174, x=1.25, pwr=7 computing order=2, npoints=6, at term=5 f^(2)(x[5])=(1/12h^2)( -10 * f(x[0]) + 61 * f(x[1]) + -156 * f(x[2]) + 214 * f(x[3]) + -154 * f(x[4]) + 45 * f(x[5] ) + O(h^4)f^(6)(z) err=0, h=0.125, deriv=31.6875, x=1.625, pwr=4 err=0, h=0.125, deriv=85.8203, x=1.625, pwr=5 err=0.133789, h=0.125, deriv=209.187, x=1.625, pwr=6 err=0.00836182, h=0.0625, deriv=89.0263, x=1.3125, pwr=6 err=1.27832, h=0.125, deriv=475.9, x=1.625, pwr=7 err=0.0692139, h=0.0625, deriv=163.586, x=1.3125, pwr=7 computing order=2, npoints=7, at term=0 f^(2)(x[0])=(1/180h^2)( 812 * f(x[0]) + -3132 * f(x[1]) + 5265 * f(x[2]) + -5080 * f(x[3]) + 2970 * f(x[4]) + -972 * f(x[5]) + 137 * f(x[6] ) + O(h^5)f^(7)(z) err=0, h=0.125, deriv=20, x=1, pwr=5 err=0, h=0.125, deriv=30, x=1, pwr=6 err=-0.107666, h=0.125, deriv=42, x=1, pwr=7 err=-0.00336456, h=0.0625, deriv=42, x=1, pwr=7 err=-1.13846, h=0.125, deriv=56, x=1, pwr=8 err=-0.0312467, h=0.0625, deriv=56, x=1, pwr=8 computing order=2, npoints=7, at term=1 f^(2)(x[1])=(1/180h^2)( 137 * f(x[0]) + -147 * f(x[1]) + -255 * f(x[2]) + 470 * f(x[3]) + -285 * f(x[4]) + 93 * f(x[5]) + -13 * f(x[6] ) + O(h^5)f^(7)(z) err=0, h=0.125, deriv=28.4766, x=1.125, pwr=5 err=0, h=0.125, deriv=48.0542, x=1.125, pwr=6 err=0.00939941, h=0.125, deriv=75.6854, x=1.125, pwr=7 err=0.000293732, h=0.0625, deriv=56.8714, x=1.0625, pwr=7 err=0.100128, h=0.125, deriv=113.528, x=1.125, pwr=8 err=0.00273943, h=0.0625, deriv=80.5678, x=1.0625, pwr=8 computing order=2, npoints=7, at term=2 f^(2)(x[2])=(1/180h^2)( -13 * f(x[0]) + 228 * f(x[1]) + -420 * f(x[2]) + 200 * f(x[3]) + 15 * f(x[4]) + -12 * f(x[5]) + 2 * f(x[6] ) + O(h^5)f^(7)(z) err=0, h=0.125, deriv=39.0625, x=1.25, pwr=5 err=0, h=0.125, deriv=73.2422, x=1.25, pwr=6 err=-0.00170898, h=0.125, deriv=128.174, x=1.25, pwr=7 err=-5.34058e-05, h=0.0625, deriv=75.6854, x=1.125, pwr=7 err=-0.018219, h=0.125, deriv=213.623, x=1.25, pwr=8 err=-0.000498295, h=0.0625, deriv=113.528, x=1.125, pwr=8 computing order=2, npoints=7, at term=3 f^(2)(x[3])=(1/180h^2)( 2 * f(x[0]) + -27 * f(x[1]) + 270 * f(x[2]) + -490 * f(x[3]) + 270 * f(x[4]) + -27 * f(x[5]) + 2 * f(x[6] ) + O(h^5)f^(7)(z) err=0, h=0.125, deriv=51.9922, x=1.375, pwr=5 err=0, h=0.125, deriv=107.234, x=1.375, pwr=6 err=0, h=0.125, deriv=206.425, x=1.375, pwr=7 err=0, h=0.0625, deriv=99.1785, x=1.1875, pwr=7 err=-0.000274658, h=0.125, deriv=378.446, x=1.375, pwr=8 err=-4.29153e-06, h=0.0625, deriv=157.033, x=1.1875, pwr=8 computing order=2, npoints=7, at term=4 f^(2)(x[4])=(1/180h^2)( 2 * f(x[0]) + -12 * f(x[1]) + 15 * f(x[2]) + 200 * f(x[3]) + -420 * f(x[4]) + 228 * f(x[5]) + -13 * f(x[6] ) + O(h^5)f^(7)(z) err=0, h=0.125, deriv=67.5, x=1.5, pwr=5 err=0, h=0.125, deriv=151.875, x=1.5, pwr=6 err=0.00170898, h=0.125, deriv=318.938, x=1.5, pwr=7 err=5.34058e-05, h=0.0625, deriv=128.174, x=1.25, pwr=7 err=0.0193787, h=0.125, deriv=637.875, x=1.5, pwr=8 err=0.000516415, h=0.0625, deriv=213.623, x=1.25, pwr=8 computing order=2, npoints=7, at term=5 f^(2)(x[5])=(1/180h^2)( -13 * f(x[0]) + 93 * f(x[1]) + -285 * f(x[2]) + 470 * f(x[3]) + -255 * f(x[4]) + -147 * f(x[5]) + 137 * f(x[6] ) + O(h^5)f^(7)(z) err=0, h=0.125, deriv=85.8203, x=1.625, pwr=5 err=0, h=0.125, deriv=209.187, x=1.625, pwr=6 err=-0.00939941, h=0.125, deriv=475.9, x=1.625, pwr=7 err=-0.000293732, h=0.0625, deriv=163.586, x=1.3125, pwr=7 err=-0.106659, h=0.125, deriv=1031.12, x=1.625, pwr=8 err=-0.00284147, h=0.0625, deriv=286.275, x=1.3125, pwr=8 computing order=2, npoints=7, at term=6 f^(2)(x[6])=(1/180h^2)( 137 * f(x[0]) + -972 * f(x[1]) + 2970 * f(x[2]) + -5080 * f(x[3]) + 5265 * f(x[4]) + -3132 * f(x[5]) + 812 * f(x[6] ) + O(h^5)f^(7)(z) err=0, h=0.125, deriv=107.188, x=1.75, pwr=5 err=0, h=0.125, deriv=281.367, x=1.75, pwr=6 err=0.107666, h=0.125, deriv=689.35, x=1.75, pwr=7 err=0.00336456, h=0.0625, deriv=206.425, x=1.375, pwr=7 err=1.23019, h=0.125, deriv=1608.48, x=1.75, pwr=8 err=0.03268, h=0.0625, deriv=378.446, x=1.375, pwr=8 computing order=2, npoints=8, at term=0 f^(2)(x[0])=(1/180h^2)( 938 * f(x[0]) + -4014 * f(x[1]) + 7911 * f(x[2]) + -9490 * f(x[3]) + 7380 * f(x[4]) + -3618 * f(x[5]) + 1019 * f(x[6]) + -126 * f(x[7] ) + O(h^6)f^(8)(z) err=0, h=0.125, deriv=30, x=1, pwr=6 err=0, h=0.125, deriv=42, x=1, pwr=7 err=0.0997009, h=0.125, deriv=56, x=1, pwr=8 err=0.00155783, h=0.0625, deriv=56, x=1, pwr=8 err=1.24146, h=0.125, deriv=72, x=1, pwr=9 err=0.0167091, h=0.0625, deriv=72, x=1, pwr=9 computing order=2, npoints=8, at term=1 f^(2)(x[1])=(1/180h^2)( 126 * f(x[0]) + -70 * f(x[1]) + -486 * f(x[2]) + 855 * f(x[3]) + -670 * f(x[4]) + 324 * f(x[5]) + -90 * f(x[6]) + 11 * f(x[7] ) + O(h^6)f^(8)(z) err=0, h=0.125, deriv=48.0542, x=1.125, pwr=6 err=0, h=0.125, deriv=75.6854, x=1.125, pwr=7 err=-0.00796509, h=0.125, deriv=113.528, x=1.125, pwr=8 err=-0.000124454, h=0.0625, deriv=80.5678, x=1.0625, pwr=8 err=-0.0998726, h=0.125, deriv=164.21, x=1.125, pwr=9 err=-0.0013403, h=0.0625, deriv=110.061, x=1.0625, pwr=9 computing order=2, npoints=8, at term=2 f^(2)(x[2])=(1/180h^2)( -11 * f(x[0]) + 214 * f(x[1]) + -378 * f(x[2]) + 130 * f(x[3]) + 85 * f(x[4]) + -54 * f(x[5]) + 16 * f(x[6]) + -2 * f(x[7] ) + O(h^6)f^(8)(z) err=0, h=0.125, deriv=73.2422, x=1.25, pwr=6 err=0, h=0.125, deriv=128.174, x=1.25, pwr=7 err=0.00143433, h=0.125, deriv=213.623, x=1.25, pwr=8 err=2.24113e-05, h=0.0625, deriv=113.528, x=1.125, pwr=8 err=0.0180588, h=0.125, deriv=343.323, x=1.25, pwr=9 err=0.000241935, h=0.0625, deriv=164.21, x=1.125, pwr=9 computing order=2, npoints=8, at term=3 f^(2)(x[3])=(1/180h^2)( 2 * f(x[0]) + -27 * f(x[1]) + 270 * f(x[2]) + -490 * f(x[3]) + 270 * f(x[4]) + -27 * f(x[5]) + 2 * f(x[6]) + 0 * f(x[7] ) + O(h^6)f^(8)(z) err=0, h=0.125, deriv=107.234, x=1.375, pwr=6 err=0, h=0.125, deriv=206.425, x=1.375, pwr=7 err=-0.000274658, h=0.125, deriv=378.446, x=1.375, pwr=8 err=-4.29153e-06, h=0.0625, deriv=157.033, x=1.1875, pwr=8 err=-0.0033989, h=0.125, deriv=669.039, x=1.375, pwr=9 err=-4.58658e-05, h=0.0625, deriv=239.755, x=1.1875, pwr=9 computing order=2, npoints=8, at term=4 f^(2)(x[4])=(1/180h^2)( 0 * f(x[0]) + 2 * f(x[1]) + -27 * f(x[2]) + 270 * f(x[3]) + -490 * f(x[4]) + 270 * f(x[5]) + -27 * f(x[6]) + 2 * f(x[7] ) + O(h^6)f^(8)(z) err=0, h=0.125, deriv=151.875, x=1.5, pwr=6 err=0, h=0.125, deriv=318.938, x=1.5, pwr=7 err=-0.000274658, h=0.125, deriv=637.875, x=1.5, pwr=8 err=-4.29153e-06, h=0.0625, deriv=213.623, x=1.25, pwr=8 err=-0.00370789, h=0.125, deriv=1230.19, x=1.5, pwr=9 err=-4.82798e-05, h=0.0625, deriv=343.323, x=1.25, pwr=9 computing order=2, npoints=8, at term=5 f^(2)(x[5])=(1/180h^2)( -2 * f(x[0]) + 16 * f(x[1]) + -54 * f(x[2]) + 85 * f(x[3]) + 130 * f(x[4]) + -378 * f(x[5]) + 214 * f(x[6]) + -11 * f(x[7] ) + O(h^6)f^(8)(z) err=0, h=0.125, deriv=209.187, x=1.625, pwr=6 err=0, h=0.125, deriv=475.9, x=1.625, pwr=7 err=0.00143433, h=0.125, deriv=1031.12, x=1.625, pwr=8 err=2.24113e-05, h=0.0625, deriv=286.275, x=1.3125, pwr=8 err=0.0190544, h=0.125, deriv=2154.3, x=1.625, pwr=9 err=0.000249714, h=0.0625, deriv=483.09, x=1.3125, pwr=9 computing order=2, npoints=8, at term=6 f^(2)(x[6])=(1/180h^2)( 11 * f(x[0]) + -90 * f(x[1]) + 324 * f(x[2]) + -670 * f(x[3]) + 855 * f(x[4]) + -486 * f(x[5]) + -70 * f(x[6]) + 126 * f(x[7] ) + O(h^6)f^(8)(z) err=0, h=0.125, deriv=281.367, x=1.75, pwr=6 err=0, h=0.125, deriv=689.35, x=1.75, pwr=7 err=-0.00796509, h=0.125, deriv=1608.48, x=1.75, pwr=8 err=-0.000124454, h=0.0625, deriv=378.446, x=1.375, pwr=8 err=-0.106224, h=0.125, deriv=3619.09, x=1.75, pwr=9 err=-0.00138992, h=0.0625, deriv=669.039, x=1.375, pwr=9 computing order=2, npoints=8, at term=7 f^(2)(x[7])=(1/180h^2)( -126 * f(x[0]) + 1019 * f(x[1]) + -3618 * f(x[2]) + 7380 * f(x[3]) + -9490 * f(x[4]) + 7911 * f(x[5]) + -4014 * f(x[6]) + 938 * f(x[7] ) + O(h^6)f^(8)(z) err=0, h=0.125, deriv=370.789, x=1.875, pwr=6 err=0, h=0.125, deriv=973.32, x=1.875, pwr=7 err=0.0997009, h=0.125, deriv=2433.3, x=1.875, pwr=8 err=0.00155783, h=0.0625, deriv=494.123, x=1.4375, pwr=8 err=1.33831, h=0.125, deriv=5865.99, x=1.875, pwr=9 err=0.0174657, h=0.0625, deriv=913.245, x=1.4375, pwr=9 computing order=3, npoints=4, at term=0 f^(3)(x[0])=(1/1h^3)( -1 * f(x[0]) + 3 * f(x[1]) + -3 * f(x[2]) + 1 * f(x[3] ) + O(h^1)f^(4)(z) err=0, h=0.125, deriv=0, x=1, pwr=2 err=0, h=0.125, deriv=6, x=1, pwr=3 err=-4.5, h=0.125, deriv=24, x=1, pwr=4 err=-2.25, h=0.0625, deriv=24, x=1, pwr=4 err=-24.8438, h=0.125, deriv=60, x=1, pwr=5 err=-11.8359, h=0.0625, deriv=60, x=1, pwr=5 computing order=3, npoints=4, at term=1 f^(3)(x[1])=(1/1h^3)( -1 * f(x[0]) + 3 * f(x[1]) + -3 * f(x[2]) + 1 * f(x[3] ) + O(h^1)f^(4)(z) err=0, h=0.125, deriv=0, x=1.125, pwr=2 err=0, h=0.125, deriv=6, x=1.125, pwr=3 err=-1.5, h=0.125, deriv=27, x=1.125, pwr=4 err=-0.75, h=0.0625, deriv=25.5, x=1.0625, pwr=4 err=-8.90625, h=0.125, deriv=75.9375, x=1.125, pwr=5 err=-4.10156, h=0.0625, deriv=67.7344, x=1.0625, pwr=5 computing order=3, npoints=4, at term=2 f^(3)(x[2])=(1/1h^3)( -1 * f(x[0]) + 3 * f(x[1]) + -3 * f(x[2]) + 1 * f(x[3] ) + O(h^1)f^(4)(z) err=0, h=0.125, deriv=0, x=1.25, pwr=2 err=0, h=0.125, deriv=6, x=1.25, pwr=3 err=1.5, h=0.125, deriv=30, x=1.25, pwr=4 err=0.75, h=0.0625, deriv=27, x=1.125, pwr=4 err=8.90625, h=0.125, deriv=93.75, x=1.25, pwr=5 err=4.10156, h=0.0625, deriv=75.9375, x=1.125, pwr=5 computing order=3, npoints=4, at term=3 f^(3)(x[3])=(1/1h^3)( -1 * f(x[0]) + 3 * f(x[1]) + -3 * f(x[2]) + 1 * f(x[3] ) + O(h^1)f^(4)(z) err=0, h=0.125, deriv=0, x=1.375, pwr=2 err=0, h=0.125, deriv=6, x=1.375, pwr=3 err=4.5, h=0.125, deriv=33, x=1.375, pwr=4 err=2.25, h=0.0625, deriv=28.5, x=1.1875, pwr=4 err=28.5938, h=0.125, deriv=113.438, x=1.375, pwr=5 err=12.7734, h=0.0625, deriv=84.6094, x=1.1875, pwr=5 computing order=3, npoints=5, at term=0 f^(3)(x[0])=(1/2h^3)( -5 * f(x[0]) + 18 * f(x[1]) + -24 * f(x[2]) + 14 * f(x[3]) + -3 * f(x[4] ) + O(h^2)f^(5)(z) err=0, h=0.125, deriv=6, x=1, pwr=3 err=0, h=0.125, deriv=24, x=1, pwr=4 err=3.28125, h=0.125, deriv=60, x=1, pwr=5 err=0.820312, h=0.0625, deriv=60, x=1, pwr=5 err=23.2031, h=0.125, deriv=120, x=1, pwr=6 err=5.36133, h=0.0625, deriv=120, x=1, pwr=6 computing order=3, npoints=5, at term=1 f^(3)(x[1])=(1/2h^3)( -3 * f(x[0]) + 10 * f(x[1]) + -12 * f(x[2]) + 6 * f(x[3]) + -1 * f(x[4] ) + O(h^2)f^(5)(z) err=0, h=0.125, deriv=6, x=1.125, pwr=3 err=0, h=0.125, deriv=27, x=1.125, pwr=4 err=0.46875, h=0.125, deriv=75.9375, x=1.125, pwr=5 err=0.117188, h=0.0625, deriv=67.7344, x=1.0625, pwr=5 err=3.51562, h=0.125, deriv=170.859, x=1.125, pwr=6 err=0.791016, h=0.0625, deriv=143.936, x=1.0625, pwr=6 computing order=3, npoints=5, at term=2 f^(3)(x[2])=(1/2h^3)( -1 * f(x[0]) + 2 * f(x[1]) + 0 * f(x[2]) + -2 * f(x[3]) + 1 * f(x[4] ) + O(h^2)f^(5)(z) err=0, h=0.125, deriv=6, x=1.25, pwr=3 err=0, h=0.125, deriv=30, x=1.25, pwr=4 err=-0.46875, h=0.125, deriv=93.75, x=1.25, pwr=5 err=-0.117188, h=0.0625, deriv=75.9375, x=1.125, pwr=5 err=-3.51562, h=0.125, deriv=234.375, x=1.25, pwr=6 err=-0.791016, h=0.0625, deriv=170.859, x=1.125, pwr=6 computing order=3, npoints=5, at term=3 f^(3)(x[3])=(1/2h^3)( 1 * f(x[0]) + -6 * f(x[1]) + 12 * f(x[2]) + -10 * f(x[3]) + 3 * f(x[4] ) + O(h^2)f^(5)(z) err=0, h=0.125, deriv=6, x=1.375, pwr=3 err=0, h=0.125, deriv=33, x=1.375, pwr=4 err=0.46875, h=0.125, deriv=113.438, x=1.375, pwr=5 err=0.117188, h=0.0625, deriv=84.6094, x=1.1875, pwr=5 err=3.51562, h=0.125, deriv=311.953, x=1.375, pwr=6 err=0.791016, h=0.0625, deriv=200.947, x=1.1875, pwr=6 computing order=3, npoints=5, at term=4 f^(3)(x[4])=(1/2h^3)( 3 * f(x[0]) + -14 * f(x[1]) + 24 * f(x[2]) + -18 * f(x[3]) + 5 * f(x[4] ) + O(h^2)f^(5)(z) err=0, h=0.125, deriv=6, x=1.5, pwr=3 err=0, h=0.125, deriv=36, x=1.5, pwr=4 err=3.28125, h=0.125, deriv=135, x=1.5, pwr=5 err=0.820312, h=0.0625, deriv=93.75, x=1.25, pwr=5 err=26.0156, h=0.125, deriv=405, x=1.5, pwr=6 err=5.71289, h=0.0625, deriv=234.375, x=1.25, pwr=6 computing order=3, npoints=6, at term=0 f^(3)(x[0])=(1/4h^3)( -17 * f(x[0]) + 71 * f(x[1]) + -118 * f(x[2]) + 98 * f(x[3]) + -41 * f(x[4]) + 7 * f(x[5] ) + O(h^3)f^(6)(z) err=0, h=0.125, deriv=24, x=1, pwr=4 err=0, h=0.125, deriv=60, x=1, pwr=5 err=-2.63672, h=0.125, deriv=120, x=1, pwr=6 err=-0.32959, h=0.0625, deriv=120, x=1, pwr=6 err=-22.9995, h=0.125, deriv=210, x=1, pwr=7 err=-2.59103, h=0.0625, deriv=210, x=1, pwr=7 computing order=3, npoints=6, at term=1 f^(3)(x[1])=(1/4h^3)( -7 * f(x[0]) + 25 * f(x[1]) + -34 * f(x[2]) + 22 * f(x[3]) + -7 * f(x[4]) + 1 * f(x[5] ) + O(h^3)f^(6)(z) err=0, h=0.125, deriv=27, x=1.125, pwr=4 err=0, h=0.125, deriv=75.9375, x=1.125, pwr=5 err=-0.175781, h=0.125, deriv=170.859, x=1.125, pwr=6 err=-0.0219727, h=0.0625, deriv=143.936, x=1.0625, pwr=6 err=-1.62012, h=0.125, deriv=336.379, x=1.125, pwr=7 err=-0.178162, h=0.0625, deriv=267.63, x=1.0625, pwr=7 computing order=3, npoints=6, at term=2 f^(3)(x[2])=(1/4h^3)( -1 * f(x[0]) + -1 * f(x[1]) + 10 * f(x[2]) + -14 * f(x[3]) + 7 * f(x[4]) + -1 * f(x[5] ) + O(h^3)f^(6)(z) err=0, h=0.125, deriv=30, x=1.25, pwr=4 err=0, h=0.125, deriv=93.75, x=1.25, pwr=5 err=0.175781, h=0.125, deriv=234.375, x=1.25, pwr=6 err=0.0219727, h=0.0625, deriv=170.859, x=1.125, pwr=6 err=1.60986, h=0.125, deriv=512.695, x=1.25, pwr=7 err=0.177521, h=0.0625, deriv=336.379, x=1.125, pwr=7 computing order=3, npoints=6, at term=3 f^(3)(x[3])=(1/4h^3)( 1 * f(x[0]) + -7 * f(x[1]) + 14 * f(x[2]) + -10 * f(x[3]) + 1 * f(x[4]) + 1 * f(x[5] ) + O(h^3)f^(6)(z) err=0, h=0.125, deriv=33, x=1.375, pwr=4 err=0, h=0.125, deriv=113.438, x=1.375, pwr=5 err=-0.175781, h=0.125, deriv=311.953, x=1.375, pwr=6 err=-0.0219727, h=0.0625, deriv=200.947, x=1.1875, pwr=6 err=-1.62012, h=0.125, deriv=750.637, x=1.375, pwr=7 err=-0.178162, h=0.0625, deriv=417.594, x=1.1875, pwr=7 computing order=3, npoints=6, at term=4 f^(3)(x[4])=(1/4h^3)( -1 * f(x[0]) + 7 * f(x[1]) + -22 * f(x[2]) + 34 * f(x[3]) + -25 * f(x[4]) + 7 * f(x[5] ) + O(h^3)f^(6)(z) err=0, h=0.125, deriv=36, x=1.5, pwr=4 err=0, h=0.125, deriv=135, x=1.5, pwr=5 err=0.175781, h=0.125, deriv=405, x=1.5, pwr=6 err=0.0219727, h=0.0625, deriv=234.375, x=1.25, pwr=6 err=1.60986, h=0.125, deriv=1063.12, x=1.5, pwr=7 err=0.177521, h=0.0625, deriv=512.695, x=1.25, pwr=7 computing order=3, npoints=6, at term=5 f^(3)(x[5])=(1/4h^3)( -7 * f(x[0]) + 41 * f(x[1]) + -98 * f(x[2]) + 118 * f(x[3]) + -71 * f(x[4]) + 17 * f(x[5] ) + O(h^3)f^(6)(z) err=0, h=0.125, deriv=39, x=1.625, pwr=4 err=0, h=0.125, deriv=158.438, x=1.625, pwr=5 err=2.63672, h=0.125, deriv=514.922, x=1.625, pwr=6 err=0.32959, h=0.0625, deriv=271.318, x=1.3125, pwr=6 err=25.4502, h=0.125, deriv=1464.31, x=1.625, pwr=7 err=2.7442, h=0.0625, deriv=623.184, x=1.3125, pwr=7 computing order=4, npoints=5, at term=0 f^(4)(x[0])=(1/1h^4)( 1 * f(x[0]) + -4 * f(x[1]) + 6 * f(x[2]) + -4 * f(x[3]) + 1 * f(x[4] ) + O(h^1)f^(5)(z) err=0, h=0.0625, deriv=0, x=1, pwr=3 err=0, h=0.0625, deriv=24, x=1, pwr=4 err=-15, h=0.0625, deriv=120, x=1, pwr=5 err=-7.5, h=0.03125, deriv=120, x=1, pwr=5 err=-96.0938, h=0.0625, deriv=360, x=1, pwr=6 err=-46.5234, h=0.03125, deriv=360, x=1, pwr=6 err=-22.8809, h=0.015625, deriv=360, x=1, pwr=6 computing order=4, npoints=5, at term=1 f^(4)(x[1])=(1/1h^4)( 1 * f(x[0]) + -4 * f(x[1]) + 6 * f(x[2]) + -4 * f(x[3]) + 1 * f(x[4] ) + O(h^1)f^(5)(z) err=0, h=0.0625, deriv=0, x=1.0625, pwr=3 err=0, h=0.0625, deriv=24, x=1.0625, pwr=4 err=-7.5, h=0.0625, deriv=127.5, x=1.0625, pwr=5 err=-3.75, h=0.03125, deriv=123.75, x=1.03125, pwr=5 err=-49.6875, h=0.0625, deriv=406.406, x=1.0625, pwr=6 err=-23.6719, h=0.03125, deriv=382.852, x=1.03125, pwr=6 err=-11.543, h=0.015625, deriv=371.338, x=1.01562, pwr=6 computing order=4, npoints=5, at term=2 f^(4)(x[2])=(1/1h^4)( 1 * f(x[0]) + -4 * f(x[1]) + 6 * f(x[2]) + -4 * f(x[3]) + 1 * f(x[4] ) + O(h^1)f^(5)(z) err=0, h=0.0625, deriv=0, x=1.125, pwr=3 err=0, h=0.0625, deriv=24, x=1.125, pwr=4 err=0, h=0.0625, deriv=135, x=1.125, pwr=5 err=0, h=0.03125, deriv=127.5, x=1.0625, pwr=5 err=-0.46875, h=0.0625, deriv=455.625, x=1.125, pwr=6 err=-0.117188, h=0.03125, deriv=406.406, x=1.0625, pwr=6 err=-0.0292969, h=0.015625, deriv=382.852, x=1.03125, pwr=6 computing order=4, npoints=5, at term=3 f^(4)(x[3])=(1/1h^4)( 1 * f(x[0]) + -4 * f(x[1]) + 6 * f(x[2]) + -4 * f(x[3]) + 1 * f(x[4] ) + O(h^1)f^(5)(z) err=0, h=0.0625, deriv=0, x=1.1875, pwr=3 err=0, h=0.0625, deriv=24, x=1.1875, pwr=4 err=7.5, h=0.0625, deriv=142.5, x=1.1875, pwr=5 err=3.75, h=0.03125, deriv=131.25, x=1.09375, pwr=5 err=51.5625, h=0.0625, deriv=507.656, x=1.1875, pwr=6 err=24.1406, h=0.03125, deriv=430.664, x=1.09375, pwr=6 err=11.6602, h=0.015625, deriv=394.541, x=1.04688, pwr=6 computing order=4, npoints=5, at term=4 f^(4)(x[4])=(1/1h^4)( 1 * f(x[0]) + -4 * f(x[1]) + 6 * f(x[2]) + -4 * f(x[3]) + 1 * f(x[4] ) + O(h^1)f^(5)(z) err=0, h=0.0625, deriv=0, x=1.25, pwr=3 err=0, h=0.0625, deriv=24, x=1.25, pwr=4 err=15, h=0.0625, deriv=150, x=1.25, pwr=5 err=7.5, h=0.03125, deriv=135, x=1.125, pwr=5 err=106.406, h=0.0625, deriv=562.5, x=1.25, pwr=6 err=49.1016, h=0.03125, deriv=455.625, x=1.125, pwr=6 err=23.5254, h=0.015625, deriv=406.406, x=1.0625, pwr=6 computing order=4, npoints=6, at term=0 f^(4)(x[0])=(1/1h^4)( 3 * f(x[0]) + -14 * f(x[1]) + 26 * f(x[2]) + -24 * f(x[3]) + 11 * f(x[4]) + -2 * f(x[5] ) + O(h^2)f^(6)(z) err=0, h=0.0625, deriv=24, x=1, pwr=4 err=0, h=0.0625, deriv=120, x=1, pwr=5 err=7.96875, h=0.0625, deriv=360, x=1, pwr=6 err=1.99219, h=0.03125, deriv=360, x=1, pwr=6 err=61.9336, h=0.0625, deriv=840, x=1, pwr=7 err=14.7144, h=0.03125, deriv=840, x=1, pwr=7 err=3.58246, h=0.015625, deriv=840, x=1, pwr=7 computing order=4, npoints=6, at term=1 f^(4)(x[1])=(1/1h^4)( 2 * f(x[0]) + -9 * f(x[1]) + 16 * f(x[2]) + -14 * f(x[3]) + 6 * f(x[4]) + -1 * f(x[5] ) + O(h^2)f^(6)(z) err=0, h=0.0625, deriv=24, x=1.0625, pwr=4 err=0, h=0.0625, deriv=127.5, x=1.0625, pwr=5 err=2.34375, h=0.0625, deriv=406.406, x=1.0625, pwr=6 err=0.585938, h=0.03125, deriv=382.852, x=1.03125, pwr=6 err=18.6621, h=0.0625, deriv=1007.55, x=1.0625, pwr=7 err=4.38354, h=0.03125, deriv=921.237, x=1.03125, pwr=7 err=1.06064, h=0.015625, deriv=879.993, x=1.01562, pwr=7 computing order=4, npoints=6, at term=2 f^(4)(x[2])=(1/1h^4)( 1 * f(x[0]) + -4 * f(x[1]) + 6 * f(x[2]) + -4 * f(x[3]) + 1 * f(x[4]) + 0 * f(x[5] ) + O(h^2)f^(6)(z) err=0, h=0.0625, deriv=24, x=1.125, pwr=4 err=0, h=0.0625, deriv=135, x=1.125, pwr=5 err=-0.46875, h=0.0625, deriv=455.625, x=1.125, pwr=6 err=-0.117188, h=0.03125, deriv=406.406, x=1.0625, pwr=6 err=-3.69141, h=0.0625, deriv=1196.02, x=1.125, pwr=7 err=-0.871582, h=0.03125, deriv=1007.55, x=1.0625, pwr=7 err=-0.211487, h=0.015625, deriv=921.237, x=1.03125, pwr=7 computing order=4, npoints=6, at term=3 f^(4)(x[3])=(1/1h^4)( 0 * f(x[0]) + 1 * f(x[1]) + -4 * f(x[2]) + 6 * f(x[3]) + -4 * f(x[4]) + 1 * f(x[5] ) + O(h^2)f^(6)(z) err=0, h=0.0625, deriv=24, x=1.1875, pwr=4 err=0, h=0.0625, deriv=142.5, x=1.1875, pwr=5 err=-0.46875, h=0.0625, deriv=507.656, x=1.1875, pwr=6 err=-0.117188, h=0.03125, deriv=430.664, x=1.09375, pwr=6 err=-3.89648, h=0.0625, deriv=1406.63, x=1.1875, pwr=7 err=-0.897217, h=0.03125, deriv=1099.09, x=1.09375, pwr=7 err=-0.214691, h=0.015625, deriv=963.749, x=1.04688, pwr=7 computing order=4, npoints=6, at term=4 f^(4)(x[4])=(1/1h^4)( -1 * f(x[0]) + 6 * f(x[1]) + -14 * f(x[2]) + 16 * f(x[3]) + -9 * f(x[4]) + 2 * f(x[5] ) + O(h^2)f^(6)(z) err=0, h=0.0625, deriv=24, x=1.25, pwr=4 err=0, h=0.0625, deriv=150, x=1.25, pwr=5 err=2.34375, h=0.0625, deriv=562.5, x=1.25, pwr=6 err=0.585938, h=0.03125, deriv=455.625, x=1.125, pwr=6 err=19.2773, h=0.0625, deriv=1640.62, x=1.25, pwr=7 err=4.46045, h=0.03125, deriv=1196.02, x=1.125, pwr=7 err=1.07025, h=0.015625, deriv=1007.55, x=1.0625, pwr=7 computing order=4, npoints=6, at term=5 f^(4)(x[5])=(1/1h^4)( -2 * f(x[0]) + 11 * f(x[1]) + -24 * f(x[2]) + 26 * f(x[3]) + -14 * f(x[4]) + 3 * f(x[5] ) + O(h^2)f^(6)(z) err=0, h=0.0625, deriv=24, x=1.3125, pwr=4 err=0, h=0.0625, deriv=157.5, x=1.3125, pwr=5 err=7.96875, h=0.0625, deriv=620.156, x=1.3125, pwr=6 err=1.99219, h=0.03125, deriv=481.289, x=1.15625, pwr=6 err=67.0605, h=0.0625, deriv=1899.23, x=1.3125, pwr=7 err=15.3552, h=0.03125, deriv=1298.48, x=1.15625, pwr=7 err=3.66257, h=0.015625, deriv=1052.66, x=1.07812, pwr=7 computing order=5, npoints=6, at term=0 f^(5)(x[0])=(1/1h^5)( -1 * f(x[0]) + 5 * f(x[1]) + -10 * f(x[2]) + 10 * f(x[3]) + -5 * f(x[4]) + 1 * f(x[5] ) + O(h^1)f^(6)(z) err=0, h=0.03125, deriv=0, x=1, pwr=4 err=0, h=0.03125, deriv=120, x=1, pwr=5 err=-56.25, h=0.03125, deriv=720, x=1, pwr=6 err=-28.125, h=0.015625, deriv=720, x=1, pwr=6 err=-410.156, h=0.03125, deriv=2520, x=1, pwr=7 err=-200.977, h=0.015625, deriv=2520, x=1, pwr=7 err=-99.4629, h=0.0078125, deriv=2520, x=1, pwr=7 computing order=5, npoints=6, at term=1 f^(5)(x[1])=(1/1h^5)( -1 * f(x[0]) + 5 * f(x[1]) + -10 * f(x[2]) + 10 * f(x[3]) + -5 * f(x[4]) + 1 * f(x[5] ) + O(h^1)f^(6)(z) err=0, h=0.03125, deriv=0, x=1.03125, pwr=4 err=0, h=0.03125, deriv=120, x=1.03125, pwr=5 err=-33.75, h=0.03125, deriv=742.5, x=1.03125, pwr=6 err=-16.875, h=0.015625, deriv=731.25, x=1.01562, pwr=6 err=-250.195, h=0.03125, deriv=2679.96, x=1.03125, pwr=7 err=-121.611, h=0.015625, deriv=2599.37, x=1.01562, pwr=7 err=-59.9341, h=0.0078125, deriv=2559.53, x=1.00781, pwr=7 computing order=5, npoints=6, at term=2 f^(5)(x[2])=(1/1h^5)( -1 * f(x[0]) + 5 * f(x[1]) + -10 * f(x[2]) + 10 * f(x[3]) + -5 * f(x[4]) + 1 * f(x[5] ) + O(h^1)f^(6)(z) err=0, h=0.03125, deriv=0, x=1.0625, pwr=4 err=0, h=0.03125, deriv=120, x=1.0625, pwr=5 err=-11.25, h=0.03125, deriv=765, x=1.0625, pwr=6 err=-5.625, h=0.015625, deriv=742.5, x=1.03125, pwr=6 err=-85.3125, h=0.03125, deriv=2844.84, x=1.0625, pwr=7 err=-41.0156, h=0.015625, deriv=2679.96, x=1.03125, pwr=7 err=-20.0977, h=0.0078125, deriv=2599.37, x=1.01562, pwr=7 computing order=5, npoints=6, at term=3 f^(5)(x[3])=(1/1h^5)( -1 * f(x[0]) + 5 * f(x[1]) + -10 * f(x[2]) + 10 * f(x[3]) + -5 * f(x[4]) + 1 * f(x[5] ) + O(h^1)f^(6)(z) err=0, h=0.03125, deriv=0, x=1.09375, pwr=4 err=0, h=0.03125, deriv=120, x=1.09375, pwr=5 err=11.25, h=0.03125, deriv=787.5, x=1.09375, pwr=6 err=5.625, h=0.015625, deriv=753.75, x=1.04688, pwr=6 err=84.4922, h=0.03125, deriv=3014.65, x=1.09375, pwr=7 err=40.8105, h=0.015625, deriv=2761.79, x=1.04688, pwr=7 err=20.0464, h=0.0078125, deriv=2639.51, x=1.02344, pwr=7 computing order=5, npoints=6, at term=4 f^(5)(x[4])=(1/1h^5)( -1 * f(x[0]) + 5 * f(x[1]) + -10 * f(x[2]) + 10 * f(x[3]) + -5 * f(x[4]) + 1 * f(x[5] ) + O(h^1)f^(6)(z) err=0, h=0.03125, deriv=0, x=1.125, pwr=4 err=0, h=0.03125, deriv=120, x=1.125, pwr=5 err=33.75, h=0.03125, deriv=810, x=1.125, pwr=6 err=16.875, h=0.015625, deriv=765, x=1.0625, pwr=6 err=259.219, h=0.03125, deriv=3189.38, x=1.125, pwr=7 err=123.867, h=0.015625, deriv=2844.84, x=1.0625, pwr=7 err=60.498, h=0.0078125, deriv=2679.96, x=1.03125, pwr=7 computing order=5, npoints=6, at term=5 f^(5)(x[5])=(1/1h^5)( -1 * f(x[0]) + 5 * f(x[1]) + -10 * f(x[2]) + 10 * f(x[3]) + -5 * f(x[4]) + 1 * f(x[5] ) + O(h^1)f^(6)(z) err=0, h=0.03125, deriv=0, x=1.15625, pwr=4 err=0, h=0.03125, deriv=120, x=1.15625, pwr=5 err=56.25, h=0.03125, deriv=832.5, x=1.15625, pwr=6 err=28.125, h=0.015625, deriv=776.25, x=1.07812, pwr=6 err=438.867, h=0.03125, deriv=3369.02, x=1.15625, pwr=7 err=208.154, h=0.015625, deriv=2929.13, x=1.07812, pwr=7 err=101.257, h=0.0078125, deriv=2720.72, x=1.03906, pwr=7 computing order=5, npoints=7, at term=0 f^(5)(x[0])=(1/2h^5)( -7 * f(x[0]) + 40 * f(x[1]) + -95 * f(x[2]) + 120 * f(x[3]) + -85 * f(x[4]) + 32 * f(x[5]) + -5 * f(x[6] ) + O(h^2)f^(7)(z) err=0, h=0.03125, deriv=120, x=1, pwr=5 err=0, h=0.03125, deriv=720, x=1, pwr=6 err=20.5078, h=0.03125, deriv=2520, x=1, pwr=7 err=5.12695, h=0.015625, deriv=2520, x=1, pwr=7 err=174.829, h=0.03125, deriv=6720, x=1, pwr=8 err=42.3615, h=0.015625, deriv=6720, x=1, pwr=8 err=10.4224, h=0.0078125, deriv=6720, x=1, pwr=8 computing order=5, npoints=7, at term=1 f^(5)(x[1])=(1/2h^5)( -5 * f(x[0]) + 28 * f(x[1]) + -65 * f(x[2]) + 80 * f(x[3]) + -55 * f(x[4]) + 20 * f(x[5]) + -3 * f(x[6] ) + O(h^2)f^(7)(z) err=0, h=0.03125, deriv=120, x=1.03125, pwr=5 err=0, h=0.03125, deriv=742.5, x=1.03125, pwr=6 err=8.20312, h=0.03125, deriv=2679.96, x=1.03125, pwr=7 err=2.05078, h=0.015625, deriv=2599.37, x=1.01562, pwr=7 err=70.752, h=0.03125, deriv=7369.89, x=1.03125, pwr=8 err=17.0471, h=0.015625, deriv=7039.95, x=1.01562, pwr=8 err=4.18173, h=0.0078125, deriv=6878.73, x=1.00781, pwr=8 computing order=5, npoints=7, at term=2 f^(5)(x[2])=(1/2h^5)( -3 * f(x[0]) + 16 * f(x[1]) + -35 * f(x[2]) + 40 * f(x[3]) + -25 * f(x[4]) + 8 * f(x[5]) + -1 * f(x[6] ) + O(h^2)f^(7)(z) err=0, h=0.03125, deriv=120, x=1.0625, pwr=5 err=0, h=0.03125, deriv=765, x=1.0625, pwr=6 err=0.820312, h=0.03125, deriv=2844.84, x=1.0625, pwr=7 err=0.205078, h=0.015625, deriv=2679.96, x=1.03125, pwr=7 err=7.28027, h=0.03125, deriv=8060.39, x=1.0625, pwr=8 err=1.73035, h=0.015625, deriv=7369.89, x=1.03125, pwr=8 err=0.421509, h=0.0078125, deriv=7039.95, x=1.01562, pwr=8 computing order=5, npoints=7, at term=3 f^(5)(x[3])=(1/2h^5)( -1 * f(x[0]) + 4 * f(x[1]) + -5 * f(x[2]) + 0 * f(x[3]) + 5 * f(x[4]) + -4 * f(x[5]) + 1 * f(x[6] ) + O(h^2)f^(7)(z) err=0, h=0.03125, deriv=120, x=1.09375, pwr=5 err=0, h=0.03125, deriv=787.5, x=1.09375, pwr=6 err=-1.64062, h=0.03125, deriv=3014.65, x=1.09375, pwr=7 err=-0.410156, h=0.015625, deriv=2761.79, x=1.04688, pwr=7 err=-14.3555, h=0.03125, deriv=8792.72, x=1.09375, pwr=8 err=-3.43506, h=0.015625, deriv=7709.99, x=1.04688, pwr=8 err=-0.839539, h=0.0078125, deriv=7203.66, x=1.02344, pwr=8 computing order=5, npoints=7, at term=4 f^(5)(x[4])=(1/2h^5)( 1 * f(x[0]) + -8 * f(x[1]) + 25 * f(x[2]) + -40 * f(x[3]) + 35 * f(x[4]) + -16 * f(x[5]) + 3 * f(x[6] ) + O(h^2)f^(7)(z) err=0, h=0.03125, deriv=120, x=1.125, pwr=5 err=0, h=0.03125, deriv=810, x=1.125, pwr=6 err=0.820312, h=0.03125, deriv=3189.38, x=1.125, pwr=7 err=0.205078, h=0.015625, deriv=2844.84, x=1.0625, pwr=7 err=7.0752, h=0.03125, deriv=9568.12, x=1.125, pwr=8 err=1.70471, h=0.015625, deriv=8060.39, x=1.0625, pwr=8 err=0.418091, h=0.0078125, deriv=7369.89, x=1.03125, pwr=8 computing order=5, npoints=7, at term=5 f^(5)(x[5])=(1/2h^5)( 3 * f(x[0]) + -20 * f(x[1]) + 55 * f(x[2]) + -80 * f(x[3]) + 65 * f(x[4]) + -28 * f(x[5]) + 5 * f(x[6] ) + O(h^2)f^(7)(z) err=0, h=0.03125, deriv=120, x=1.15625, pwr=5 err=0, h=0.03125, deriv=832.5, x=1.15625, pwr=6 err=8.20312, h=0.03125, deriv=3369.02, x=1.15625, pwr=7 err=2.05078, h=0.015625, deriv=2929.13, x=1.07812, pwr=7 err=72.8027, h=0.03125, deriv=10387.8, x=1.15625, pwr=8 err=17.3035, h=0.015625, deriv=8421.25, x=1.07812, pwr=8 err=4.21347, h=0.0078125, deriv=7538.66, x=1.03906, pwr=8 computing order=5, npoints=7, at term=6 f^(5)(x[6])=(1/2h^5)( 5 * f(x[0]) + -32 * f(x[1]) + 85 * f(x[2]) + -120 * f(x[3]) + 95 * f(x[4]) + -40 * f(x[5]) + 7 * f(x[6] ) + O(h^2)f^(7)(z) err=0, h=0.03125, deriv=120, x=1.1875, pwr=5 err=0, h=0.03125, deriv=855, x=1.1875, pwr=6 err=20.5078, h=0.03125, deriv=3553.59, x=1.1875, pwr=7 err=5.12695, h=0.015625, deriv=3014.65, x=1.09375, pwr=7 err=184.058, h=0.03125, deriv=11253, x=1.1875, pwr=8 err=43.515, h=0.015625, deriv=8792.72, x=1.09375, pwr=8 err=10.5662, h=0.0078125, deriv=7709.99, x=1.04688, pwr=8 computing order=6, npoints=7, at term=0 f^(6)(x[0])=(1/1h^6)( 1 * f(x[0]) + -6 * f(x[1]) + 15 * f(x[2]) + -20 * f(x[3]) + 15 * f(x[4]) + -6 * f(x[5]) + 1 * f(x[6] ) + O(h^1)f^(7)(z) err=0, h=0.015625, deriv=0, x=1, pwr=5 err=0, h=0.015625, deriv=720, x=1, pwr=6 err=-236.25, h=0.015625, deriv=5040, x=1, pwr=7 err=-118.125, h=0.0078125, deriv=5040, x=1, pwr=7 err=-1936.76, h=0.015625, deriv=20160, x=1, pwr=8 err=-956.703, h=0.0078125, deriv=20160, x=1, pwr=8 err=-475.375, h=0.00390625, deriv=20160, x=1, pwr=8 computing order=6, npoints=7, at term=1 f^(6)(x[1])=(1/1h^6)( 1 * f(x[0]) + -6 * f(x[1]) + 15 * f(x[2]) + -20 * f(x[3]) + 15 * f(x[4]) + -6 * f(x[5]) + 1 * f(x[6] ) + O(h^1)f^(7)(z) err=0, h=0.015625, deriv=0, x=1.01562, pwr=5 err=0, h=0.015625, deriv=720, x=1.01562, pwr=6 err=-157.5, h=0.015625, deriv=5118.75, x=1.01562, pwr=7 err=-78.75, h=0.0078125, deriv=5079.38, x=1.00781, pwr=7 err=-1301.84, h=0.015625, deriv=20794.9, x=1.01562, pwr=8 err=-640.473, h=0.0078125, deriv=20476.2, x=1.00781, pwr=8 err=-317.567, h=0.00390625, deriv=20317.8, x=1.00391, pwr=8 computing order=6, npoints=7, at term=2 f^(6)(x[2])=(1/1h^6)( 1 * f(x[0]) + -6 * f(x[1]) + 15 * f(x[2]) + -20 * f(x[3]) + 15 * f(x[4]) + -6 * f(x[5]) + 1 * f(x[6] ) + O(h^1)f^(7)(z) err=0, h=0.015625, deriv=0, x=1.03125, pwr=5 err=0, h=0.015625, deriv=720, x=1.03125, pwr=6 err=-78.75, h=0.015625, deriv=5197.5, x=1.03125, pwr=7 err=-39.375, h=0.0078125, deriv=5118.75, x=1.01562, pwr=7 err=-657.07, h=0.015625, deriv=21439.7, x=1.03125, pwr=8 err=-321.781, h=0.0078125, deriv=20794.9, x=1.01562, pwr=8 err=-159.145, h=0.00390625, deriv=20476.2, x=1.00781, pwr=8 computing order=6, npoints=7, at term=3 f^(6)(x[3])=(1/1h^6)( 1 * f(x[0]) + -6 * f(x[1]) + 15 * f(x[2]) + -20 * f(x[3]) + 15 * f(x[4]) + -6 * f(x[5]) + 1 * f(x[6] ) + O(h^1)f^(7)(z) err=0, h=0.015625, deriv=0, x=1.04688, pwr=5 err=0, h=0.015625, deriv=720, x=1.04688, pwr=6 err=0, h=0.015625, deriv=5276.25, x=1.04688, pwr=7 err=0, h=0.0078125, deriv=5158.12, x=1.02344, pwr=7 err=-2.46094, h=0.015625, deriv=22094.3, x=1.04688, pwr=8 err=-0.628906, h=0.0078125, deriv=21116.1, x=1.02344, pwr=8 err=-0.106445, h=0.00390625, deriv=20635.3, x=1.01172, pwr=8 computing order=6, npoints=7, at term=4 f^(6)(x[4])=(1/1h^6)( 1 * f(x[0]) + -6 * f(x[1]) + 15 * f(x[2]) + -20 * f(x[3]) + 15 * f(x[4]) + -6 * f(x[5]) + 1 * f(x[6] ) + O(h^1)f^(7)(z) err=0, h=0.015625, deriv=0, x=1.0625, pwr=5 err=0, h=0.015625, deriv=720, x=1.0625, pwr=6 err=78.75, h=0.015625, deriv=5355, x=1.0625, pwr=7 err=39.375, h=0.0078125, deriv=5197.5, x=1.03125, pwr=7 err=661.992, h=0.015625, deriv=22758.8, x=1.0625, pwr=8 err=322.984, h=0.0078125, deriv=21439.7, x=1.03125, pwr=8 err=159.547, h=0.00390625, deriv=20794.9, x=1.01562, pwr=8 computing order=6, npoints=7, at term=5 f^(6)(x[5])=(1/1h^6)( 1 * f(x[0]) + -6 * f(x[1]) + 15 * f(x[2]) + -20 * f(x[3]) + 15 * f(x[4]) + -6 * f(x[5]) + 1 * f(x[6] ) + O(h^1)f^(7)(z) err=0, h=0.015625, deriv=0, x=1.07812, pwr=5 err=0, h=0.015625, deriv=720, x=1.07812, pwr=6 err=157.5, h=0.015625, deriv=5433.75, x=1.07812, pwr=7 err=78.75, h=0.0078125, deriv=5236.88, x=1.03906, pwr=7 err=1336.29, h=0.015625, deriv=23433, x=1.07812, pwr=8 err=649.059, h=0.0078125, deriv=21765.8, x=1.03906, pwr=8 err=319.815, h=0.00390625, deriv=20955.2, x=1.01953, pwr=8 computing order=6, npoints=7, at term=6 f^(6)(x[6])=(1/1h^6)( 1 * f(x[0]) + -6 * f(x[1]) + 15 * f(x[2]) + -20 * f(x[3]) + 15 * f(x[4]) + -6 * f(x[5]) + 1 * f(x[6] ) + O(h^1)f^(7)(z) err=0, h=0.015625, deriv=0, x=1.09375, pwr=5 err=0, h=0.015625, deriv=720, x=1.09375, pwr=6 err=236.25, h=0.015625, deriv=5512.5, x=1.09375, pwr=7 err=118.125, h=0.0078125, deriv=5276.25, x=1.04688, pwr=7 err=2020.43, h=0.015625, deriv=24117.2, x=1.09375, pwr=8 err=977.594, h=0.0078125, deriv=22094.3, x=1.04688, pwr=8 err=480.699, h=0.00390625, deriv=21116.1, x=1.02344, pwr=8 computing order=6, npoints=8, at term=0 f^(6)(x[0])=(1/1h^6)( 4 * f(x[0]) + -27 * f(x[1]) + 78 * f(x[2]) + -125 * f(x[3]) + 120 * f(x[4]) + -69 * f(x[5]) + 22 * f(x[6]) + -3 * f(x[7] ) + O(h^2)f^(8)(z) err=0, h=0.015625, deriv=720, x=1, pwr=6 err=0, h=0.015625, deriv=5040, x=1, pwr=7 err=56.6023, h=0.015625, deriv=20160, x=1, pwr=8 err=14.0977, h=0.0078125, deriv=20160, x=1, pwr=8 err=528.795, h=0.015625, deriv=60480, x=1, pwr=9 err=129.73, h=0.0078125, deriv=60480, x=1, pwr=9 err=38.25, h=0.00390625, deriv=60480, x=1, pwr=9 computing order=6, npoints=8, at term=1 f^(6)(x[1])=(1/1h^6)( 3 * f(x[0]) + -20 * f(x[1]) + 57 * f(x[2]) + -90 * f(x[3]) + 85 * f(x[4]) + -48 * f(x[5]) + 15 * f(x[6]) + -2 * f(x[7] ) + O(h^2)f^(8)(z) err=0, h=0.015625, deriv=720, x=1.01562, pwr=6 err=0, h=0.015625, deriv=5118.75, x=1.01562, pwr=7 err=27.0698, h=0.015625, deriv=20794.9, x=1.01562, pwr=8 err=6.76953, h=0.0078125, deriv=20476.2, x=1.00781, pwr=8 err=254.361, h=0.015625, deriv=63359.5, x=1.01562, pwr=9 err=62.2359, h=0.0078125, deriv=61908.6, x=1.00781, pwr=9 err=20.0222, h=0.00390625, deriv=61191.5, x=1.00391, pwr=9 computing order=6, npoints=8, at term=2 f^(6)(x[2])=(1/1h^6)( 2 * f(x[0]) + -13 * f(x[1]) + 36 * f(x[2]) + -55 * f(x[3]) + 50 * f(x[4]) + -27 * f(x[5]) + 8 * f(x[6]) + -1 * f(x[7] ) + O(h^2)f^(8)(z) err=0, h=0.015625, deriv=720, x=1.03125, pwr=6 err=0, h=0.015625, deriv=5197.5, x=1.03125, pwr=7 err=7.38306, h=0.015625, deriv=21439.7, x=1.03125, pwr=8 err=1.80859, h=0.0078125, deriv=20794.9, x=1.01562, pwr=8 err=69.9063, h=0.015625, deriv=66329, x=1.03125, pwr=9 err=17.0237, h=0.0078125, deriv=63359.5, x=1.01562, pwr=9 err=7.35306, h=0.00390625, deriv=61908.6, x=1.00781, pwr=9 computing order=6, npoints=8, at term=3 f^(6)(x[3])=(1/1h^6)( 1 * f(x[0]) + -6 * f(x[1]) + 15 * f(x[2]) + -20 * f(x[3]) + 15 * f(x[4]) + -6 * f(x[5]) + 1 * f(x[6]) + 0 * f(x[7] ) + O(h^2)f^(8)(z) err=0, h=0.015625, deriv=720, x=1.04688, pwr=6 err=0, h=0.015625, deriv=5276.25, x=1.04688, pwr=7 err=-2.46094, h=0.015625, deriv=22094.3, x=1.04688, pwr=8 err=-0.628906, h=0.0078125, deriv=21116.1, x=1.02344, pwr=8 err=-23.1866, h=0.015625, deriv=69389.9, x=1.04688, pwr=9 err=-5.67154, h=0.0078125, deriv=64832.9, x=1.02344, pwr=9 err=-0.235676, h=0.00390625, deriv=62631.3, x=1.01172, pwr=9 computing order=6, npoints=8, at term=4 f^(6)(x[4])=(1/1h^6)( 0 * f(x[0]) + 1 * f(x[1]) + -6 * f(x[2]) + 15 * f(x[3]) + -20 * f(x[4]) + 15 * f(x[5]) + -6 * f(x[6]) + 1 * f(x[7] ) + O(h^2)f^(8)(z) err=0, h=0.015625, deriv=720, x=1.0625, pwr=6 err=0, h=0.015625, deriv=5355, x=1.0625, pwr=7 err=-2.46094, h=0.015625, deriv=22758.8, x=1.0625, pwr=8 err=-0.605469, h=0.0078125, deriv=21439.7, x=1.03125, pwr=8 err=-23.5325, h=0.015625, deriv=72543.5, x=1.0625, pwr=9 err=-5.71289, h=0.0078125, deriv=66329, x=1.03125, pwr=9 err=-2.72241, h=0.00390625, deriv=63359.5, x=1.01562, pwr=9 computing order=6, npoints=8, at term=5 f^(6)(x[5])=(1/1h^6)( -1 * f(x[0]) + 8 * f(x[1]) + -27 * f(x[2]) + 50 * f(x[3]) + -55 * f(x[4]) + 36 * f(x[5]) + -13 * f(x[6]) + 2 * f(x[7] ) + O(h^2)f^(8)(z) err=0, h=0.015625, deriv=720, x=1.07812, pwr=6 err=0, h=0.015625, deriv=5433.75, x=1.07812, pwr=7 err=7.3833, h=0.015625, deriv=23433, x=1.07812, pwr=8 err=1.81641, h=0.0078125, deriv=21765.8, x=1.03906, pwr=8 err=70.2521, h=0.015625, deriv=75791.3, x=1.07812, pwr=9 err=17.0932, h=0.0078125, deriv=67848, x=1.03906, pwr=9 err=3.91448, h=0.00390625, deriv=64093.4, x=1.01953, pwr=9 computing order=6, npoints=8, at term=6 f^(6)(x[6])=(1/1h^6)( -2 * f(x[0]) + 15 * f(x[1]) + -48 * f(x[2]) + 85 * f(x[3]) + -90 * f(x[4]) + 57 * f(x[5]) + -20 * f(x[6]) + 3 * f(x[7] ) + O(h^2)f^(8)(z) err=0, h=0.015625, deriv=720, x=1.09375, pwr=6 err=0, h=0.015625, deriv=5512.5, x=1.09375, pwr=7 err=27.0706, h=0.015625, deriv=24117.2, x=1.09375, pwr=8 err=6.73047, h=0.0078125, deriv=22094.3, x=1.04688, pwr=8 err=259.552, h=0.015625, deriv=79134.5, x=1.09375, pwr=9 err=62.8738, h=0.0078125, deriv=69389.9, x=1.04688, pwr=9 err=9.69662, h=0.00390625, deriv=64832.9, x=1.02344, pwr=9 computing order=6, npoints=8, at term=7 f^(6)(x[7])=(1/1h^6)( -3 * f(x[0]) + 22 * f(x[1]) + -69 * f(x[2]) + 120 * f(x[3]) + -125 * f(x[4]) + 78 * f(x[5]) + -27 * f(x[6]) + 4 * f(x[7] ) + O(h^2)f^(8)(z) err=0, h=0.015625, deriv=720, x=1.10938, pwr=6 err=0, h=0.015625, deriv=5591.25, x=1.10938, pwr=7 err=56.6006, h=0.015625, deriv=24811.2, x=1.10938, pwr=8 err=14.1836, h=0.0078125, deriv=22425.3, x=1.05469, pwr=8 err=545.752, h=0.015625, deriv=82574.7, x=1.10938, pwr=9 err=131.888, h=0.0078125, deriv=70955, x=1.05469, pwr=9 err=25.1457, h=0.00390625, deriv=65578.1, x=1.02734, pwr=9 many formulas checked against Abramowitz and Stegun, Handbook of Mathematical Functions Table 25.2