用Fortran或者VB语言写出用高斯积分法计算下式的程序段。

************线性方程组得解法********
'
Sub XXFCZ(ByRef A() As Double, ByRef B() As Double)
'高斯消元法
'A(i,j)是系数
'B(i)是右端项
Dim n As Long
Dim k As Long
Dim l As Double
Dim J As Long
Dim sum As Double
n = UBound(B)
For k = 1 To n - 1
For i = k + 1 To n
l = A(i, k) / A(k, k)
For J = k + 1 To n
A(i, J) = A(i, J) - l * A(k, J)
Next J
B(i) = B(i) - l * B(k)
Next i
Next k '以上是消元过程
B(n) = B(n) / A(n, n)
For i = n - 1 To 1 Step -1
sum = 0
For J = i + 1 To n
sum = sum + A(i, J) * B(J)
Next J
B(i) = (B(i) - sum) / A(i, i)
Next i '以上是回代过程
End Sub

下面的 高斯消元法就可以解决 方程系数为零的情况。它可以求解任意 线性方程。A（） 为方程系数，N为方程的未知数，B为方程右侧的常熟。
Function