ef341626-30a8-473d-a2e6-523d781c30fe 17 6 9 decimal

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

t x D μ x 1 el * x 2 el - x 1 el (2^x 2 el (2^+ (x 1 el * - μ x 2 el * x 1 el + x 1 el (2^x 2 el (2^+ (x 2 el * - 21 mat :

μ 0.2 - : t0 0 : t1 20 : N 1000 : 41 line

IC 0121 mat :

Sol IC t0 t1 N t x D Rkadapt :

x Sol 1 col : X Sol 2 col : Y Sol 3 col : 31 line

U x X augment : V x Y augment : Spiral Y X augment : 31 line

Spiral U V 2 1 sys

=============================================================================

α t y t 2 diff * k t y t diff (* + β t y * + p ω t * cos * ≡

t y D y 2 el p ω t * cos * k y 2 el * - β y 1 el * - α / 21 mat :

β 0.03 : p 0.25 : ω 0.9 : α 1.5 : k 1.5 : 51 line

t0 0 : t1 20 : N 400 : 31 line

IC 0.5 - 121 mat :

Either 'rkfixed' or 'Rkadapt' s IC 020400 t y D Rkadapt :

u s 1 col : XX s 2 col : YY s 3 col : 31 line

vx u XX augment : vy u YY augment : 21 line

vx vy 2 1 sys

k 0.025 : g 9.807 - : 21 line

t0 0 : t1 10 : N 200 : 31 line

t y f y 1 el g - k y 2 el * y 1 el (2^y 2 el (2^+ sqrt * y 1 el 1 + (* - :

init 0 0.125 21 mat :

sol init t0 t1 N t y f Rkadapt :

Smath bug in there. Does not recognise this style DE f(t,y)

x sol 1 col : Y sol 2 col : V x Y augment : 31 line