<p> Draghilev Jacobian/rkfixed Solver </p>

5dedc6aa-5d9b-4ac2-8aba-e29a846bdb93 2 6 5 decimal

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

X0 tmin tmax N D n f rows 1 + : k 1 n range u k el x k el : for b 1 - (f u Jacobian 1 n 1 - n n submatrix * : матрица свободных членов ≡ A f u Jacobian 1 n 1 - 1 n 1 - submatrix : основная(квадратная) матрица ≡ n 2 > det 1 el b A 1 n 1 - 2 n 1 - submatrix augment det - : определителu дополнит.матриц ≡ det n 1 - el A 1 n 1 - 1 n 2 - submatrix b augment det - : det n el A det - : k 2 : k n 2 - ≤ k k 1 + : det k el A 1 n 1 - 1 k 1 - submatrix b A 1 n 1 - k 1 + n 1 - submatrix augment det - : for 41 line det 1 el b : det n el A det : 21 line if dS 1 : i 1 n range D2 i el det i el dS / : for t x D1 D2 : Матрица  правых частей  системы ОДУ ≡ X0 tmin tmax N t x D1 rkfixed 91 line :

x y f x 4^2 x 2^* y * + 3 x * y 2^* - y 3^+ :

1.66235897862237 0.562279512062302 i * -

function for Draghilev. f 1 el x 1 el 4^2 x 1 el 2^* x 2 el * + 3 x 1 el * x 2 el 2^* - x 2 el 3^+ :

init1 1.66235897862237 0.562279512062302 i * - 21 mat :

Contours are defined by β ... solutions are hyper-reflexive.
'N' defines the fineness ... be generous.
Each contour closes on ≠ length of respective solutions
The best closed contours depend upon 't1max, t2max'. 
Adjust these two manually from the dual observer. tmin 0 : t1max 0.33 : N 100 : 31 sys

B1 init1 tmin t1max N D : 11 sys

The parametrization Human eye by Draghilev method

B1 2 col B1 3 col augment 1 1 sys

x y f x 4^2 x 2^* y * + 3 x * y 2^* - y 3^+ :

R 1 y f 0 ≡ y solve maple : 0.324717957244745 - 1.66235897862237 0.562279512062302 i * - 1.66235897862237 0.562279512062302 i * + 31 mat