<p> ∫Adaptive(f,a,b,N,ε) </p>

2db2e12d-7b48-4733-8f0a-069c4243b8ca 3 15 12 decimal

&[DATE] &[TIME] - &[FILENAME]

f a b N ε ∫Adaptive sum 0 : i 1 : tol i el ε : aa i el a : h i el b a - (2 / : 51 line .... fa i el a f : fc i el a h i el + f : fb i el b f : s i el h i el fa i el 4 fc i el * + fb i el + (3 / * : L i el 1 : 51 line 13 mat above => initialise the while loop i 0 > fd aa i el h i el 2 / + f : fe aa i el 3 h i el * 2 / + f : 12 mat s1 h i el fa i el 4 fd * + fc i el + 6 / (* : s2 h i el fc i el 4 fe * + fb i el + 6 / (* : 12 mat v1 aa i el : v2 fa i el : v3 fc i el : v4 fb i el : v5 h i el : v6 ε : v7 s i el : v8 L i el : 18 mat i i 1 - : s1 s2 + v7 - abs v6 < sum sum s1 s2 + (+ : v8 N ≥ Solution Insufficient Number of Intervals : i i 1 + : aa i el v1 v5 + : fa i el v3 : fc i el fe : fb i el v4 : 14 mat h i el v5 2 / : tol i el v6 2 / : s i el s2 : L i el v8 1 + : 14 mat i i 1 + : aa i el v1 : fa i el v2 : fc i el fd : fb i el v3 : h i el h i 1 - el : s i el s1 : 16 mat L i el L i 1 - el : 61 line if if 51 line while sum 41 line eval :

γ 0.5772156649015 :

x c 11 - x ≤ cases :

z digamma γ - z n n z + (* / (n 120 sum + 11 line :

0.5 - digamma 1.93882242848689 -

x digamma x c *

x f x digamma :

'N' iterations
precision N 1000 : TOL 106 -^: 21 sys

User inputs

lower limit
upper limit a 1 - 1012 -^+ : b e : 21 sys

t0 1 time :

Sanity Mathcad .....   -26.33 f a b N TOL ∫Adaptive 26.33 -

1 time t0 - 28 s *