07-Mar-2024 11:44:55 conte_deboor_test(): MATLAB/Octave version 5.2.0. Test conte_deboor() bisect_test Test bisect(): X = 2.000000e+00, F(X) = 5.000000e+00 X = 3.000000e+00, F(X) = 2.300000e+01 F is of the same sign at the two endpoints 2.000000e+00, 3.000000e+00 X = 1.000000e+00, F(X) = -1.000000e+00 X = 2.000000e+00, F(X) = 5.000000e+00 X = 1.500000e+00, F(X) = 8.750000e-01 X = 1.250000e+00, F(X) = -2.968750e-01 X = 1.375000e+00, F(X) = 2.246094e-01 X = 1.312500e+00, F(X) = -5.151367e-02 X = 1.343750e+00, F(X) = 8.261108e-02 X = 1.328125e+00, F(X) = 1.457596e-02 X = 1.320312e+00, F(X) = -1.871061e-02 X = 1.324219e+00, F(X) = -2.127945e-03 X = 1.326172e+00, F(X) = 6.208830e-03 X = 1.325195e+00, F(X) = 2.036651e-03 X = 1.324707e+00, F(X) = -4.659488e-05 X = 1.324951e+00, F(X) = 9.947910e-04 X = 1.324829e+00, F(X) = 4.740388e-04 X = 1.324768e+00, F(X) = 2.137072e-04 X = 1.324738e+00, F(X) = 8.355244e-05 X = 1.324722e+00, F(X) = 1.847785e-05 X = 1.324715e+00, F(X) = -1.405875e-05 X = 1.324718e+00, F(X) = 2.209495e-06 X = 1.324717e+00, F(X) = -5.924640e-06 X = 1.324718e+00, F(X) = -1.857576e-06 X = 1.324718e+00, F(X) = 1.759583e-07 X = 1.324718e+00, F(X) = -8.408093e-07 X = 1.324718e+00, F(X) = -3.324255e-07 X = 1.324718e+00, F(X) = -7.823363e-08 X = 1.324718e+00, F(X) = 4.886233e-08 X = 1.324718e+00, F(X) = -1.468565e-08 X = 1.324718e+00, F(X) = 1.708834e-08 X = 1.324718e+00, F(X) = 1.201345e-09 X = 1.324718e+00, F(X) = -6.742154e-09 X = 1.324718e+00, F(X) = -2.770405e-09 X = 1.324718e+00, F(X) = -7.845300e-10 X = 1.324718e+00, F(X) = 2.084073e-10 X = 1.324718e+00, F(X) = -2.880616e-10 X = 1.324718e+00, F(X) = -3.982714e-11 X = 1.324718e+00, F(X) = 8.429013e-11 X = 1.324718e+00, F(X) = 2.223133e-11 X = 1.324718e+00, F(X) = -8.797740e-12 X = 1.324718e+00, F(X) = 6.716849e-12 X = 1.324718e+00, F(X) = -1.040501e-12 X = 1.324718e+00, F(X) = 2.838174e-12 X = 1.324718e+00, F(X) = 8.988366e-13 X = 1.324718e+00, F(X) = -7.083223e-14 X = 1.324718e+00, F(X) = 4.138911e-13 X = 1.324718e+00, F(X) = 1.714184e-13 X = 1.324718e+00, F(X) = 5.040413e-14 X = 1.324718e+00, F(X) = -1.043610e-14 X = 1.324718e+00, F(X) = 1.998401e-14 X = 1.324718e+00, F(X) = 4.884981e-15 X = 1.324718e+00, F(X) = -2.553513e-15 X = 1.324718e+00, F(X) = 1.110223e-15 X = 1.324718e+00, F(X) = -7.771561e-16 X = 1.324718e+00, F(X) = 2.220446e-16 IFLAG = 1 (should be 1) The zero is 1.324718e+00 plus/minus 1.110223e-16 X = 1.000000e+00, F(X) = -1.000000e+00 X = 2.000000e+00, F(X) = 5.000000e+00 X = 1.500000e+00, F(X) = 8.750000e-01 X = 1.250000e+00, F(X) = -2.968750e-01 X = 1.375000e+00, F(X) = 2.246094e-01 X = 1.312500e+00, F(X) = -5.151367e-02 X = 1.343750e+00, F(X) = 8.261108e-02 X = 1.328125e+00, F(X) = 1.457596e-02 X = 1.320312e+00, F(X) = -1.871061e-02 X = 1.324219e+00, F(X) = -2.127945e-03 X = 1.326172e+00, F(X) = 6.208830e-03 X = 1.325195e+00, F(X) = 2.036651e-03 X = 1.324707e+00, F(X) = -4.659488e-05 X = 1.324951e+00, F(X) = 9.947910e-04 X = 1.324829e+00, F(X) = 4.740388e-04 X = 1.324768e+00, F(X) = 2.137072e-04 X = 1.324738e+00, F(X) = 8.355244e-05 X = 1.324722e+00, F(X) = 1.847785e-05 X = 1.324715e+00, F(X) = -1.405875e-05 X = 1.324718e+00, F(X) = 2.209495e-06 X = 1.324717e+00, F(X) = -5.924640e-06 The zero is 1.324718e+00 plus/minus 9.536743e-07 calccf_test Test calccf() 0 2.51 0.5 2.91 1 3.3 1.5 3.68 2 4.04 2.5 4.38 3 4.7 3.5 4.98 4 5.22 4.5 5.4 5 5.54 5.5 5.67 6 5.78 6.5 5.44 7 5.4 7.5 5.47 8 5.57 8.5 5.65 9 5.7 9.5 5.74 10 5.78 10.5 5.81 11 5.83 11.5 5.84 12 5.84 12.5 5.83 13 5.81 13.5 5.79 14 5.75 14.5 5.71 15 5.65 15.5 5.57 16 5.48 16.5 5.36 17 5.22 17.5 5.07 18 4.9 18.5 4.72 19 4.53 19.5 4.34 20 4.16 Graphics saved as "calccf_test.png" cheb_test: Test cheb() X cheb(exp)(x) exp(x) 0.000000 1 1 0.500000 1.64872 1.64872 1.000000 2.71828 2.71828 1.500000 4.48169 4.48169 2.000000 7.38906 7.38906 2.500000 12.1825 12.1825 3.000000 20.0855 20.0855 3.500000 33.1152 33.1155 4.000000 54.5966 54.5982 4.500000 90.0076 90.0171 5.000000 148.366 148.413 cheb_test: Normal end of execution. factor_cd_test: Test factor_cd() A = Diagonal Matrix 1 0 0 0 1 0 0 0 1 W = 1 0 0 0 1 0 0 0 1 IPIVOT = 1 2 3 IFLAG = 1 ans = 0 0 0 0 0 0 0 0 0 A = 1 1 1 1 1 1 1 1 1 W = 1 1 1 1 0 0 1 0 0 IPIVOT = 1 2 3 IFLAG = 0 ans = 0 0 0 0 0 0 0 0 0 A = 1 1 1 1 1 1 1 1 1 W = 1 1 1 1 0 0 1 0 0 IPIVOT = 1 2 3 IFLAG = 0 ans = 0 0 0 0 0 0 0 0 0 A = 0 1 0 0 1 1 1 1 1 W = 1 1 1 0 1 1 0 1 -1 IPIVOT = 3 2 1 IFLAG = -1 ans = 0 0 0 0 0 0 0 0 0 A = 0.632400 0.741535 0.908128 0.057138 0.715639 0.629997 0.336090 0.717292 0.471943 W = 0.632400 0.741535 0.908128 0.090352 0.648641 0.547946 0.531451 0.498277 -0.283712 IPIVOT = 1 2 3 IFLAG = 1 ans = 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 5.5511e-17 0.0000e+00 0.0000e+00 W = 1 1 1 0 1 1 0 0 1 IPIVOT = 1 2 3 IFLAG = 1 ans = 1.1102e-16 -4.1633e-17 0.0000e+00 W = 1.00000 1.00000 0.00000 0.50000 -0.50000 0.00000 1.00000 -0.00000 1.00000 IPIVOT = 2 1 3 IFLAG = -1 ans = 0.0000e+00 -4.1633e-17 1.1102e-16 W = 1 0 0 1 1 0 0 1 1 IPIVOT = 1 2 3 IFLAG = 1 ans = 0.0000e+00 -4.1633e-17 0.0000e+00 W = 0.920507 0.452659 0.424381 0.039542 0.779837 0.490970 0.927429 -0.041781 0.514046 IPIVOT = 3 2 1 IFLAG = -1 ans = 2.2204e-16 1.1102e-16 -2.2204e-16 W = 1 0 0 0 1 0 0 0 1 IPIVOT = 3 1 2 IFLAG = 1 ans = 0 0 0 fft_cd_test: Test fft_cd() horner_test Test horner() X P(X) Tan(X) 1.5000 14.100580 14.101420 0.7500 0.931176 0.931596 0.0000 -0.000163 0.000000 -0.7500 -0.931625 -0.931596 -1.5000 -14.101400 -14.101420 inverse_test: Demonstrate computation of inverse matrix using factor_cd() and subst() the computed inverse is AINV = 2.5000e-01 5.5511e-17 -2.5000e-01 5.0000e-01 -2.0000e-01 1.0000e-01 1.0000e+00 -6.0000e-01 -2.0000e-01 invitr_test: Test invitr() EIGENVALUE EIGENVECTOR COMPONENTS 10.828860 0.156204 0.029183 -0.152965 0.036365 0.158713 0.271722 -0.615101 0.332062 1.127093 0.661906 -0.892925 0.561759 0.852585 -0.439612 -1.150160 -1.130634 1.381377 0.587930 -0.865712 0.475976 1.015111 1.315620 -1.485739 0.593547 0.862426 -0.485905 -0.932922 -1.385090 1.515425 0.594726 -0.858708 0.489359 0.891927 1.412937 -1.523862 0.594972 0.856717 -0.490762 -0.872554 -1.424718 1.526538 0.595022 -0.855810 0.491409 0.863580 1.429889 -1.527472 0.595033 0.855424 -0.491725 -0.859455 -1.432201 1.527822 0.595035 -0.855267 0.491883 0.857566 1.433244 -1.527959 0.595036 0.855204 -0.491962 -0.856701 -1.433717 1.528014 0.595036 -0.855179 0.492002 0.856305 1.433932 -1.528036 invitr() returned iflag = -1 eigenvalue estimate e is 0.595036 eigenvector estimate v is: -0.855179 0.492002 0.856305 1.43393 -1.52804 Norm of residual A*v-v*e = 0.000183432 lgndre_test: Test lgndre() N Integral estimate 1 16 2 34.2785862786 3 58.6383250916 4 13.4796679462 5 22.3889300383 6 37.8299190562 7 32.7214135442 8 24.4540386939 9 26.9911822132 10 32.8239686913 lgndre_test: Normal end of execution. mrgfls_test Test mrgfls() X = 2.000000e+00, F(X) = 5.000000e+00 X = 3.000000e+00, F(X) = 2.300000e+01 F is of the same sign at the two endpoints 2.000000e+00, 3.000000e+00 The zero is NaN within [2.000000e+00 .. 3.000000e+00], IFLAG = -1 X = 1.000000e+00, F(X) = -1.000000e+00 X = 2.000000e+00, F(X) = 5.000000e+00 X = 1.166667e+00, F(X) = -5.787037e-01 X = 1.323308e+00, F(X) = -6.003901e-03 X = 1.326543e+00, F(X) = 7.796236e-03 The zero is 1.326543e+00 within [1.323308e+00 .. 1.326543e+00], IFLAG = 0 X = 1.000000e+00, F(X) = -1.000000e+00 X = 2.000000e+00, F(X) = 5.000000e+00 X = 1.166667e+00, F(X) = -5.787037e-01 X = 1.323308e+00, F(X) = -6.003901e-03 X = 1.326543e+00, F(X) = 7.796236e-03 X = 1.324716e+00, F(X) = -1.022135e-05 X = 1.324718e+00, F(X) = -1.736197e-08 The zero is 1.324718e+00 within [1.324718e+00 .. 1.326543e+00], IFLAG = 1 X = 1.000000e+00, F(X) = -1.000000e+00 X = 2.000000e+00, F(X) = 5.000000e+00 X = 1.166667e+00, F(X) = -5.787037e-01 X = 1.323308e+00, F(X) = -6.003901e-03 X = 1.326543e+00, F(X) = 7.796236e-03 X = 1.324716e+00, F(X) = -1.022135e-05 no convergence in 4 iterations The zero is 1.324716e+00 within [1.324716e+00 .. 1.326543e+00], IFLAG = 2 muller_test: Test muller() Example 3.13 ROOT F(ROOT) REAL IMAGINARY REAL IMAGINARY PART PART PART PART 2.4048256e+00 0.0000000e+00 1.4846886e-15 0.0000000e+00 5.5200781e+00 0.0000000e+00 -4.7081011e-14 0.0000000e+00 8.6537279e+00 0.0000000e+00 2.0233708e-12 0.0000000e+00 Example 3.14 ROOT F(ROOT) REAL IMAGINARY REAL IMAGINARY PART PART PART PART 1.2134117e+00 0.0000000e+00 -1.7763568e-15 0.0000000e+00 -6.0670583e-01 1.4506122e+00 -2.2204460e-15 4.4408921e-16 -6.0670583e-01 -1.4506122e+00 -4.4408921e-16 6.6613381e-15 Example 3.15 ROOT F(ROOT) REAL IMAGINARY REAL IMAGINARY PART PART PART PART 1.0000000e+00 -1.0000000e+00 1.7763568e-15 -8.8817842e-16 1.0000000e+00 1.0000000e+00 1.0569323e-13 -2.0383695e-13 -9.0941228e-12 1.4142136e+00 -8.5338758e-10 -1.2774036e-09 8.3273031e-11 -1.4142136e+00 -8.5326501e-10 -1.2778039e-09 1.7000000e+00 1.3032728e-10 6.3396577e-10 9.4957759e-10 Example 3.16 ROOT F(ROOT) REAL IMAGINARY REAL IMAGINARY PART PART PART PART 2.0000000e+00 -3.5073458e-16 0.0000000e+00 4.2088149e-14 3.0000000e+00 7.4072987e-11 -5.2023097e-10 3.5555034e-09 1.0000000e+00 -1.7020175e-16 0.0000000e+00 -1.2254526e-13 6.0000000e+00 -5.4934304e-15 3.8198777e-11 6.5921165e-13 5.0000000e+00 -7.6380671e-13 -1.5370460e-10 -3.6662722e-11 4.0000000e+00 -2.3674029e-10 -1.2669261e-09 8.5226503e-09 7.0000000e+00 -1.3038021e-20 3.7289283e-11 -9.3873750e-18 Example 3.17 ROOT F(ROOT) REAL IMAGINARY REAL IMAGINARY PART PART PART PART -1.4142136e+00 0.0000000e+00 -1.4551915e-11 0.0000000e+00 1.4142136e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 -2.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 2.0000000e+00 0.0000000e+00 -4.3655746e-11 0.0000000e+00 -8.0000000e+00 0.0000000e+00 1.4973921e-08 0.0000000e+00 8.0000000e+00 0.0000000e+00 1.9688741e-08 0.0000000e+00 -1.0000000e+01 0.0000000e+00 -1.4333636e-09 0.0000000e+00 1.0000000e+01 0.0000000e+00 0.0000000e+00 0.0000000e+00 ortpol_test Test ortpol() X P(X)-F(X) -1.0000, 6.667520e-03 -0.9000, 2.766170e-04 -0.8000, -3.292904e-03 -0.7000, -4.674874e-03 -0.6000, -4.458116e-03 -0.5000, -3.181712e-03 -0.4000, -1.329771e-03 -0.3000, 6.743520e-04 -0.2000, 2.474446e-03 -0.1000, 3.788506e-03 0.0000, 4.416536e-03 0.1000, 4.249174e-03 0.2000, 3.277226e-03 0.3000, 1.602199e-03 0.4000, -5.520572e-04 0.5000, -2.826479e-03 0.6000, -4.712575e-03 0.7000, -5.536703e-03 0.8000, -4.442708e-03 0.9000, -3.727226e-04 1.0000, 7.954046e-03 Orthogonal polynomial coefficients: 1.1937 1.1140 0.5402 0.1770 ortval_test Test ortval() X P(X)-F(X) -1.0000, 6.667520e-03 -0.9000, 2.766170e-04 -0.8000, -3.292904e-03 -0.7000, -4.674874e-03 -0.6000, -4.458116e-03 -0.5000, -3.181712e-03 -0.4000, -1.329771e-03 -0.3000, 6.743520e-04 -0.2000, 2.474446e-03 -0.1000, 3.788506e-03 0.0000, 4.416536e-03 0.1000, 4.249174e-03 0.2000, 3.277226e-03 0.3000, 1.602199e-03 0.4000, -5.520572e-04 0.5000, -2.826479e-03 0.6000, -4.712575e-03 0.7000, -5.536703e-03 0.8000, -4.442708e-03 0.9000, -3.727226e-04 1.0000, 7.954046e-03 Orthogonal polynomial coefficients: 1 1.1937 2 1.1140 3 0.5402 4 0.1770 Compare approximating polynomial to f(x): x p(x) f(x) -1.0000 0.3612 0.3679 -0.9333 0.3912 0.3932 -0.8667 0.4215 0.4204 -0.8000 0.4526 0.4493 -0.7333 0.4847 0.4803 -0.6667 0.5182 0.5134 -0.6000 0.5533 0.5488 -0.5333 0.5903 0.5866 -0.4667 0.6297 0.6271 -0.4000 0.6716 0.6703 -0.3333 0.7165 0.7165 -0.2667 0.7646 0.7659 -0.2000 0.8163 0.8187 -0.1333 0.8718 0.8752 -0.0667 0.9314 0.9355 0.0000 0.9956 1.0000 0.0667 1.0645 1.0689 0.1333 1.1386 1.1426 0.2000 1.2181 1.2214 0.2667 1.3034 1.3056 0.3333 1.3947 1.3956 0.4000 1.4924 1.4918 0.4667 1.5968 1.5947 0.5333 1.7081 1.7046 0.6000 1.8268 1.8221 0.6667 1.9532 1.9477 0.7333 2.0874 2.0820 0.8000 2.2300 2.2255 0.8667 2.3811 2.3790 0.9333 2.5411 2.5430 1.0000 2.7103 2.7183 pcubic_test Test pcubic() 0 2.51 0.5 2.91 1 3.3 1.5 3.68 2 4.04 2.5 4.38 3 4.7 3.5 4.98 4 5.22 4.5 5.4 5 5.54 5.5 5.67 6 5.78 6.5 5.44 7 5.4 7.5 5.47 8 5.57 8.5 5.65 9 5.7 9.5 5.74 10 5.78 10.5 5.81 11 5.83 11.5 5.84 12 5.84 12.5 5.83 13 5.81 13.5 5.79 14 5.75 14.5 5.71 15 5.65 15.5 5.57 16 5.48 16.5 5.36 17 5.22 17.5 5.07 18 4.9 18.5 4.72 19 4.53 19.5 4.34 20 4.16 Graphics saved as "pcubic_test.png" polint_test: Test polint() Graphics saved as "polint_test.png" note that the error has a double zero at 4 and a simple zero at 5 and appears to have a triple zero at 2, as expected rgfls_test Test rgfls() which uses the regula falsi method to seek a solution of a nonlinear equation. X = 2.000000e+00, F(X) = 5.000000e+00 X = 3.000000e+00, F(X) = 2.300000e+01 F is of the same sign at the two endpoints 2.000000e+00, 3.000000e+00 The zero is NaN within [2.000000e+00 .. 3.000000e+00], IFLAG = -1 X = 1.000000e+00, F(X) = -1.000000e+00 X = 2.000000e+00, F(X) = 5.000000e+00 X = 1.166667e+00, F(X) = -5.787037e-01 X = 1.253112e+00, F(X) = -2.853630e-01 X = 1.293437e+00, F(X) = -1.295421e-01 X = 1.311281e+00, F(X) = -5.658849e-02 X = 1.318989e+00, F(X) = -2.430375e-02 X = 1.322283e+00, F(X) = -1.036185e-02 X = 1.323684e+00, F(X) = -4.403950e-03 X = 1.324279e+00, F(X) = -1.869258e-03 X = 1.324532e+00, F(X) = -7.929592e-04 X = 1.324639e+00, F(X) = -3.363010e-04 no convergence in 10 iterations The zero is 1.324639e+00 within [1.324639e+00 .. 2.000000e+00], IFLAG = 2 X = 1.000000e+00, F(X) = -1.000000e+00 X = 2.000000e+00, F(X) = 5.000000e+00 X = 1.166667e+00, F(X) = -5.787037e-01 X = 1.253112e+00, F(X) = -2.853630e-01 X = 1.293437e+00, F(X) = -1.295421e-01 X = 1.311281e+00, F(X) = -5.658849e-02 X = 1.318989e+00, F(X) = -2.430375e-02 X = 1.322283e+00, F(X) = -1.036185e-02 X = 1.323684e+00, F(X) = -4.403950e-03 X = 1.324279e+00, F(X) = -1.869258e-03 X = 1.324532e+00, F(X) = -7.929592e-04 X = 1.324639e+00, F(X) = -3.363010e-04 no convergence in 10 iterations The zero is 1.324639e+00 within [1.324639e+00 .. 2.000000e+00], IFLAG = 2 X = 1.000000e+00, F(X) = -1.000000e+00 X = 2.000000e+00, F(X) = 5.000000e+00 X = 1.166667e+00, F(X) = -5.787037e-01 X = 1.253112e+00, F(X) = -2.853630e-01 X = 1.293437e+00, F(X) = -1.295421e-01 X = 1.311281e+00, F(X) = -5.658849e-02 no convergence in 4 iterations The zero is 1.311281e+00 within [1.311281e+00 .. 2.000000e+00], IFLAG = 2 rk2_test Test rk2() Demonstrate second order Runge Kutta method. Method 2 -- Simplified Runge-Kutta order 2 H Y(1.5) YPRIM(1.5) Y(2) YPRIM(2) 6.25000000e-02 -6.65527271e-01 4.45202743e-01 -4.98224194e-01 2.50884750e-01 3.12500000e-02 -6.66377044e-01 4.44637442e-01 -4.99548659e-01 2.50225467e-01 1.56250000e-02 -6.66593713e-01 4.44493075e-01 -4.99886330e-01 2.50056822e-01 7.81250000e-03 -6.66648363e-01 4.44456647e-01 -4.99971484e-01 2.50014257e-01 spline_cd_test Test spline_cd() 0 2.51 0.5 2.91 1 3.3 1.5 3.68 2 4.04 2.5 4.38 3 4.7 3.5 4.98 4 5.22 4.5 5.4 5 5.54 5.5 5.67 6 5.78 6.5 5.44 7 5.4 7.5 5.47 8 5.57 8.5 5.65 9 5.7 9.5 5.74 10 5.78 10.5 5.81 11 5.83 11.5 5.84 12 5.84 12.5 5.83 13 5.81 13.5 5.79 14 5.75 14.5 5.71 15 5.65 15.5 5.57 16 5.48 16.5 5.36 17 5.22 17.5 5.07 18 4.9 18.5 4.72 19 4.53 19.5 4.34 20 4.16 Graphics saved as "spline_cd_test.png" subst_test: Test subst() W = 1 1 1 0 1 1 0 0 1 IPIVOT = 1 2 3 IFLAG = 1 ans = -1.1102e-16 1.1102e-16 0.0000e+00 W = 1.00000 1.00000 0.00000 0.50000 -0.50000 0.00000 1.00000 -0.00000 1.00000 IPIVOT = 2 1 3 IFLAG = -1 ans = 0.0000e+00 0.0000e+00 -5.5511e-17 W = 1 0 0 1 1 0 0 1 1 IPIVOT = 1 2 3 IFLAG = 1 ans = 0.0000e+00 0.0000e+00 5.5511e-17 W = 0.76705 0.83613 0.40900 0.11626 0.80676 0.20760 0.69684 0.17232 0.32101 IPIVOT = 2 3 1 IFLAG = 1 ans = 0.0000e+00 -2.2204e-16 3.3307e-16 W = 1 0 0 0 1 0 0 0 1 IPIVOT = 3 1 2 IFLAG = 1 ans = 0 0 0 table_cd_test: Test table_cd() KP1 = 2, XK(KP1) = 0.700000, FXBAR = 0.564642, ERROR = 0.023873 KP1 = 3, XK(KP1) = 0.500000, FXBAR = 0.588515, ERROR = 0.000592 value = 5.8910742e-01, IFLAG = 1, error = -3.7339695e-05 KP1 = 2, XK(KP1) = 0.700000, FXBAR = 0.564642, ERROR = 0.023873 KP1 = 3, XK(KP1) = 0.500000, FXBAR = 0.588515, ERROR = 0.000592 KP1 = 4, XK(KP1) = 0.800000, FXBAR = 0.589107, ERROR = 0.000036 KP1 = 5, XK(KP1) = 0.400000, FXBAR = 0.589144, ERROR = 0.000001 KP1 = 6, XK(KP1) = 0.900000, FXBAR = 0.589145, ERROR = 0.000000 KP1 = 7, XK(KP1) = 0.300000, FXBAR = 0.589145, ERROR = 0.000000 KP1 = 8, XK(KP1) = 1.000000, FXBAR = 0.589145, ERROR = 0.000000 KP1 = 9, XK(KP1) = 0.200000, FXBAR = 0.589145, ERROR = 0.000000 KP1 = 10, XK(KP1) = 0.100000, FXBAR = 0.589145, ERROR = 0.000000 KP1 = 11, XK(KP1) = 0.000000, FXBAR = 0.589145, ERROR = -0.000000 KP1 = 12, XK(KP1) = -0.100000, FXBAR = 0.589145, ERROR = -0.000000 KP1 = 13, XK(KP1) = -0.200000, FXBAR = 0.589145, ERROR = 0.000000 KP1 = 14, XK(KP1) = -0.300000, FXBAR = 0.589145, ERROR = 0.000000 KP1 = 15, XK(KP1) = -0.400000, FXBAR = 0.589145, ERROR = -0.000000 value = 5.8914476e-01, IFLAG = 1, error = 1.1102230e-16 1.200000 not in table range value = 8.4147098e-01, IFLAG = 3, error = -9.0568101e-02 trap_test: Test trap() N Integral estimate 2 16.1606490872 4 16.1751521298 8 30.0632333826 16 28.9120538071 32 29.3062382835 64 29.3209883235 128 29.3249060238 256 29.3258867692 512 29.32613204 1024 29.3261933629 trap_test: Normal end of execution. example2p4: N MAXIMUM ERROR 2 6.4615385e-01 4 4.3813387e-01 6 6.1666759e-01 8 1.0451739e+00 10 1.9156431e+00 12 3.6052745e+00 14 7.1920080e+00 16 1.4051542e+01 example3p11 Newton's method for finding a real zero of a polynomial I X AP APP 1 1.0000000e+00 -2.1000000e+00 1.6000000e+00 2 2.3125000e+00 1.2253150e+01 3.9531570e+01 3 2.0025414e+00 3.6459060e+00 1.8126400e+01 4 1.8014035e+00 8.7345726e-01 1.0066133e+01 5 1.7146316e+00 1.0919285e-01 7.6418992e+00 6 1.7003429e+00 2.5000177e-03 7.2942739e+00 7 1.7000002e+00 1.4001847e-06 7.2861046e+00 8 1.7000000e+00 4.3787196e-13 7.2861000e+00 conte_deboor_test(): Normal end of execution. 07-Mar-2024 11:44:56