Найти окружность проходящую через три точки на плоскости

data = [(1, 1), (-1, 3), (5, 3), (-3, -4), (6, 0), (1, 5), (-2, 0), (3, 3), (0, 0)]
minimum = {}
points = [[data[i], data[j], data[k]] for i in range(len(data)-2) for j in range(i+1, len(data)-1) for k in range(j+1, len(data))]
for z in points:
    x1, y1 = z[0]
    x2, y2 = z[1]
    x3, y3 = z[2]
    try:
        y0 = (2*(x1-x3)*(y2**2 - y1**2 - x1**2 + x2**2) + 2*(x1-x2)*(x1**2 - x3**2 - y3**2 + y1**2)) / (4*(y2-y1)*(x1-x3) - 4*(x1-x2)*(y3-y1))
    except ZeroDivisionError:
        continue
    if x1 - x3 != 0:
        x0 = ((x1**2 - x3**2 - y3**2 + y1**2) + 2*y0*(y3-y1))/ (2*(x1-x3))
    if x1 - x2 != 0:
        x0 = ((x1**2 - x2**2 - y2**2 + y1**2) + 2*y0*(y2-y1))/ (2*(x1-x2))
    R = ((x0 - x1)** 2 + (y0 - y1)** 2) ** (1/2)
 
    n1, n2, n3 = 0, 0, 0
 
    for t in data:
        if (t[0] - x0) ** 2 + (t[1] - y0) ** 2 < R ** 2:
            n1 += 1
        elif (t[0] - x0) ** 2 + (t[1] - y0) ** 2 == R ** 2:
            n3 += 1
        else:
            n2 += 1
        minimum[abs(n1 - n2)] = ['(x - {:.1f}) ** 2 + (y - {:.1f}) ** 2 = {:.1f} ** 2'.format(x0, y0, R).replace('- -', '+ '), z[0], z[1], z[2], n1, n2, n3]
 
print('Окружность:', minimum[min(minimum)][0])
print('Три точки, определяющие окружность:', str(minimum[min(minimum)][1]) + ',',  str(minimum[min(minimum)][2]) + ',', minimum[min(minimum)][3])
print('Всего кол-во точек, принадлежащих окружности:', minimum[min(minimum)][6])
print('Кол-во точек внутри круга:', minimum[min(minimum)][4])
print('Кол-во точек вне круга:', minimum[min(minimum)][5])

data = [(1, 1), (-1, 3), (5, 3), (-3, -4), (6, 0), (1, 5), (-2, 0), (3, 3), (0, 0), (100, -100), (1, 100)]
minimum = {}
points = [[data[i], data[j], data[k]] for i in range(len(data)-2) for j in range(i+1, len(data)-1) for k in range(j+1, len(data))]
for z in points:
    x1, y1 = z[0]
    x2, y2 = z[1]
    x3, y3 = z[2]
    try:
        y0 = (2*(x1-x3)*(y2**2 - y1**2 - x1**2 + x2**2) + 2*(x1-x2)*(x1**2 - x3**2 - y3**2 + y1**2)) / (4*(y2-y1)*(x1-x3) - 4*(x1-x2)*(y3-y1))
    except ZeroDivisionError:
        continue
    if x1 - x3 != 0:
        x0 = ((x1**2 - x3**2 - y3**2 + y1**2) + 2*y0*(y3-y1))/ (2*(x1-x3))
    if x1 - x2 != 0:
        x0 = ((x1**2 - x2**2 - y2**2 + y1**2) + 2*y0*(y2-y1))/ (2*(x1-x2))
    R = ((x0 - x1)** 2 + (y0 - y1)** 2) ** (1/2)
 
    n1, n2, n3 = 0, 0, 0
 
    for t in data:
        if round((t[0] - x0) ** 2 + (t[1] - y0) ** 2, 3) < round(R ** 2, 3):
            n1 += 1
        elif round((t[0] - x0) ** 2 + (t[1] - y0) ** 2, 3) > round(R ** 2, 3):
            n2 += 1
        else:
            n3 += 1
    minimum[abs(n1 - n2)] = ['(x - {:.1f}) ** 2 + (y - {:.1f}) ** 2 = {:.1f} ** 2'.format(x0, y0, R).replace('- -', '+ '), z[0], z[1], z[2], n1, n2, n3]
 
print('Окружность:', minimum[min(minimum)][0])
print('Три точки, определяющие окружность:', str(minimum[min(minimum)][1]) + ',',  str(minimum[min(minimum)][2]) + ',', minimum[min(minimum)][3])
print('Всего кол-во точек, принадлежащих окружности:', minimum[min(minimum)][6])
print('Кол-во точек внутри круга:', minimum[min(minimum)][4])
print('Кол-во точек вне круга:', minimum[min(minimum)][5])