<p> Int(f,a,b,n) ... Simpson/Romberg Kernel </p>

0764a218-6786-4c57-972c-e77745f9db42 2 6 12 decimal

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

ε 1015 -^:

x f x sin x ε + / :

x f <p> Int(f,a,b,n) ... Simpson/Romberg Kernel </p>

Simpson Integrator in the Romberg Transform Domain [0,1] <= [a,b] ...  
... integrates  Smath ainterp(X,Y,x), cinterp(X,Y,x), cubic(x) from  l_p_c_spline f a b n Int w s b a - (34 / 2 w * 1 - (* 14 / 2 w * 1 - (3^* - 12 / + (* a + : w G w s f 12 w * 1 - (2^- (* : α 0 : β 1 : h β α - n / : 14 mat Simpson 13 / h 2 / α G β G + (* h α k h * + G k 1 n 1 - sum * + (* 23 / h * α h 2 / - k h * + G k 1 n sum * + eval : 32 / b a - (* Simpson * 41 line :

π 2 / 1.570796

f 04000040000 Int 1.570788