Rank: Advanced Member Groups: Registered, Advanced Member Joined: 06/04/2023(UTC) Posts: 348 Was thanked: 28 time(s) in 27 post(s)
|
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,106kb) downloaded 29 time(s).Edited by user 08 April 2023 22:12:41(UTC)
| Reason: Not specified
|
1 user thanked OchkovVF for this useful post.
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,990 Was thanked: 1128 time(s) in 724 post(s)
|
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 (22kb) downloaded 10 time(s).And this is the modified version. Newton2Da.sm (20kb) downloaded 14 time(s). |
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866 Was thanked: 981 time(s) in 809 post(s)
|
Originally Posted by: Jean Giraud ... can't find 2nd root 2nd root confirmed exact. Solve(f,q,x0,x1,TOL) Copy.sm (163kb) downloaded 7 time(s).
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 23/07/2013(UTC) Posts: 1,129 Was thanked: 510 time(s) in 340 post(s)
|
Tried to solve your functions with al_nleqsolve(). Regards Newton.sm (14kb) downloaded 12 time(s).
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 28/08/2014(UTC) Posts: 1,368 Was thanked: 829 time(s) in 525 post(s)
|
Another way, keeping the mathematical formulation and the iteration algorithm separate. As here. Newton2Da.sm (14kb) downloaded 11 time(s).Best regards. Alvaro. Edited by user 10 April 2023 00:27:45(UTC)
| Reason: Not specified
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866 Was thanked: 981 time(s) in 809 post(s)
|
Thanks Alvaro, added to my next Samples visit. Take care ... Jean
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 23/07/2013(UTC) Posts: 1,129 Was thanked: 510 time(s) in 340 post(s)
|
Originally Posted by: OchkovVF I cannot use the roots function Actually, you can. Regards
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866 Was thanked: 981 time(s) in 809 post(s)
|
Originally Posted by: OchkovVF 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 (28kb) downloaded 10 time(s).
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 28/08/2014(UTC) Posts: 1,368 Was thanked: 829 time(s) in 525 post(s)
|
|
1 user thanked Razonar for this useful post.
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 06/04/2023(UTC) Posts: 348 Was thanked: 28 time(s) in 27 post(s)
|
Originally Posted by: Razonar From my future article Если уравнений не два, а три, то матрица Якоби будет состоять из трёх строк и трёх столбцов. При этом графически отображать решение придётся не двумя кривыми, а тремя поверхностями на 3D-графике, причём их взаимное пересечение и будет отмечать решение (задание читателю). Слабо!
|
1 user thanked OchkovVF for this useful post.
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866 Was thanked: 981 time(s) in 809 post(s)
|
|
|
|
|
Forum Jump
You cannot post new topics in this forum.
You cannot reply to topics in this forum.
You cannot delete your posts in this forum.
You cannot edit your posts in this forum.
You cannot create polls in this forum.
You cannot vote in polls in this forum.