草草计算了一下,不知道有没算错,你试试吧
var
d, a,b,t1,t2, x0, y0, x,y: Double
d := sin(c/2)
t1 := ((sqr(1/a))-1)*(sqr(x2-x1)+sqr(y2-y1))
t2 := sqe(1/a)*(sqr(x2-x1)+sqr(y2-y1))
a := (t1-t2)/2
b := (sqr(x1)+sqr(x2)+sqr(y1)+sqr(y2)-(t1+t2))/2
x0 := (((b+a)*y2)-((b-1)*y1))/(2*(x1*y2-x2*y1))
y0 := (((b+a)*x2)-((b-a)*x1))/(2*(x2*y1-x1*y2))
x :=(x0+x1)/2
y := (y0+y1)/2
(x,y)就是圆心坐标