円の交点を求めるプログラムが出来ない。

ある座標系に、中心A(a,b) 半径r1の円アと中心B(c,d) 半径r2の円イがある。
円、アとイの交差するポイントを求めよ（交差しない時も判断せよ）

という問題が、Pythonでどう書いたらいいものかまったく検討がつかない。

円の式は(上記円アの場合)

(x-a)^2 + (y-b)^2 = r^2

となるので、コレを０＝の式にして、円ア ＝ 円イの式にして
といていけばいいのだが、手と紙を使えばなんとも無い問題だが
プログラミングでどうやって書けばいいのかがわからない。

引き続き考える。