Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,494 Was thanked: 1274 time(s) in 745 post(s)
|
Very strange behavior. |
Russia ☭ forever Viacheslav N. Mezentsev |
1 user thanked uni for this useful post.
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 28/08/2014(UTC) Posts: 1,356 Was thanked: 815 time(s) in 516 post(s)
|
Hi Viacheslav. Numerical procedures: You can always find examples where they fail. I find this variation of your example particularly worrisome since the function changes sign in the interval where the zeros are requested. It seems that it takes as a tolerance something proportional to the width of the interval It is true that apparently the SMath algorithm can be improved in several ways. However, something good must have. Here SMath succeeds where Maple, Maxima and FindRoot do not. solve_bug.sm (10kb) downloaded 7 time(s).Best regards. Alvaro.
|
1 user thanked Razonar for this useful post.
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 23/07/2013(UTC) Posts: 1,128 Was thanked: 508 time(s) in 338 post(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: uni Very strange behavior. Often, strange behavior results from numerical instability. I have observed that fact between MCD 11 XP 32 bits vs same SS 6179 32 bits. The best example is the Mathematica 4.0 famous function. Solving the Goat radius is another example ... no way for 15 decimals ! Mathsoft Minerr uses a specific technique to cope with this fact, i.e: fool two similar results by injecting a small random. As well, Smath solvers turn unstable, pedal within their own end iterate. You forgot to show the function. Cheers ... Jean.
|
|
|
|
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.