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

abc83e55-fbeb-499c-8062-f84451b71c8a 9 15 12 decimal

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

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 :

α 1523043253 - 27354.7 - 43256.75 25388.3 - 40061 mat :

Fréchet: extreme value-II x f α 1 el α 2 el x α 4 el - (/ (α 5 el^* α 2 el x α 4 el - (/ (α 3 el^- exp * α 6 el + :

end 12 :

x c 10 x ≤ (x end ≤ (& cases :

0 f 403

x Y D x f :

init 0 :

sol init 0 end 100 x Y D dn_GearsBDF :

in-situ Simpson acc

x Frechet f 0 x 10 Int :

t0 1 time :

x f x c * sol 2 1 sys x f x c * x Frechet x c * 2 1 sys

1 time t0 - 19 s *

data X 0121 range : i 1 X rows range F i el X i el Frechet : for X F augment 31 line :

data 0014702186035727411494517110621429724388826339927642102855611292451229807132 mat

Result is above max 
allowed positive number x f x 0 end int #

12 Frechet 29807

sol sol rows row 122978612 mat

=============================================================

α 1.1077 142.40751 1.01751 - 0.07235 13.58616 - 51 mat :

x α Φ x α 1 el α 2 el x α 4 el / - exp * + α 3 el x α 5 el / - exp * + (* :

x c 10 x ≤ (x 1.25 ≤ (& cases :

x Y D x α Φ :

init 0 :

end 1.25 :

sol init 0 end 100 x Y D dn_GearsBDF :

sol sol rows row 1.25 0.766 12 mat

x α Φ x 0 1.25 int 0.765

x α Φ x c * sol 2 1 sys