c89eae30-78dc-428a-9a0d-a2d4f23b9fe0 65 4 5 false false 0 decimal

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

The intersection points of two curves

x0,y0-initial approximation,n-number of solving x0 y0 n inter sol StepMax 0 : ε 1014 -^: X0 1 el x0 : X0 2 el y0 71 Random * 2 - (1 el + : 14 mat u 1 el x 1 el : u 2 el x 2 el : x Jac u f u Jacobian : 13 mat X0 StepMax ε x f x Jac al_nleqsolve 31 line : k 1 n range a k el sol transpose : δ k el a k el 1 col 1 el a k el 2 col 1 el f1 abs : (accuracy of solution ≡ k 1 ≡ r a k el δ k el augment : r r a k el δ k el augment stack : if 31 line for rc r 3 csort : (sorting by δ ≡ k 1 n range rc 3 col k el 102 -^< root k el rc 1 n 12 submatrix k row : break if k 1 ≡ res root k el . 6 red 13 mat augment : res res root k el . 6 red 13 mat augment stack : if 21 line for res 51 line :

x f x 2 el abs cos - 2 x 2 el * x 2 el 2^abs cos * 2 x 1 el * abs cos * + x 1 el abs sin - 2 x 2 el 2^abs sin * 2 x 1 el * abs sin * + 21 mat :

x y f1 y abs cos - 2 y * y 2^abs cos * 2 x * abs cos * + :

x y f2 x abs sin - 2 y 2^abs sin * 2 x * abs sin * + :

A1 3.5 - 1.1 - 40 inter :

A2 1 - 1 - 40 inter :

A3 1.3 1 - 30 inter :

A5 31 - 40 inter :

A4 21 - 40 inter :

f1 f2 A1 A2 A3 A4 A5 7 1 sys