Есть множество точек в которое необходимо вписать окружность и определить центр окружности и радиус.

Private Sq As Double, Par(1 To 3) As Double
Private I As Integer, J As Integer, N As Integer
Private XY() As Variant
Sub BegI()
Dim R As Double, A As Double, B As Double, K As Integer, L As Integer
Dim Mt(1 To 3, 1 To 4) As Double, Fn() As Double, Fpl() As Double
Dim Xmin As Double, Ymin As Double, Xmax As Double, Ymax As Double
R = 6: A = 4: B = 4: L = 0
N = 26
ReDim XY(2, N)
'формируем модельный пример окружности со случайными отклонениями от нее
For I = -2 To 10
  XY(1, L + 1) = I: XY(1, L + 2) = I
  Sq = Sqr(R ^ 2 - (I - A) ^ 2)
  XY(2, L + 1) = (B + Sq) * (1 + (0.2 * (Rnd - 0.5)))
  XY(2, L + 2) = (B - Sq) * (1 + (0.2 * (Rnd - 0.5)))
  L = L + 2
Next
Xmin = 1000000: Ymin = 1000000: Xmax = 0: Ymax = 0
For I = 1 To N
  If Xmin > XY(1, I) Then Xmin = XY(1, I)
  If Ymin > XY(2, I) Then Ymin = XY(2, I)
  If Xmax < XY(1, I) Then Xmax = XY(1, I)
  If Ymax < XY(2, I) Then Ymax = XY(2, I)
Next
Par(2) = (Xmin + Xmax) / 2: Par(3) = (Ymin + Ymax) / 2
 
'начальное приближение для радиуса - половина длинной стороны
Xmin = (Xmax - Xmin) / 2: Ymin = (Ymax - Ymin) / 2
Par(1) = (IIf(Xmin > Ymin, Xmin, Ymin)) ^ 2
'Если вы координаты точек введете в лист Экселя, то
'вместо модельного примера поставьте следующие операторы
'N = число ваших точек
'XY = Range("A1:B") & N
'
 
L = 3 'ЧИСЛО УРАВНЕНИЙ
'формирование матрицы 3 на 4 для расчета по Ньютону-Рафсону
Beg:
Fn = FuncI
For I = 1 To L  'считаем производные
  Par(I) = Par(I) + 0.001
  Fpl = FuncI
  For J = 1 To L
    Mt(J, I) = 1000 * (Fpl(J) - Fn(J))
  Next
  Par(I) = Par(I) - 0.001
  Mt(I, L + 1) = Fn(I)
Next
'решение системы уравнений для уточнения значений параметров
For I = 1 To L
  Sq = Mt(I, I)
  For J = I To L + 1
    Mt(I, J) = Mt(I, J) / Sq
  Next
  For K = I + 1 To L
    For J = I + 1 To L + 1
      Mt(K, J) = Mt(K, J) - Mt(K, I) * Mt(I, J)
    Next
  Next
Next
For I = L To 1 Step -1
  For J = I - 1 To 1 Step -1
    Mt(J, L + 1) = Mt(J, L + 1) - Mt(I, L + 1) * Mt(J, I)
  Next
  Par(I) = Par(I) - Mt(I, L + 1)
Next
MsgBox "R=" & Sqr(Par(1)) & " a=" & Par(2) & " b=" & Par(3)
If Abs(Mt(1, L + 1)) > 0.0001 Then GoTo Beg
End Sub
Function FuncI() As Double()
Dim F(1 To 3) As Double, J As Integer
F(1) = 0: F(2) = 0: F(3) = 0
For J = 1 To N
  Sq = Par(1) - (XY(1, J) - Par(2)) ^ 2 - (XY(2, J) - Par(3)) ^ 2
  F(1) = F(1) + Sq: F(2) = F(2) + Sq * XY(1, J): F(3) = F(3) + Sq * XY(2, J)
Next
FuncI = F
End Function