Tue Oct 19 11:57:22 2021 line_ncc_rule_test(): Python version: 3.6.9 line_ncc_rule() computes closed Newton Cotes quadrature rules for a line segment. line_ncc_rule_test01(): line_ncc_rule() computes the closed Newton-Cotes rules using N equally spaced points for an interval [A,B]. Closed Newton-Cotes Rule 1 I X(I) W(I) 0 -1 2 Sum(|W)|) = 2.0 Closed Newton-Cotes Rule 2 I X(I) W(I) 0 -1 1 1 1 1 Sum(|W)|) = 2.0 Closed Newton-Cotes Rule 3 I X(I) W(I) 0 -1 0.333333 1 0 1.33333 2 1 0.333333 Sum(|W)|) = 2.0 Closed Newton-Cotes Rule 4 I X(I) W(I) 0 -1 0.25 1 -0.333333 0.75 2 0.333333 0.75 3 1 0.25 Sum(|W)|) = 1.9999999999999996 Closed Newton-Cotes Rule 5 I X(I) W(I) 0 -1 0.155556 1 -0.5 0.711111 2 0 0.266667 3 0.5 0.711111 4 1 0.155556 Sum(|W)|) = 2.0 Closed Newton-Cotes Rule 6 I X(I) W(I) 0 -1 0.131944 1 -0.6 0.520833 2 -0.2 0.347222 3 0.2 0.347222 4 0.6 0.520833 5 1 0.131944 Sum(|W)|) = 2.0000000000000044 Closed Newton-Cotes Rule 7 I X(I) W(I) 0 -1 0.097619 1 -0.666667 0.514286 2 -0.333333 0.0642857 3 0 0.647619 4 0.333333 0.0642857 5 0.666667 0.514286 6 1 0.097619 Sum(|W)|) = 2.0000000000000058 Closed Newton-Cotes Rule 8 I X(I) W(I) 0 -1 0.0869213 1 -0.714286 0.414005 2 -0.428571 0.153125 3 -0.142857 0.345949 4 0.142857 0.345949 5 0.428571 0.153125 6 0.714286 0.414005 7 1 0.0869213 Sum(|W)|) = 2.000000000000016 Closed Newton-Cotes Rule 9 I X(I) W(I) 0 -1 0.0697707 1 -0.75 0.415379 2 -0.5 -0.0654674 3 -0.25 0.740459 4 0 -0.320282 5 0.25 0.740459 6 0.5 -0.0654674 7 0.75 0.415379 8 1 0.0697707 Sum(|W)|) = 2.902433862433886 Closed Newton-Cotes Rule 10 I X(I) W(I) 0 -1 0.0637723 1 -0.777778 0.351362 2 -0.555556 0.0241071 3 -0.333333 0.431786 4 -0.111111 0.128973 5 0.111111 0.128973 6 0.333333 0.431786 7 0.555556 0.0241071 8 0.777778 0.351362 9 1 0.0637723 Sum(|W)|) = 1.999999999999812 Closed Newton-Cotes Rule 11 I X(I) W(I) 0 -1 0.0536683 1 -0.8 0.355072 2 -0.6 -0.162087 3 -0.4 0.909893 4 -0.2 -0.87031 5 0 1.42753 6 0.2 -0.87031 7 0.4 0.909893 8 0.6 -0.162087 9 0.8 0.355072 10 1 0.0536683 Sum(|W)|) = 6.129589546255551 Closed Newton-Cotes Rule 12 I X(I) W(I) 0 -1 0.0498665 1 -0.818182 0.309711 2 -0.636364 -0.0743385 3 -0.454545 0.579317 4 -0.272727 -0.220356 5 -0.0909091 0.355801 6 0.0909091 0.355801 7 0.272727 -0.220356 8 0.454545 0.579317 9 0.636364 -0.0743385 10 0.818182 0.309711 11 1 0.0498665 Sum(|W)|) = 3.1787785677584375 line_ncc_rule_test02(): Use a sequence of Closed Newton Cotes rules for an estimate Q of the integral: I = integral ( 0 <= x <= 1 ) exp(x) dx. The exact value is: I = 1.718281828459045 N Q |Q-I| 1 1 0.718282 2 1.85914 0.140859 3 1.71886 0.000579323 4 1.71854 0.000258325 5 1.71828 8.59466e-07 6 1.71828 4.84531e-07 7 1.71828 1.05866e-09 8 1.71828 6.47825e-10 9 1.71828 1.50013e-11 10 1.71828 9.1511e-11 11 1.71828 4.75218e-10 12 1.71828 1.41685e-09 13 1.71828 1.31693e-08 14 1.71828 2.22914e-08 15 1.71828 1.33324e-06 16 1.71828 2.23641e-06 17 1.71825 2.68474e-05 18 1.7183 1.40187e-05 19 1.72013 0.00184715 20 1.72148 0.00319798 21 1.76478 0.0464972 22 2.16025 0.441972 line_ncc_rule_test(): Normal end of execution. Tue Oct 19 11:57:22 2021