1f24c4a6-6ee9-45ca-9c7b-15b880367196 22 3 5 decimal

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

Smath ODE librairies

Matlab solvers

Intel solvers

Useless solvers. # # # # # al_rkckadapt # # # # # dn_AdamsMoulton # # # # # dn_ExplicitRK45 # # # # # dn_GearsBDF # # # # # dn_ImplicitRK5 51 line

ONE  possible limited solver. # # # # # mwode113 # # # # # mwode15s # # # # # mwode23 # # # # # mwode23s # # # # # mwode45 51 line

# # # # # Rkadapt # # # # # rkfixed # # # # # rkm9mka # # # # # rkm9mkn # # # # # rkm9st 51 line

# # # # # mk52lfa # # # # # mk52lfn 21 line

x f x x 2^cos * 2 x * - 3 + :

data spec pts i 1 data rows range p i el a b c spec : for v p 1 el : Unest 'p' j 2 p rows range v v p j el stack : for data v augment 51 line :

t0 1 time :

a b 12 mat 1512 mat :

N 200 :

x y D x f :

M 011 mat a b N x y D rkm9mka :

1 time t0 - 0.5 s *

prepare the project
items are invariant vx M 1 col : vy M 2 col : 12 mat t q vx vy t cinterp : 21 sys

user population δ 0.1 : X a b a δ + (range : 12 mat

t0 1 time :

data i 1 X rows range Ω i el X i el q eval : for X Ω augment 21 line :

user points spec
... NUMERIC ... a b c spec o 6 black 13 mat : pts data a b c spec pts : 21 sys

x r 1 a x ≤ (x b ≤ (& cases :

pts M 2 1 sys

1 time t0 - 0.7 s *

crappy ∫(z) over cinterp
plots on canvas: very slow t I u q u a t int :

pts M 21 sys

u f u a b int 12.48691 -

5 I 14.907 -

M M rows 2 el 12.48691 -

x map x f x int maple : x 2^sin 2 x 2^* - 6 x * + 2 /

Observe fidelity of solver(s)

x map 1 map - x q - ( x r * x map 1 map - ( x r * pts 2 1 sys

Parametric application . . . that you may have to disambiguate

]]>

t μ α D α 1 el μ 2 el μ 1 el - (* μ 1 el α 3 el μ 3 el - (* μ 1 el - μ 1 el μ 2 el * α 2 el μ 3 el * - 31 mat :

init 0.1 0.1 0.1 31 mat :

Example collaboration Smath Martin, Viacheslav, Radovan

x y XY x y / - :

t.0 0 :

t.max 50 :

n 2000 :

adapt init t.0 t.max n t μ 2 1.5 2031 mat D Rkadapt :

X adapt 1 col : Y adapt 2 col : Z adapt 3 col : 31 line

user init t.0 t.max n t α 2 1.5 2031 mat D dn_GearsBDF :

U user 1 col : V user 2 col : W user 3 col : 31 line

Y Z augment XY 2 1 sys V W augment XY 2 1 sys