I cannot use the roots function and must use one old good simple method by Newton.
What can you say about it?

Newton-Eng.docx (1,07 MiB) descargado 989 vez(es).

Very nice.

I had problems when reproducing your 2D version. Depending on initial guess I got errors "requested matrix element does not exist". The situation seemed to improve a little by pre-allocating the vectors X and Y yet it didn't become really robust.

So instead of assigning to j+1 elements I used augment() to add new points to the iteration history.

This is the original version:


Newton2D.sm (21,4 KiB) descargado 1559 vez(es).

And this is the modified version.


Newton2Da.sm (19,03 KiB) descargado 1680 vez(es).

Wrote

... can't find 2nd root

2nd root confirmed exact.

Solve(f,q,x0,x1,TOL) Copy.sm (162,62 KiB) descargado 1679 vez(es).

Tried to solve your functions with al_nleqsolve().

Regards

Newton.sm (13,85 KiB) descargado 1597 vez(es).

Another way, keeping the mathematical formulation and the iteration algorithm separate. As here.

Newton2Da.sm (13,34 KiB) descargado 1725 vez(es).



Best regards.
Alvaro.

Thanks Alvaro, added to my next Samples visit.
Take care ... Jean

Wrote

I cannot use the roots function

Actually, you can.

Regards

Wrote

I cannot use the roots function and must use one old good simple method by Newton.
What can you say about it?

An old vintage in Smath ... two examples.
Take care Valery ... Jean.

Solve_roost f1(x,y). f2(x,y).sm (27,42 KiB) descargado 1672 vez(es).

Two variations added.

Newton2Da. v2.sm (38,35 KiB) descargado 1669 vez(es).

Newton2Da. v2.pdf (152,11 KiB) descargado 1143 vez(es).

Best regards.
Alvaro.

Wrote

Two variations added.

Newton2Da. v2.sm (38,35 KiB) descargado 1669 vez(es).

Newton2Da. v2.pdf (152,11 KiB) descargado 1143 vez(es).

Best regards.
Alvaro.

From my future article
Если уравнений не два, а три, то матрица Якоби будет состоять из трёх строк и трёх столбцов. При этом графически отображать решение придётся не двумя кривыми, а тремя поверхностями на 3D-графике, причём их взаимное пересечение и будет отмечать решение (задание читателю).
Слабо!

You may find this Treasury useful/productive.


Solve(f,q,x0,x1,TOL) Copy.sm (300,16 KiB) descargado 1703 vez(es).