bbd6f987-0d47-4d9f-8a42-2fb56c0a2e74 21 6 9 decimal

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

Jacobian stiff demo

x y z f x y * x 3^z * - x y * sin z * 21 mat :

f x y z JacobiStiff x y z f 1 el (x diff x y z f 1 el (y diff x y z f 1 el (z diff x y z f 2 el (x diff x y z f 2 el (y diff x y z f 2 el (z diff 23 mat :

Expand symbolic f x y z JacobiStiff y - 3 x 2^* z * + (- x x 3^- y x y * cos * z * x x y * cos * z * x y * sin 23 mat

Substitute by local definitions t Y J x t : y Y 1 el : z Y 2 el : f x y z JacobiStiff 41 line :

t Y J 3 t 2^* Y 2 el * Y 1 el - (- t t 3^- Y 1 el t Y 1 el * cos * Y 2 el * t t Y 1 el * cos * Y 2 el * t Y 1 el * sin 23 mat

t0 0 : t1 3 : N 200 : 31 line

IC 5021 mat :

S IC t0 t1 N t Y J Rkadapt :

t S 1 col : U S 2 col : V S 3 col : 31 line

SOL_1 t U augment : SOL_2 t V augment : 21 line

SOL_1 SOL_2 2 1 sys

step i 1 U rows range U i el 32 ≤ U i el U i el 0 : if for t U augment 21 line :

step SOL_2 2 1 sys