On 9 bře, 14:57, "[email protected]" <[email protected]> wrote:
> I got the following message from Andrej Vodopivec:
>
> You can solve this system in maxima with to_poly_solve:
>
> (%i2) to_poly_solve([(x-5)^2+y^2-16, (y-3)^2+x^2-9], [x,y],
> use_grobner=true);
> (%o2) %union([x=-(9*sqrt(55)-135)/68,y=-(3*5^(3/2)*sqrt(11)-123)/68],
> [x=(9*sqrt(55)+135)/68,y=(3*5^(3/2)*sqrt(11)+123)/68])
>
> HTH, Andrej
And this is Andrej's solution in Sage
Andrej, many thanks for pointing to use_grobner.
R.
----------------------------------------------------------------------
| Sage Version 4.3.3, Release Date: 2010-02-21 |
| Type notebook() for the GUI, and license() for information. |
----------------------------------------------------------------------
sage: x,y=var('x y')
sage: c1(x,y)=(x-5)^2+y^2-16; c2(x,y)=(y-3)^2+x^2-9
sage: M=c1._maxima_().parent()
sage: M.to_poly_solve([c1,c2],[x,y],'use_grobner=true').sage()
[[x == -9/68*sqrt(55) + 135/68, y == -15/68*sqrt(5)*sqrt(11) +
123/68], [x == 9/68*sqrt(55) + 135/68, y == 15/68*sqrt(5)*sqrt(11) +
123/68]]
sage:
--
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org
