def estimation(rectangle, point1, point2, rayon1, rayon2, num_p=100000): p1 = point1 p2 = point2 c1 = Circle(p1, rayon1) c2 = Circle(p2, rayon2) frame = rectangle.framing_rectangle(c1, c2) cpt = 0 for i in range(num_p + 1): p = rectangle.random_point() # print(f"Does Point: {p} belongs to Circle :{c1} & Circle :{c2}") #print((c1.belongs(p) and c2.belongs(p))) if (c1.belongs(p) & c2.belongs(p)): cpt += 1 return (cpt / num_p) * r.surface()