Accuracy of Integral with epsilon

10092021, 02:01 PM
Post: #15




RE: Accuracy of Integral with epsilon
(10092021 01:41 AM)robve Wrote:(10082021 11:28 PM)lrdheat Wrote: Did not think the integral was remarkable, but was fascinated by how different calculators handled it as written. My Prime in home shows 40.0000000143 (!). Impressed with the WP 34S. WP 34S results with a IEEE 754 based implementation matching the Python code given by the author: eps=10^4: abs error = 0.0051939; est. error = 0.0155524; evals = 7 eps=10^5: abs error = 0.000130107; est. error = 0.0156835; evals = 9 eps=10^6: abs error = 0.000121986; est. error = 6.30233e06; evals = 17 eps=10^7: abs error = 2.14894e06; est. error = 3.30639e06; evals = 19 eps=10^8: abs error = 2.14894e06; est. error = 3.30639e06; evals = 19 eps=10^9: abs error = 5.81537e09; est. error = 3.35997e06; evals = 21 Note that the given eps is not necessarily smaller than the abs error obtained ("under delivers"). Only the last result delivers a decent estimate. Whether or not the convergence check is satisfied with the estimated error is a detail that may differ in the implementation(s) of TanhSinh, thus eps=10^4 may produce a better answer for some, but to get the better answer 19, 21 or more evaluations are needed to get close enough. The qthsh implementation (discussed on the HP forum thread) of the same algorithm with some differences has 27 evaluations for eps=10^4 with small abs error 5.63299e09. Qthsh approaches integration boundaries more consistently, evaluating additional boundary points by handling singularities instead of terminating early. Otherwise the algorithms are the same. This explains why qthsh has 27 versus WP 34S 21 evaluations. Like the WP 34S, qthsh was verified and compared to WP 34S and other algorithms, including a verification of accuracy with over 800 integrals.  Rob "I count on old friends"  HP 71B,PrimeTi VOY200,Nspire CXII CASCasio fxCG50...Sharp PCG850,E500,2500,1500,14xx,13xx,12xx... 

« Next Oldest  Next Newest »

User(s) browsing this thread: 1 Guest(s)