564cf393-9f74-42a7-b1c8-74cabbc75291 74 3 5 decimal

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

t0 1 time :

dn_GearsBDF is 1/1 MCD The MCD solver "Radau" introduced Mathcad 2001i

http://oai.cwi.nl/oai/asset/4619/04619D.pdf

t Y f 7100171 Y 1 el * - 43 Y 2 el * + 832 Y 3 el * + (* + 10000 / 171 Y 1 el * 875 Y 2 el * - 100 / 21003 Y 3 el * - 43 Y 4 el * + (* 7 Y 5 el * + 200 / 1652 Y 2 el * 7 Y 4 el * - (* 171 Y 3 el * + 100 / 349 Y 5 el * - 86 Y 6 el Y 7 el + (* + 200 / Y 6 el 4328000 Y 8 el * + (* - 323 Y 4 el * 57 Y 5 el * + 23 Y 7 el * + (* + 100 / 28000 Y 6 el * Y 8 el * 181 Y 7 el * - 100 / 28000 Y 6 el * Y 8 el * - 181 Y 7 el * + 100 / 81 mat :

IC 1000000 0.0057 81 mat :

3 possible robust solvers. # # # # # dn_AdamsMoulton 3 min ≡ # # # # # dn_ExplicitRK45 2 min ≡ # # # # # dn_GearsBDF 10 sec ≡ 31 line

5 possible robust solvers. # # # # # mwode113 2.5 min ≡ # # # # # mwode15s 40 sec ≡ # # # # # mwode23 2.5 min ≡ # # # # # mwode23s 50sec ≡ # # # # # mwode45 3 min ≡ 51 line

sol IC 03221000 t Y f dn_GearsBDF :

X sol 1 col : C sol 6 col : D sol 7 col : 31 line

X C augment X D augment 2 1 sys

Rkadapt: adequate for these 3 pairs reactions sol IC 0101000 t Y f dn_GearsBDF :

X sol 1 col : U sol 2 col : V sol 5 col : W sol 3 col : Z sol 4 col : A sol 8 col : B sol 9 col : 71 line

X U augment X V augment 2 1 sys X W augment X Z augment 2 1 sys X A augment X B augment 2 1 sys

1 time t0 - 10.5 s *

n x L 1 - (k^2 n * 2 k * - (! (x n 2 k * -^* (* 2 n^k! (* n k - (! (n 2 k * - (! (* (* / (k 0 n 2 / Floor sum 11 line :

n x Radau n 1 - x L n x L + x 1 + / :

x y axes x y / - :

ε 1012 -^:

5 x L 6 x L axes 3 1 sys 5 x ε - Radau axes 2 1 sys

k 16 range Legendre k el k x L : for

k 15 range Radau k el k x Radau : for

Legendre x 1 - 3 x 2^* + 2 / x 3 - 5 x 2^* + (* 2 / 35 x 2^* 6 - 7 x 2^* + (* + 8 / x 157 x 2^* 10 - 9 x 2^* + (* + (* 8 / 5 - 21 x 2^* 5 x 2^15 - 11 x 2^* + (* + (* + 16 / 61 mat

Radau 11 - x 23 x * + (* + 21 x + (* / 1 x 2^35 x * + (* - 3 x * + 21 x + (* / - 3 x 223 - 5 x 2^* + (* 15 x * - (* 35 x 3^* + (* + 81 x + (* / 3121 x 5^* + (* 5 x 2^* 67 x 2^* - 14 x * + (* - 15 x * + 81 x + (* / 51 mat