Calculating the Center of A Circle



-----maxima-bounces at math.utexas.edu wrote: -----

>Hello,
>    s: solve([(x-x0)^2+(y-y0)^2=r^2,(x-x1)^2+(y-y1)^2=r^2],[x,y]);
>This is a relatively simple problem that my Mathcad can solve in less
>than a second.
>Maxima takes about a minute to come up with a [ ] answer.
>The others on sci.math.symbolic agree this is a bug or should be.

Yes, this is a bug. If you could, report it to our bug list. The
function algsys isn't nearly as good as it should be. Also:

  (%i10) eq : [y0^2-2*y*y0+y^2+x0^2-2*x*x0+x^2=r^2,y1^2-2*y*y1+y^2
  +x1^2-2*x*x1+x^2=r^2]$

OK:

  (%i11) algsys(eq,[x,y,r])$

Not OK:

  (%i12) algsys(eq,[r,x,y]);
  Maxima encountered a Lisp error:

Set algebraic to true and try again:

  (%i13) algebraic : true$

No Lisp error, but no solution either :(

 (%i14) algsys(eq,[r,x,y]);
 (%o14) []

Barton