# -*- coding: utf-8 -*-import mathimport numpy as npdef insec(p1,r1,p2,r2): x = p1[0] y = p1[1] R = r1 a = p2[0] b = p2[1] S = r2 d = math.sqrt((abs(a-x))**2 + (abs(b-y))**2) if d > (R+S) or d < (abs(R-S)): print ("Two circles have no intersection") return elif d == 0 and R==S : print ("Two circles have same center!") return else: A = (R**2 - S**2 + d**2) / (2 * d) h = math.sqrt(R**2 - A**2) x2 = x + A * (a-x)/d y2 = y + A * (b-y)/d x3 = round(x2 - h * (b - y) / d,2) y3 = round(y2 + h * (a - x) / d,2) x4 = round(x2 + h * (b - y) / d,2) y4 = round(y2 - h * (a - x) / d,2) print (x3, y3) print (x4, y4) c1=np.array([x3, y3]) c2=np.array([x4, y4]) return c1,c2 P1=np.array([-5,0])R1=10P2=np.array([5,0])R2=5C=insec(P1,R1,P2,R2)C1=C[0]C2=C[1] 
