Определить координаты центра и минимальный радиус круга, который содержит все заданные треугольники

import time
import math
start = time.time()
                        
N = 7
maxdistance = 0
x = [3, -2,  1, -1, 8,  6, 3, 4, 7, -1, -3,  1,  3, -4, -6, -5, -8, 0, 0, -3, -3]
y = [1, -1, -4,  4, 1, -3, 7, 5, 6, -4, -6, -8, -7,  2, -2, -5,  1, 6, 7,  6,  9]
 
for i in range (0,N*3):
    for j in range (0,N*3):
        distance = math.sqrt((x[i] - x[j])**2 + (y[i] - y[j])**2)
        if distance > maxdistance:
            maxdistance = distance
            x1 = x[i]
            y1 = y[i]
            x2 = x[j]
            y2 = y[j]
x0 = (x1 + x2)/2
y0 = (y1 + y2)/2
R = maxdistance/2
print("Точка 1: [",x1," ; ",y1,"]")
print("Точка 2: [",x2," ; ",y2,"]")
print("Максимальная дистанция: ",R*2)
for i in range (0,N*3):
    distance = math.sqrt((x[i] - x0)**2 + (y[i] - y0)**2)
    if distance > R:
        x3 = x[i]
        y3 = y[i]
        R = distance
        x0 = (-((y1 - y2)*((x3**2) + (y3**2)) + (y2 - y3)*((x1**2) + (y1**2)) + (y3 - y1)*((x2**2) + (y2**2)))
                    /(2*((x1 - x2)*(y3 - y1) - (y1 - y2)*(x3 - x1))))
        y0 = (((x1 - x2)*((x3**2) + (y3**2)) + (x2 - x3)*((x1**2) + (y1**2)) + (x3 - x1)*((x2**2) + (y2**2)))
                    /(2*((x1 - x2)*(y3 - y1) - (y1 - y2)*(x3 - x1))))
        i = -1
 
#Проверка
for i in range (0,N*3):
    RA = math.sqrt((x[i] - x0)**2 + (y[i] - y0)**2)
    if RA > R:
        print("Ошибка!!!")
print("Минимальный радиус окружности равен: ",R)
print("Центр окружности: [",x0," ; ",y0,"]")
finish = time.time()
print("Время работы программы составило: ",finish - start,"c")