[VBA]各位老大我改了一段矩阵求逆的程序~请帮忙改正一下谢谢

chengxu_84

Dim N As Integer
Dim ch As Double, Z As Double, zz As Double
Dim P As Integer
Dim qn(1 To 6, 1 To 6) As Double
Dim nb(1 To 6, 1 To 6) As Double '返回的逆阵

 '对mc(1to 6,1 to 6)求逆

Dim E(1 To 6, 1 To 12) As Double
 '应用的原理是: [ A I]经过矩阵行变换后成为[ I ivt_A], E阵(N,2N)阶 既 [ A I]
For i = 1 To 6
  For j = 1 To 6
E(i, j) = mc(i, j)
  Next j
Next i
For i = 1 To 6
  For j = 7 To 12
If i = j - 6 Then
E(i, j) = 1
Else
E(i, j) = 0
End If
  Next j
Next i
'以上步骤实现用A，I初始化E阵

For j = 1 To 6
If E(j, j) = 0 Then
For i = 1 To 6
If E(i, j) <> 0 Then Exit For
Next i
For P = 1 To 12
ch = E(i, P)
E(i, P) = E(j, P)
E(j, P) = ch
Debug.Print E(i, P)
Next P '将对角线元素为0的行与对应列上元素非零的行，进行行变换

End If

Z = E(j, j)
For P = 1 To 12
E(j, P) = E(j, P) / Z '将矩阵E （A，I） A矩阵每行都除以对角线元素,使对角线元素变为1
Next P

For i = 1 To 6
If i <> j And E(i, j) <> 0 Then
zz = E(i, j)

For P = 1 To 12
E(i, P) = E(i, P) - E(j, P) * zz
Next P
End If
 Next i
Next j

For i = 1 To 6
 For j = 1 To 6
nb(i, j) = E(i, j + 6)
 
 Next j
Next i
 
  ' qn(i, j) = nb(i, j) '返回求出的A阵的逆阵B阵