[Maxima] solution of system of quadratic eqns.
fateman at eecs.berkeley.edu
Sat Jun 23 16:09:47 CDT 2012
On 6/23/2012 9:55 AM, Barton Willis wrote:
> Surely the bug is related to the fact that the last call to resultant returns 256*(sqrt(7)^2-7) instead of 0.
Ah, yes that would account for it.
I just try typing
and got 0.
So the result in CRE ("rat") form is simplified according to its lights.
Here's a fix:
:lisp (setf (symbol-function 'oldres)(symbol-function '$resultant))
:lisp (defun $resultant(a b c)($ratsimp(oldres a b c)))
now solve works.
More information about the Maxima