<p>Draghilev's method</p>

8541ce43-7ebb-4fab-b57e-08818739186d 17 4 5 false false 0 decimal

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

Draghilev's method. Examples. Viacheslav N. Mezentsev, viacheslavmezentsev@ya.ru

#n #v NameVec #v 1 #n range el vectorize 11 line :

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

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 :

1 F X.0 t.min t.max N Draghilev out 0 : n X.0 F eval rows : 12 mat X 2 n * 1 +_u NameVec : X1 X2 12 mat X 1 n range el X n 2 + (2 n * 1 + (range el 12 mat : S X1 F_u n 1 + el X2 F * - evalm : k 1 n 1 + (range Y X1 : k n 1 + (< Y k el_u n 1 + el : if jac S vector convert Y vector convert jacobian simplify det evalf maple : k n 1 + ≡ out k el jac : out k el jac - : if 41 line for k n 2 + (2 n * 1 + (range out k el 0 : for Xo X.0 1 X.0 stack eval : result Xo t.min t.max N t_u D out : rkfixed : V result n 2 + col : out 0 : 12 mat Roots V Search : k 1 Roots length range m 2 n 1 + (range out k m 1 - el result m col Roots k el el result m col Roots k el 1 + el V Roots k el el V Roots k el 1 + el Interpol eval : for for out 121 line :

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

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

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

4 appVersion 0.99.6858.3232

X.0 01 stack :

t.min 0 :

t.max 1.6 π * :

N 100 :

result X F X.0 t.min t.max N Draghilev :

k result rows :

k 12

m 1 k range o_O m el o : for

m 1 k range o_# m el m num2str : for

XY result 1 k range 12 range el :

x y f1 x y f2 XY o_O augment XY o_# augment 4 1 sys

XY 0.4552 1.1463 0.6541 1.1508 1.1002 0.8526 1.0989 0.244 0.8316 0.2772 - 0.1656 0.9063 - 0.4547 - 1.1462 - 0.6531 - 1.1508 - 1.1002 - 0.8529 - 1.099 - 0.2439 - 0.8312 - 0.2782 0.1663 - 0.9061 122 mat