vb 获取两圆交点

Function Pow(a As Double, b As Double) As Double
    Pow = Format((a ^ b), "0.0000000000")
End Function
Sub circle_circle(ByVal xx0 As Double, ByVal yy0 As Double, ByVal r0 As Double, ByVal xx1 As Double, ByVal yy1 As Double, ByVal r1 As Double, ByRef X_1 As Double, ByRef Y_1 As Double, ByRef x_2 As Double, ByRef y_2 As Double)
Dim k As Double, k1 As Double
Dim intL As Double
Dim l0 As Double, l1 As Double, lh As Double
Dim center_x As Double, center_y As Double
Dim angl As Double
Dim x0 As Double, y0 As Double, x1 As Double, y1 As Double
x0 = CDbl(xx0)
y0 = CDbl(yy0)
x1 = CDbl(xx1)
y1 = CDbl(yy1)
intL = Sqr(Pow(y1 - y0, 2) + Pow(x1 - x0, 2))
l1 = (Pow(r0, 2) - Pow(r1, 2) - Pow(intL, 2)) / (-2 * intL)
l0 = intL - l1
Dim p1 As Double, p2 As Double, p3 As Double
p1 = Pow(r1, 2)
p2 = Pow(l1, 2)
p3 = p1 - p2
lh = Sqr(p3)
If (xx1 = xx0) Then '圆心垂直，交点水平对称
angl = 90
If (y1 > y0) Then
center_x = x0
center_y = y0 + l0
x_2 = CDbl(center_x - lh)
y_2 = CDbl(center_y)
X_1 = CDbl(center_x + lh)