求助!!改进的平方根方法求对称正定线性方程组
求助!!改进的平方根方法求对称正定线性方程组
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
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