d8r7.txt

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

Private Sub Command1_Click()
    'PROGRAM D8R7
    'Driver for routine HQR
    NP = 5
    Dim A(5, 5), WR(5), WI(5)
    A(1, 1) = 1#:  A(1, 2) = 2#: A(1, 3) = 0#:   A(1, 4) = 0#: A(1, 5) = 0#
    A(2, 1) = -2#: A(2, 2) = 3#: A(2, 3) = 0#:   A(2, 4) = 0#: A(2, 5) = 0#
    A(3, 1) = 3#:  A(3, 2) = 4#: A(3, 3) = 50#:  A(3, 4) = 0#: A(3, 5) = 0#
    A(4, 1) = -4#: A(4, 2) = 5#: A(4, 3) = -60#: A(4, 4) = 7#: A(4, 5) = 0#
    A(5, 1) = -5#: A(5, 2) = 6#: A(5, 3) = -70#: A(5, 4) = 8#: A(5, 5) = -9#
    Print Tab(5); " Matrix:"
    Print Tab(5)
    For I = 1 To NP
        Print Tab(5);
        For J = 1 To NP
            Print Format$(A(I, J), "##.00"),
        Next J
        Print
    Next I
    Call BALANC(A(), NP)
    Call ELMHES(A(), NP)
    Call HQR(A(), NP, WR(), WI())
    Print Tab(5)
    Print Tab(5); "Eigenvalues:"
    Print Tab(5)
    Print Tab(5); "    Real          Imag."
    For I = 1 To NP
        Print Tab(5); Format$(WR(I), ".000000E+00");
        Print Tab(19); Format$(WI(I), ".000000E+00")
    Next I
End Sub