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

f2d5c58c-29a8-4f20-8e0e-92a808230cce 8 15 12 decimal

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

x Γ x Gamma :

α 1.063276612 :

β 190.68821559 :

12 / Γ 1.772

1 β α^α Γ * / 0.003889

α Γ 0.967

x g 1 β α^α Γ * / x α 1 -^e x - β /^* (* :

x Q u g u 0 x int :

x g x Q <p>Int(f,a,b,n) ... Simpson/Romberg Kernel [Jean] </p>

Basic Romberg-Simpson ... integration acc from options ... NUMERIC f a b Romberg w s b a - (34 / 2 w * 1 - (* 14 / 2 w * 1 - (3^* - 12 / + (* a + : w G w s f 12 w * 1 - (2^- (* : 32 / b a - (* w G w 01 int * 31 line eval :

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 
Integration accuracy from user command line 'n' ... NUMERIC 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 / : 15 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 :

x f 1 β α^α Γ * / x α 1 -^e x - β /^* (* :

x Q u f u 0 x int :

x cumulative f 0 x 100 Int :

x Y D x f :

init 0 :

sol init 02001000 x Y D rkfixed :

x cumulative

Smath 6179 Simpson [acc 50] sanity rkfixed

x Q sol 2 1 sys