mathmatic求两圆交点

求两圆的交点，Solve[(x - a1)*(b2 - b1) == (y - b1)*(a2 - a1) && (x - c1)*(d2 -
                                           d1) == (y - d1)*(c2 - c1), {x, y}]，当然这个很恐怖

即使化简也很恐怖{{x -> (1/(
   2 ((a1 - a2)^2 + (b1 - b2)^2)))(a1^3 - a1^2 a2 + a2^3 +
     a1 (-a2^2 + (b1 - b2)^2 - r1^2 + r2^2) -
     Sqrt[-(b1 - b2)^2 (((a1 - a2)^2 + (b1 - b2)^2 - r1^2)^2 -
        2 ((a1 - a2)^2 + (b1 - b2)^2 + r1^2) r2^2 + r2^4)] +
     a2 ((b1 - b2)^2 + (r1 - r2) (r1 + r2))),
  y -> (1/(2 ((a1 - a2)^2 + (b1 - b2)^2) (b1 -
      b2)))(a1^2 (b1 - b2) (b1 + b2) +
     a2^2 (b1 - b2) (b1 + b2) + (b1 - b2)^2 (b1^2 - b2^2 - r1^2 +
        r2^2) - a2 Sqrt[-(b1 -
          b2)^2 (((a1 - a2)^2 + (b1 - b2)^2 - r1^2)^2 -
         2 ((a1 - a2)^2 + (b1 - b2)^2 + r1^2) r2^2 + r2^4)] +
     a1 (2 a2 (-b1^2 + b2^2) +
        Sqrt[-(b1 - b2)^2 (((a1 - a2)^2 + (b1 - b2)^2 - r1^2)^2 -
           2 ((a1 - a2)^2 + (b1 - b2)^2 + r1^2) r2^2 +
           r2^4)]))}, {x -> (1/(
   2 ((a1 - a2)^2 + (b1 - b2)^2)))(a1^3 - a1^2 a2 + a2^3 +
     a1 (-a2^2 + (b1 - b2)^2 - r1^2 + r2^2) +
     Sqrt[-(b1 - b2)^2 (((a1 - a2)^2 + (b1 - b2)^2 - r1^2)^2 -
        2 ((a1 - a2)^2 + (b1 - b2)^2 + r1^2) r2^2 + r2^4)] +
     a2 ((b1 - b2)^2 + (r1 - r2) (r1 + r2))),
  y -> (1/(
   2 ((a1 - a2)^2 + (b1 - b2)^2) (b1 -
      b2)))(a1^2 (b1 - b2) (b1 + b2) +
     a2^2 (b1 - b2) (b1 + b2) + (b1 - b2)^2 (b1^2 - b2^2 - r1^2 +
        r2^2) + a2 Sqrt[-(b1 -
          b2)^2 (((a1 - a2)^2 + (b1 - b2)^2 - r1^2)^2 -
         2 ((a1 - a2)^2 + (b1 - b2)^2 + r1^2) r2^2 + r2^4)] -
     a1 (2 a2 (b1 - b2) (b1 + b2) +
        Sqrt[-(b1 - b2)^2 (((a1 - a2)^2 + (b1 - b2)^2 - r1^2)^2 -
           2 ((a1 - a2)^2 + (b1 - b2)^2 + r1^2) r2^2 + r2^4)]))}}，

束手无策，

但是人类智慧是无穷的，把圆圆心相对于{a1,b1}偏移就简单多了

Solve[(x)^2 + (y)^2 == r1^2 && (x - a)^2 + (y - b)^2 == r2^2, {x,
  y}]   a=a2-a1,b=b2-b1

{{x -> (a^3 + a (b^2 + r1^2 - r2^2) - Sqrt[
    b^2 (-(a^2 + b^2 - r1^2)^2 + 2 (a^2 + b^2 + r1^2) r2^2 - r2^4)])/(
   2 (a^2 + b^2)),
  y -> (b^4 + b^2 (a^2 + r1^2 - r2^2) +
    a Sqrt[b^2 (-(a^2 + b^2 - r1^2)^2 + 2 (a^2 + b^2 + r1^2) r2^2 -
        r2^4)])/(2 b (a^2 + b^2))}, {x -> (
   a^3 + a (b^2 + r1^2 - r2^2) + Sqrt[
    b^2 (-(a^2 + b^2 - r1^2)^2 + 2 (a^2 + b^2 + r1^2) r2^2 - r2^4)])/(
   2 (a^2 + b^2)),
  y -> (b^4 + b^2 (a^2 + r1^2 - r2^2) -
    a Sqrt[b^2 (-(a^2 + b^2 - r1^2)^2 + 2 (a^2 + b^2 + r1^2) r2^2 -
        r2^4)])/(2 b (a^2 + b^2))}}

化简  c=a^2+b^2

  x (a(c+ r1^2 - r2^2)+-b((r1 + r2)^2-c))/(2c)
  y (b(c+ r1^2 - r2^2)-+a((r1 + r2)^2-c))/(2c)

x=x+a1

y=y+b1