<p>RK23</p>

7446336e-14b1-4edd-aa79-bfdb9158050a 5 15 12 decimal

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

m 140 :

g 9.8 :

V 16 :

p 1.42 :

I.xG 6.7 :

h.G 0.40 :

I.zG 36.4 :

b.G 0.703 :

NOT needed t Clear 1

m V * ψ' * b.G V ψ' / / atan sin p V ψ' / / atan sin / * p * (m V * ψ' * b.G * b.G V ψ' / / atan sin * (- I.zG m b.G * + (ψ'' * (- 0 ≡

t Ψ D Ψ 2 el m V * Ψ 1 el * b.G V Ψ 1 el / / atan sin p V Ψ 1 el / / atan sin / * p * m V * Ψ 1 el * b.G * b.G V Ψ 1 el / / atan sin * - I.zG m b.G 2^* + / 21 mat :

t.o t.end 12 mat 0912 mat :

ψ'.o ψ''.o 12 mat 0012 mat :

IC 0.1 021 mat :

you can switch from  Ψ  to U ψsol IC 01100 t U D dn_GearsBDF :

t ψ' ψsol 1 col ψsol 2 col t cinterp :

t ψ'' ψsol 1 col ψsol 3 col t cinterp :

x ψ''

disabled

Runge - Kutta 2° e 3° ordine

2 D\032E\ to\032E\ xo\032E\ N\032E\ E\032E\ RK23 x\032E\ xu\032E\ UIn x\032E\ xu\032E\ UnitsOf * : x\032E\ UOut x\032E\ x\032E\ UnitsOf / vectorize : 12 mat k\0324\ 1 xo\032E\ length range : x\0324\ xo\032E\ k\0324\ el : X\0324\ x\0324\ transpose : 13 mat n\0324\ 1 : t\0324\ to\032E\ 1 el : T\0324\ t\0324\ 11 mat : h\0324\ to\032E\ 2 el t\0324\ - N\032E\ / eval : ho\0324\ h\0324\ : 15 mat t\0324\ x\0324\ k\0324\ h\0324\ K\0324\ h\0324\ t\0324\ x\0324\ D\032E\ * eval h\0324\ E\032E\ ho\0324\ UIn x\0324\ E\032E\ x\0324\ UIn + D\032E\ * eval try 11 line : t\0324\ to\032E\ 2 el ≤ k\0324\.1 t\0324\ x\0324\ x\0324\ h\0324\ K\0324\ : \0324\k.2 t\0324\ h\0324\ + x\0324\ k\0324\.1 + k\0324\.1 h\0324\ K\0324\ : \0324\k.3 t\0324\ 0.5 h\0324\ * + x\0324\ 0.25 k\0324\.1 * + 0.25 \0324\k.2 * + \0324\k.2 h\0324\ K\0324\ : t\0324\ trace d\0324\ k\0324\.1 2 \0324\k.3 * - \0324\k.2 + 3 / UOut normi : \0324\q E\032E\ x\0324\ UOut normi ho\0324\ UOut 12 mat max * : d\0324\ \0324\q ≤ t\0324\ t\0324\ h\0324\ + : T\0324\ n\0324\ n\0324\ 1 + : el t\0324\ : 12 mat x\0324\ x\0324\ k\0324\.1 4 \0324\k.3 * + \0324\k.2 + 6 / + eval : X\0324\ n\0324\ k\0324\ el x\0324\ k\0324\ el : 12 mat 21 line if d\0324\ 0 ≠ h\0324\ ho\0324\ 0.9 h\0324\ * \0324\q d\0324\ / 3 nthroot * 12 mat min eval : if 81 line while T\0324\ X\0324\ augment eval 61 line :