d1r7.txt

VB常用数值算法集2

Private Sub Command1_Click()
    'program D1R7
    'Driver program for routine TOEPLZ
    N = 5
    N2 = 2 * N
    Dim X(5), Y(5), R(10)
    '输入已知的方程组的系数矩阵的R
    For I = 1 To 2 * N - 1
    R(I) = 1 / I
    Next I
    '输入已知的方程组的右端向量
    For I = 1 To N
        Y(I) = 0.1 * I
    Next I
    Print Tab(5); "已知的方程组的右端向量"
    Print Tab(12); Format$(Y(1), "##.##00")
    Print Tab(12); Format$(Y(2), "##.##00")
    Print Tab(12); Format$(Y(3), "##.##00")
    Print Tab(12); Format$(Y(4), "##.##00")
    Print Tab(12); Format$(Y(5), "##.##00")
    Call TOEPLZ(R(), X(), Y(), N)
    '输出方程组的解X
    Print
    Print Tab(5); "计算出的方程组的的解"
    Print Tab(9); Format$(X(1), ".#####0e+00")
    Print Tab(9); Format$(X(2), ".#####0e+00")
    Print Tab(9); Format$(X(3), ".#####0e+00")
    Print Tab(9); Format$(X(4), ".#####0e+00")
    Print Tab(9); Format$(X(5), ".#####0e+00")
    '将计算出的解X乘以系数矩阵，以验证计算结果正确
    Print
    Print Tab(4); "将计算出的解乘以系数矩阵，以验证计算结果正确"
    Print
    Print Tab(5); "解乘以系数矩阵"; Tab(24); "方程组的右端向量"
    For I = 1 To N
        SUM1 = 0
        For J = 1 To N
           SUM1 = SUM1 + R(N + I - J) * X(J)
        Next J
        Print Tab(10); Format$(SUM1, "#.####0"); Tab(30); Format$(Y(I), "#.###0")
    Next I
End Sub