Algsys Not Solving A Simultaneous Polynomial I Would Expect It To
Subject: Algsys Not Solving A Simultaneous Polynomial I Would Expect It To
From: Tawny Owl
Date: Fri, 03 Nov 2006 00:12:09 +0000
Hi Everyone,
I am a big fan of Maxima, but I have stumbled across a simultaneous
polynomial which algsys (and hence solve) doesn't seem to be able to
resolve, and I'm wondering if anyone can tell me what the problem is,
please. Here's the original command:
(%i7) assume (s>0)$
(%i8) algsys([x^2+y^2=4,(s-x)^2+y^2=16,x^2+(s-y)^2=25],[s,x,y]);
(%o8) []
I then tried rewriting the 3 polynomials with their factors expanded:
(%i9) e1: y^2+x^2-4$
(%i10) e2: y^2+x^2-2*s*x+s^2-16$
(%i11) e3: y^2-2*s*y+x^2+s^2-25$
(%i12) algsys([e1, e2, e3], [s,x,y]);
(%o12) []
This would imply that there are no solutions. However - this is not true -
these simultaneous polynomials can be solved. If I set values for s, x and
y, and then calculate the expressions, I can obtain the correct results:
(%i57) s: sqrt(sqrt(511)/2+41/2)$
(%i58) x: sqrt(33*sqrt(511)/520+857/520)$
(%i59) y: sqrt(1223/520-33*sqrt(511)/520)$
(%i61) x^2 + y^2;
(%o61) 4
(%i63) y^2+x^2-2*s*x+s^2,numer;
(%o63) 16.0
(%i70) y^2-2*s*y+x^2+s^2,numer;
(%o70) 25.0
I'm not sure if this is related to the algsys problem, but I noticed just
now that to calculate those expressions, line %i61 above was fine, wheras
for each of the lines %i63 and %i70, I had to append "numer" to get Maxima
to calculate the value.
Thanks for any advice!
Get an INSTANT game of Mastermind at
http://mysite.wanadoo-members.co.uk/grahamlaight/jscript/index.html
_________________________________________________________________
Windows Live? Messenger has arrived. Click here to download it for free!
http://imagine-msn.com/messenger/launch80/?locale=en-gb