<p>Draghilev's method</p>

88b47f81-0219-4cf9-8e76-64af627fd560 20 4 5 decimal

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

Draghilev's method. Examples. Viacheslav N. Mezentsev, viacheslavmezentsev@gmail.com

x y f1 5 x 2^* sin 4 y 2^* sin + :

x y f2 x 2^y x * - y 2^+ 1 - :

X F x X 1 el : y X 2 el : 12 mat x y f1 x y f2 21 mat 21 line :

<p>Draghilev's method</p>

n_v NameVec V 0 : ii 1 n range V ii el_v ii el : for V 31 line :

v Search out 0 : n 1 : k 1 v length 1 - range v k el 0 > (v k 1 + el 0 < (& v k el 0 < (v k 1 + el 0 > (& | v k el 0 ≡ (| out n el k : n n 1 + : 21 line if 11 line for out 41 line :

p1 p2 v1 v2 Interpol p1 v1 v2 v1 - / p2 p1 - (* - :

1 F X0 tmin tmax N Draghilev out 0 : n X0 F eval rows : 12 mat X 2 n * 1 +_u NameVec : X1 X 1 n 11 submatrix : X2 X n 2 + 2 n * 1 + 11 submatrix : S X1 F_u n 1 + el X2 F * - evalm : ii 12 n * 1 + range ii n 1 + (> out ii el 0 : continue 21 line if Y X1 : ii n 1 + (< Y ii el_u n 1 + el : if jac S vector convert Y vector convert jacobian simplify det evalf maple : ii n 1 + ≡ out ii el jac : out ii el jac - : if 51 line for t_u D out : X0 X0 1 X0 stack eval : 12 mat result X0 tmin tmax N t_u D rkfixed : V result n 2 + col : out 0 : 12 mat Roots V Search : ii 1 Roots length range jj 2 n 1 + range out ii jj 1 - el result jj col Roots ii el el result jj col Roots ii el 1 + el V Roots ii el el V Roots ii el 1 + el Interpol eval : for for out 121 line :

X0 01 stack :

tmin 0 :

tmax 4.5 :

N 40 :

result X F X0 tmin tmax N Draghilev :

k result rows :

k 12

Points result 1 col result 2 col augment :

eq f1 f2 Points 31 sys :