Hi Damiano,
By the way, the formulas in the script you attached to your
previous email are slightly different from your original post.
> i' ve found another problem with maxima , under some circumstances
> the mnewton script crashes . see the attachment .
This output --
(%i9) mnewton([F1(N,F),F2(N,F)],[N,F],[1.0,1.0]);
Maxima encountered a Lisp error:
Error in MACSYMA-TOP-LEVEL [or a callee]:
Error in CONDITIONS::CLCS-UNIVERSAL-ERROR-HANDLER [or a callee]:
Can't print a non-number.
Fast links are on: do (use-fast-links nil) for debugging
Broken at PCL::PRINT-STD-INSTANCE. Type :H for Help.
1 (Abort) Return to top level.
dbl:MAXIMA>>
appears to be GCL's reaction to a floating point overflow.
(Lisp in general, including GCL apparently, isn't
comfortable with the IEEE 754 special values
for not-a-number and infinity.)
The floating point overflow means only that mnewton diverged.
So there's really nothing to do but try another starting point.
I wonder why "Can't print a non-number" wasn't caught by
the usual error-trapping mechanism. Surely we want to
return to the Maxima input prompt in this case.
For what it's worth,
Robert Dodier
____________________________________________________
Yahoo! Sports
Rekindle the Rivalries. Sign up for Fantasy Football
http://football.fantasysports.yahoo.com