そういえばそうだな…

円の方程式を
x^2+y^2+lx+my+n=0
という形で表して（円の方程式をばらしただけ、数Ⅱの教科書参照）
ばらした式に３つの点代入して３元連立１次方程式を解いてlとmとnを求めて
元の式(x-a)^2+.....と係数比較すれば強引に求まるかも

しかし検索したら解法のアルゴリズムとかありそうだなあ…