Sub МНК_n()
Const n = 9, m = 5
Dim x(n) As Double, y(n) As Double, i As Integer, a(m, m) As Double, AA(m) As Double, p As Integer, q As Integer, _
b(m) As Double
For i = 1 To n
x(i) = Cells(i, 1)
y(i) = Cells(i, 2)
Next i
For q = 0 To m
For i = 1 To n
For p = 0 To m
a(q, p) = a(q, p) + x(i) ^ (q + p)
Next p
b(q) = b(q) + y(i) * x(i) ^ q
Next i
Next q
метод_Гаусса a, b, m, AA
For i = 1 To n
y(0) = 0
For p = 0 To m
y(0) = y(0) + AA(p) + x(i) ^ p
If i = 1 Then Debug.Print AA(p)
Next p
Cells(i, 6) = y(0)
Next i
End Sub

Sub метод_Гаусса(a() As Double, b() As Double, n As Integer, x() As Double)
Dim i As Integer, j As Integer, A1 As Double, k As Integer
For i = 0 To n
A1 = a(i, i)
b(i) = b(i) / A1
For j = 0 To n
a(i, j) = a(i, j) / A1
Next j
For k = i + 1 To n
A1 = a(k, i)
b(k) = b(k) - b(i) * A1
For j = 0 To n
a(k, j) = a(k, j) - a(i, j) * A1
Next j
Next k
Next i
x(n) = b(n)
For i = n - 1 To 0 Step -1
x(i) = b(i)
For j = n To i + 1 Step -1
x(i) = x(i) - x(j) * a(i, j)
Next j
Next i
End Sub