c7fc7f15-cc26-413e-b246-a305299d5cc6 4 4 5 decimal

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

g 9.81 :

k 80 :

m 8 :

ω k m / sqrt :

L 0.5 :

IC 41 matrix :

IC 1 el 1 :

IC 2 el 2 :

t φ D φ 3 el φ 4 el L φ 4 el t diff (* φ 2 el sin * L φ 4 el 2^* φ 2 el cos * + g + ω 2^φ 1 el * - φ 3 el t diff (φ 2 el sin * g φ 2 el sin * - L / 41 mat :

φ 123441 mat :

n t.1 Δt / :

t.1 3 :

Δt 0.01 :

S IC 0 t.1 n t φ D mwode45 :

S 01200 0.01 1.03 2.04 0.0352 - 0.1784 - 0.02 1.06 2.08 0.0704 - 0.3568 - 0.03 1.09 2.12 0.1056 - 0.5352 - 0.04 1.12 2.16 0.1408 - 0.7136 - 0.05 1.15 2.2 0.176 - 0.892 - 0.06 1.18 2.24 0.2112 - 1.0704 - 0.07 1.21 2.28 0.2463 - 1.2488 - 0.08 1.24 2.32 0.2815 - 1.4272 - 0.09 1.27 2.36 0.3167 - 1.6056 - 0.1 1.3 2.4 0.3519 - 1.784 - 0.11 1.33 2.44 0.3871 - 1.9624 - 0.12 1.36 2.48 0.4223 - 2.1408 - 0.13 1.39 2.52 0.4575 - 2.3193 - 0.14 1.42 2.56 0.4927 - 2.4977 - 0.15 1.45 2.6 0.5279 - 2.6761 - 0.16 1.48 2.64 0.5631 - 2.8545 - 0.17 1.51 2.68 0.5983 - 3.0329 - 0.18 1.54 2.72 0.6335 - 3.2113 - 0.19 1.57 2.76 0.6686 - 3.3897 - 0.2 1.6 2.8 0.7038 - 3.5681 - 0.21 1.63 2.84 0.739 - 3.7465 - 0.22 1.66 2.88 0.7742 - 3.9249 - 235 mat

S 1 col S 4 col augment

τ 0 n range :

t θ S 3 col t 1 + el :

t λ S 2 col t 1 + el :

t x L t θ sin * eval :

t y t λ L t θ cos * + eval :

0 0 o 10 1 4 mat 0 0 + 10 1 4 mat 0 t λ - . 15 1 4 mat t x - t y - . 20 1 4 mat 0 0 0 t λ - 2 2 mat 0 t λ - t x - t y - 2 2 mat t 0 ≡ path t x - t y - 1 2 mat : path path t x - t y - 1 2 mat stack : if 7 1 sys