Для каждого треугольника триангуляции Делоне нарисовать описанную окружность

import matplotlib.pyplot as plt
import scipy.spatial as sp
import numpy as np
 
# Generate some random points
points = np.random.rand(10, 2)
 
# Create a Delaunay triangulation object
tri = sp.Delaunay(points)
 
# Plot the triangulation
plt.triplot(points[:,0], points[:,1], tri.simplices)
 
# Get the circumcenters of each triangle
circumcenters = tri.circumcenters
 
# Plot the circumcircles
for c in circumcenters:
    # Calculate the radius of the circle
    r = np.sqrt(np.sum((points[tri.simplices[0]] - c)**2, axis=1))[0]
    # Draw the circle
    circle = plt.Circle(c, r, fill=False)
    plt.gca().add_patch(circle)
 
plt.show()

import matplotlib.pyplot as plt
import scipy.spatial as sp
import numpy as np
 
# Generate some random points
points = np.random.rand(10, 2)
 
# Create a Delaunay triangulation object
tri = sp.Delaunay(points)
 
# Plot the triangulation
plt.triplot(points[:,0], points[:,1], tri.simplices)
 
# Calculate the circumcenters manually
circumcenters = []
for simplex in tri.simplices:
    # Get the coordinates of the three vertices of the simplex
    vertices = points[simplex]
    # Calculate the midpoint and perpendicular bisector of each side
    midpoints = [(vertices[i] + vertices[j])/2 for i, j in [(0, 1), (1, 2), (2, 0)]]
    bisectors = [-1*np.array([mp[1]-mp[0], mp[0]-mp[1]]) for mp in midpoints]
    # Calculate the intersection of two bisectors to get the circumcenter
    intersections = []
    for i, j in [(0, 1), (1, 2), (2, 0)]:
        A, B = bisectors[i], bisectors[j]
        b = np.dot(B,B)
        A_dot_B = np.dot(A,B)
        A_dot_A = np.dot(A,A)
        D = A_dot_A*b - A_dot_B**2
        if D != 0:
            C = [(A_dot_A*np.dot(A,B) - A_dot_B*np.dot(A,A))/D, (A_dot_A*np.dot(B,A) - A_dot_B*np.dot(B,B))/D]
            intersections.append(C)
    if intersections:
        circumcenters.append(intersections[0])
 
# Plot the circumcircles
for c in circumcenters:
    # Calculate the radius of the circle
    r = np.sqrt(np.sum((points[tri.simplices[0]] - c)**2, axis=1))[0]
    # Draw the circle
    circle = plt.Circle(c, r, fill=False)
    plt.gca().add_patch(circle)
 
plt.show()