The same results also happen in Sage 4.7 on 64 bit Debian. var ('P X Y') solve ([(1/2) == -1/2*(2*P - 1)/Y + P/X, Y == 2*P*X/(2*P + 1), Y == P + 1/2], P,X,Y)
On Jul 29, 10:52 pm, maor <maor....@gmail.com> wrote: > I tried to solve the following simple 3 equations with 3 variables: > > sage: var ('P X Y') > (P, X, Y) > > sage: eq > [(1/2) == -1/2*(2*P - 1)/Y + P/X, Y == 2*P*X/(2*P + 1), Y == P + 1/2] > > sage: solve(eq,[P,X,Y]) > [[P == (-1/2), X == 0, Y == 0]] > > I know that there are 2 solutions, but solve returned a bad solution > (X==0??? it is in the denominator) > > while playing with it a bit I discovered that the following (strange) > definition works: > > sage: eq2 > [Y == P + 1/2, [(1/2) == -1/2*(2*P - 1)/Y + P/X, Y == 2*P*X/(2*P + > 1)]] > > sage: solve(eq2,[P,X,Y]) > [[P == -1/2*sqrt(2) - 1/2, X == -sqrt(2) + 1, Y == -1/2*sqrt(2)], [P > == 1/2*sqrt(2) - 1/2, X == sqrt(2) + 1, Y == 1/2*sqrt(2)]] > > I am working with sage 4.5.2 on Ubuntu 10.04 64 bit > > Thanks, > Maor -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org