mnewton divide by zero
- Subject: mnewton divide by zero
- From: Edwin Woollett
- Date: Thu, 3 Apr 2008 11:46:44 -0700
lisp error message with mnewton:
(%i18) eqns : [ x + y = 3, x^2 + y^2 = 9 ]$
(%i19) mn(x1, x2) := mnewton( eqns, [x, y], [x1, x2] )$
(%i20) solve( eqns, [x, y] );
(%o20) [[x = 3, y = 0], [x = 0, y = 3]]
(%i21) mn(1, 2);
(%o21) [ [x = 6.37714165E-17, y = 3.0] ]
......
(%i24) mn(0, 0);
Maxima encountered a Lisp error:
Error in FUNCALL [or a callee]: Zero divisor.
Automatically continuing.
To reenable the Lisp debugger set *debugger-hook* to nil.
(%i25) mn(0, 1);
(%o25) [ [x = 1.13978659E-16, y = 3.0] ]
Is this a bug?
Ted Woollett