Algoritmo De Thomas

Const max = 20
Public fila As Integer
Public columna As Integer

Sub Cargar_Matriz(mat() As Double, vect() As Double, m As Integer)
Dim i As Integer, j As Integer
m = InputBox("Deme laDimencion de la matriz cuadrada(M x M):")
For i = 1 To m
fila = fila + 1
columna = 0
For j = 1 To m
mat(i, j) = InputBox("A[" + Str(i) + "," + Str(j) + "]=")columna = columna + 1
Cells(fila, columna) = mat(i, j)
Next j
columna = columna + 1
Cells(fila, columna) = "="
vect(i) = InputBox("b[" + Str(i) +"]=")
columna = columna + 1
Cells(fila, columna) = vect(i)
Next i
End Sub

Sub mostrar_Upper(u() As Double, m As Integer)
Dim i As Integer, j As Integer
For i = 1 To mfila = fila + 1
columna = 0
For j = 1 To m
If j >= i Then
columna = columna + 1
Cells(fila, columna) = u(i, j)
Elsecolumna = columna + 1
End If
Next j
Next i
End Sub

Sub Sustitucion_Arr_Abajo(l() As Double, vect() As Double, c() As Double, m As Integer)
Dim i As Integer, k As IntegerDim suma As Double
c(1) = vect(1)
For k = 2 To m
i = k - 1
c(k) = vect(k) - l(k, i) * c(k - 1)
Next i
End Sub

Sub Sustitucion_Abajo_Arr(mat() As Double, c() As Double,x() As Double, m As Integer)
Dim i As Integer, k As Integer
Dim suma As Double
x(m) = c(m) / mat(m, m)
For i = m - 1 To 1 Step -1

x(i) = (c(i) + x(i - 1)) / mat(i, i)Next i
End Sub

Sub Mostrar_solucion(x() As Double, m As Integer)
Dim i As Integer
For i = 1 To m
fila = fila + 1
Cells(fila, 1) = "x" + Str(i) + "=" + Str(x(i)) + ""
Next i
End Sub

Sub mostrar_thomas(mat() As Double, m As Integer)
Dim i As Integer, j As Integer
For i = 1 To m
fila = fila + 1
columna = 0
For j = 1 To i...