e73fbdeb-2ad8-477d-9a56-dbbd0f030aba 9 4 5 decimal

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

Setup a parametric ODE system

...  an invented ODE demo system  ... t X t diff k 1 el t X * - ≡ t Y t diff k 1 el t X * k 2 el t Y * - k 3 el t Z * + ≡ t Z t diff k 2 el t Y * k 3 el t Z * - ≡ 31 mat

... the executable numerical ODE ... t Y k D k 1 el Y 1 el * - k 1 el Y 1 el * k 2 el Y 2 el * - k 3 el Y 3 el * + k 2 el Y 2 el * k 3 el Y 3 el * - 31 mat :

Smath executable ... fill-in  [IC, k] sol # # # 31 mat t0 t1 N t Y # # # 31 mat D Rkadapt ≡

fill-in the demand t0 0 : t1 5 : N 200 : 13 mat

Ω, Y, U ...as you wish

sol 21031 mat t0 t1 N t Ω 12331 mat D Rkadapt :

header t0 IC1 IC2 IC3 14 mat :

IC's

parameters

sol header sol stack :

sol t0 IC1 IC2 IC3 0210 0.02 1.95 1 0.05 0.05 1.9 1 0.09 0.08 1.86 1.01 0.13 0.1 1.81 1.02 0.17 0.12 1.76 1.02 0.21 0.15 1.72 1.03 0.24 0.18 1.68 1.04 0.28 0.2 1.64 1.05 0.31 0.22 1.6 1.07 0.34 0.25 1.56 1.08 0.36 0.28 1.52 1.09 0.39 0.3 1.48 1.1 0.41 0.32 1.45 1.12 0.44 0.35 1.41 1.13 0.46 0.38 1.37 1.14 0.48 0.4 1.34 1.16 0.5 0.42 1.31 1.17 0.52 0.45 1.28 1.18 0.54 0.48 1.24 1.2 0.56 0.5 1.21 1.21 0.58 0.52 1.18 1.22 0.59 234 mat

sol sol 2 sol rows range 14 range el :

sol1 sol 1 col sol 2 col augment : sol2 sol 1 col sol 3 col augment : sol3 sol 1 col sol 4 col augment : 31 sys

vx sol 1 col : vy sol 4 col : 12 mat x f vx vy x cinterp : 21 sys

1.23456789 f 0.9086

sol1 sol2 sol3 3 1 sys

NOT sanity checked x f x 02 int 1.4953

sol rows 201