d1r6.txt

VB常用数值算法集2

Private Sub Command1_Click()
    'program D1R6
    'Driver program for routine VANDER
    N = 5
    Dim X(5), Q(5), W(5), TERM(5)
    '输入已知的方程组的系数矩阵的X
    X(1) = 1#: X(2) = 1.5: X(3) = 2#: X(4) = 2.5: X(5) = 3#
    '输入已知的方程组的右端向量B
    Q(1) = 1#: Q(2) = 1.5: Q(3) = 2#: Q(4) = 2.5: Q(5) = 3#
    Print Tab(5); " 已知的方程组的右端向量"
    Print Tab(13); Format$(Q(1), "##.##")
    Print Tab(13); Format$(Q(2), "##.##")
    Print Tab(13); Format$(Q(3), "##.##")
    Print Tab(13); Format$(Q(4), "##.##")
    Print Tab(13); Format$(Q(5), "##.##")
    Call VANDER(X(), W(), Q(), N)
    '输出方程组的解W
    Print
    Print Tab(5); " 计算出的方程组的的解"
    Print Tab(11); Format$(W(1), "##.######")
    Print Tab(11); Format$(W(2), "##.######")
    Print Tab(11); Format$(W(3), "##.######")
    Print Tab(11); Format$(W(4), "##.######")
    Print Tab(11); Format$(W(5), "##.######")
    '将计算出的解乘以系数矩阵，以验证计算结果正确
    Print
    Print Tab(4); "将计算出的解乘以系数矩阵，以验证计算结果正确"
    Print
    Print Tab(7); "解乘以系数矩阵"; Tab(26); "方程组的右端向量"
    SUM1 = 0
    For I = 1 To N
        TERM(I) = W(I)
        SUM1 = SUM1 + W(I)
    Next I
    Print Tab(12); Format$(SUM1, "#.####"); Tab(32); Format$(Q(1), "#.####")
    For I = 2 To N
        SUM1 = 0
        For J = 1 To N
           TERM(J) = TERM(J) * X(J)
           SUM1 = SUM1 + TERM(J)
        Next J
        Print Tab(12); Format$(SUM1, "#.####"); Tab(32); Format$(Q(I), "#.####")
    Next I
End Sub