Решение СЛАУ методом Гаусса

Dim a(10, 10) As Double, b(10) As Double, x(10) As Double, tmp As Double
Dim i As Integer, j As Integer, k As Integer, l As Integer, n As Integer
Private Sub Command1_Click()
'Метод Гаусса
'Прямой проход
For k = 1 To n - 1
For i = k + 1 To n
tmp = a(i, k)
b(i) = b(i) - b(k) * tmp / a(k, k)
For j = 1 To n    '33
a(i, j) = a(i, j) - a(k, j) * tmp / a(k, k)
Next
Next
'Проверка на 0 диагонального элемента
If a(k + 1, k + 1) = 0 Then
'Поиск не нулевого a(k + 1, k + 1)
For i = k + 2 To n
If a(i, k + 1) <> 0 Then
l = i
Exit For
End If
Next
'Перестановка строк
For j = 1 To n
tmp = a(k + 1, j)
a(k + 1, j) = a(l, j)
a(l, j) = tmp
tmp = b(k + 1)
b(k + 1) = b(l)
b(l) = tmp
Next
End If
Next
'Обратный проход
x(n) = b(n) / a(n, n)
For i = n - 1 To 1 Step -1
x(i) = b(i)
For j = i + 1 To n
x(i) = x(i) - a(i, j) * x(j)
Next
x(i) = x(i) / a(i, i)
Next
For j = 1 To n
MsgBox "x(" & Str(j) & ")=" & Str(x(j))
Next
End Sub
 
Private Sub Form_Load()
n = 4
'Матрица коэфициентов
a(1, 1) = 2
a(1, 2) = -4
a(1, 3) = 1
a(1, 4) = -5
 
a(2, 1) = 4
a(2, 2) = -7
a(2, 3) = -1
a(2, 4) = -8
 
a(3, 1) = 10
a(3, 2) = -18
a(3, 3) = 2
a(3, 4) = -23
 
a(4, 1) = 2
a(4, 2) = -3
a(4, 3) = 1
a(4, 4) = -1
 
b(1) = 2
b(2) = -5
b(3) = 3
b(4) = 0
End Sub