d8r6.txt

来自「VB常用数值算法集 内含有解线性代数方程组」· 文本 代码 · 共 45 行

Private Sub Command1_Click()
    'PROGRAM D8R6
    'Driver for routine ELMHES
    NP = 5
    Dim A(5, 5), R(5), C(5)
    A(1, 1) = 1#: A(1, 2) = 2#: A(1, 3) = 300#: A(1, 4) = 4#: A(1, 5) = 5#
    A(2, 1) = 2#: A(2, 2) = 3#: A(2, 3) = 400#: A(2, 4) = 5#: A(2, 5) = 6#
    A(3, 1) = 3#: A(3, 2) = 4#: A(3, 3) = 5#:   A(3, 4) = 6#: A(3, 5) = 7#
    A(4, 1) = 4#: A(4, 2) = 5#: A(4, 3) = 600#: A(4, 4) = 7#: A(4, 5) = 8#
    A(5, 1) = 5#: A(5, 2) = 6#: A(5, 3) = 700#: A(5, 4) = 8#: A(5, 5) = 9#
    Print Tab(5)
    Print Tab(5); " ***** Original Matrix *****"
    Print Tab(5)
    For I = 1 To NP
        For J = 1 To NP
            Print Format$(A(I, J), "#.00"),
        Next J
        Print
    Next I
    Print Tab(5)
    Print Tab(5); "***** Balance Matrix *****"
    Print Tab(5)
    Call BALANC(A(), NP)
    For I = 1 To NP
        For J = 1 To NP
            Print Format$(A(I, J), "#.00"),
        Next J
        Print
    Next I
    Print Tab(5)
    Print Tab(5); " *****Reduce to Hessenherg Form *****"
    Print Tab(5)
    Call ELMHES(A(), NP)
    For J = 1 To NP - 2
        For I = J + 2 To NP
            A(I, J) = 0#
        Next I
    Next J
    For I = 1 To NP
        For J = 1 To NP
            Print Format$(A(I, J), ".0000E+00"),
        Next J
        Print
    Next I
End Sub