<p />

3c6aaaaf-4d34-4425-a77c-f2f17890c368 53 4 5 decimal

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

Kоординаты начальных точек

X0 1 el 1 - :

X0 2 el 1 :

X0 3 el 1 - :

Уравнения:

f 1 el x 2 el (2^x 3 el (2^+ 2 - :

f 2 el x 1 el (2^x 2 el (2^+ 2 - :

<p />

angle Ω1 1000 angle cos angle sin - 0 angle sin angle cos 33 mat :

angle Ω2 angle cos 0 angle sin 010 angle sin - 0 angle cos 33 mat :

E 0111 0.4 1.2095 0.7328 1.2095 0.8 1.3581 0.3944 1.3581 1.2 1.4142 0.0049 1.4142 1.6 1.3607 0.3853 - 1.3607 2 1.214 0.7253 - 1.214 2.4 1.0057 0.9943 - 1.0057 2.8 0.7605 1.1923 - 0.7605 3.2 0.494 1.3251 - 0.494 3.6 0.2161 1.3976 - 0.2161 4 0.0279 1.4171 - 0.0279 114 mat :

angle Ω3 angle cos angle sin - 0 angle sin angle cos 000133 mat :

Оси координат XYZ k OX 10013 mat (γ * k col eval : k OY 01013 mat (γ * k col eval : k OZ 00113 mat (γ * k col eval : 41 sys XYZ 01 OX * - 03 OX * - 21 OX * 23 OX * 0021 OY * 23 OY * 01 OY * - 03 OY * - 0021 OZ * 23 OZ * 01 OZ * - 03 OZ * - 82 mat : TextXYZ 00.10 21 OX * 23 OX * X 10 21 OY * 23 OY * y 10 21 OZ * 23 OZ * Z 10 45 mat : 31 sys

α π - 5 / :

φ π 8 / :

ψ 0.85 π * :

γ α Ω1 φ Ω2 * ψ Ω3 * eval :

<p>Dragilev's Method</p>

X0 tmin tmax N D n f rows 1 + : f n el 0 : j 1 : j n ≤ j j 1 + : k 1 : k n ≤ k k 1 + : m j k el f j el x k el diff : for for Матрица Якоби 21 line x Jacobi m : Матрица  правых частей алгебраических уравнений b 1 - (x Jacobi 1 n 1 - n n submatrix * : Определители A x Jacobi 1 n 1 - 1 n 1 - submatrix : n 2 > Δ 1 el b A 1 n 1 - 2 n 1 - submatrix augment det - : Δ n 1 - el A 1 n 1 - 1 n 2 - submatrix b augment det - : Δ n el A det - : k 2 : k n 2 - ≤ k k 1 + : Δ k el A 1 n 1 - 1 k 1 - submatrix b A 1 n 1 - k 1 + n 1 - submatrix augment det - : 11 line for 41 line Δ 1 el b : Δ n el A det : 21 line if dS Δ 1 el (2^Δ 2 el (2^+ Δ 3 el (2^+ sqrt : i 1 : i n ≤ i i 1 + : D2 i el Δ i el dS / : for 11 line t x D1 D2 : X0 tmin tmax N t x D1 rkfixed 131 line :

tmin 0 :

tmax 6 :

Δt 0.1 :

N tmax Δt / :

N 60

Матрица координат B X0 tmin tmax N D :

B1 B 1 N 24 submatrix :

B11 B1 γ * eval :

E1 E 1 E rows 24 submatrix :

E11 E1 γ * eval :

C u v F 2 sqrt v cos * 2 sqrt v sin * u 31 mat : u v F 1.5 - 2.5 02 π * 2020 CreateMesh 21 line :

C2 u v F u 2 sqrt v cos * 2 sqrt v sin * 31 mat : u v F 3 - 302 π * 1520 CreateMesh 21 line :

Trajectory1 B11 1 col B11 3 col augment eval :

Trajectory2 E11 1 col E11 3 col augment eval :

C2 γ * 1 col C2 γ * 3 col augment C γ * 1 col C γ * 3 col augment Trajectory2 Trajectory1 4 1 sys

τ 0 N 1 - range :

Trajectory2 Trajectory1 t 0 ≡ path B11 1 col t 1 + el B11 3 col t 1 + el . 15 red 1 5 mat : path path B11 1 col t 1 + el B11 3 col t 1 + el . 10 red 1 5 mat stack : if B11 1 col 1 el 0.8 - B11 3 col 1 el 0.2 + -1,1,-1 10 1 5 mat E11 1 col 1 el 0.5 - E11 3 col 1 el 1,1,1 10 1 5 mat XYZ TextXYZ 2.4 - 2 x2^2+x3^2=2 10 1 5 mat 2.4 - 1.8 x1^2+x2^2=2 10 1 5 mat 9 1 sys