<p>Draghilev's method</p>

4302de35-c0ae-4b63-be11-9bdedd5e64cd 12 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 z f1 3 x * y z * cos - 12 / - :

x y z f2 x 2^81 y 0.1 + (2^* - z sin + 1.06 + :

x y z f3 20 z * e x y * -^+ 13 / 3 - 10 π * + (* + :

X F x X 1 el : y X 2 el : z X 3 el : 13 mat x y z f1 x y z f2 x y z f3 31 mat 21 line :

4 appVersion 0.99.6858.3232

X.0 0.489 0.5 0.513 - stack :

t.min 0 :

t.max 0.0008 :

N 50 :

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

k result rows :

k 2

XY result 1 k range 12 range el :

m 1 k range o_O m el o : for

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

XYZ result 1 k range 13 range el :

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

XYZ 0.5 0.0009 - 0.5236 - 0.4982 0.1988 - 0.5288 - 23 mat