d8r2.txt

VB常用数值算法集 内含有解线性代数方程组

Private Sub Command1_Click()
    'PROGRAM D8R2
    'Driver for routine EIGSRT
    NP = 3
    Dim D(3), V(3, 3), A(3, 3)
    A(1, 1) = 1#: A(1, 2) = 2#: A(1, 3) = 3#
    A(2, 1) = 2#: A(2, 2) = 2#: A(2, 3) = 3#
    A(3, 1) = 3#: A(3, 2) = 3#: A(3, 3) = 3#
    Call JACOBI(A(), NP, D(), V(), NROT)
    Print Tab(5); "Unsorted Eigenvectors:"
    For I = 1 To NP
        Print Tab(5)
        Print Tab(5); " Eigenvalue "; I; " = ";
        Print Format$(D(I), ".#####0")
        Print Tab(5); " Eigenvector:"
        For J = 1 To NP
            Print Format$(V(J, I), "#.######"),
        Next J
    Next I
    Print Tab(5)
    Print Tab(5); "********** sorting **********"
    Print
    Call EIGSRT(D(), V(), NP)
    Print Tab(5); " Sorted Eigenvectors:"
    For I = 1 To NP
        Print Tab(5)
        Print Tab(5); " Eigenvalue "; I; " = ";
        Print Format$(D(I), ".#####0")
        Print Tab(5); " Eigenvector:"
        For J = 1 To NP
            Print Format$(V(J, I), ".#####0"),
        Next J
    Next I
End Sub