求助！！改进的平方根方法求对称正定线性方程组  
Private Sub gjpfgff(a() As Double, b() As Double, n As Integer)
'改进的平方根方法求对称正定线性方程组
'n是未知量的个数
'a()数组存放方程组的系数矩阵,b()数组存放方程组的右端常数项
'程序执行完毕后，b()数组中的n个数就是方程组的解,a()数组的数已不再是原方程组的系数矩阵
    Dim v() As Double
    ReDim v(n)
    Dim i As Integer, p As Integer, k As Integer, j As Integer
    
    For j = 1 To n
        For k = 1 To j - 1
            v(k) = a(k, k) * a(j, k)
        Next k
        For k = 1 To j - 1
            a(j, j) = a(j, j) - a(j, k) * v(k)
        Next k
        For i = j + 1 To n
            For k = 1 To j - 1
                a(i, j) = a(i, j) - a(i, k) * v(k)
            Next k
            a(i, j) = a(i, j) / a(j, j)
        Next i
    Next j
    
    For i = 1 To n
        For k = i + 1 To n
            a(i, k) = a(k, i) * a(i, i)
        Next k
    Next i
    
    For i = 1 To n
        For j = 1 To i - 1
            b(i) = b(i) - b(j) * a(i, j)
        Next j
    Next i
        
    For i = n To 1 Step -1
        For j = i + 1 To n
            b(i) = b(i) - b(j) * a(i, j)
        Next j
        b(i) = b(i) / a(i, i)
    Next i
End Sub
以上就是，具体的资料了，VB我是不会了，也看不懂，懂的大哥，大姐，帮忙下，帮我改成fortran90地发出来或e-mail:dodozo@163.com