<p>Draghilev's method</p>

5375561a-7fb3-4b0b-b545-d02421478b65 29 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 :

4 appVersion 0.99.6858.3232

X F x X 1 el : y X 2 el : 12 mat x 2^y 2^+ 8 - x y * sin x y * exp sin * 21 mat 21 line :

X.0 2.784 0.5 - stack :

t.min 0 :

t.max 13 :

N 500 :

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

k result rows :

k 80

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 * sin x y * exp sin * x 2 ^ y 2 ^ + 8 - XY o_O augment XY o_# augment 4 1 sys