d8r5.txt

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

Private Sub Command1_Click()
    'PROGRAM D8R5
    'Driver for routine BALANC
    NP = 5
    Dim A(5, 5), R(5), C(5)
    A(1, 1) = 1#: A(1, 2) = 100#: A(1, 3) = 1#: A(1, 4) = 100#: A(1, 5) = 1#
    A(2, 1) = 1#: A(2, 2) = 1#:   A(2, 3) = 1#: A(2, 4) = 1#:   A(2, 5) = 1#
    A(3, 1) = 1#: A(3, 2) = 100#: A(3, 3) = 1#: A(3, 4) = 100#: A(3, 5) = 1#
    A(4, 1) = 1#: A(4, 2) = 1#:   A(4, 3) = 1#: A(4, 4) = 1#:   A(4, 5) = 1#
    A(5, 1) = 1#: A(5, 2) = 100#: A(5, 3) = 1#: A(5, 4) = 100#: A(5, 5) = 1#
    'Print norms
    For I = 1 To NP
        R(I) = 0#
        C(I) = 0#
        For J = 1 To NP
            R(I) = R(I) + Abs(A(I, J))
            C(I) = C(I) + Abs(A(J, I))
        Next J
    Next I
    Print
    Print Tab(5); "Rows:"
    Print
    For I = 1 To NP
        Print Format$(R(I), "####.00000"),
    Next I
    Print Tab(5)
    Print Tab(5); " Columns:"
    Print
    For I = 1 To NP
        Print Format$(C(I), "####.00000"),
    Next I
    Print Tab(5)
    Print Tab(5); " ********** Balancing Matrix **********"
    Call BALANC(A(), NP)
    'Print norms
    For I = 1 To NP
        R(I) = 0#
        C(I) = 0#
        For J = 1 To NP
            R(I) = R(I) + Abs(A(I, J))
            C(I) = C(I) + Abs(A(J, I))
        Next J
    Next I
    Print Tab(5)
    Print Tab(5); "Rows:"
    Print
    For I = 1 To NP
        Print Format$(R(I), "####.00000"),
    Next I
    Print Tab(5)
    Print Tab(5); " Columns:"
    Print
    For I = 1 To NP
        Print Format$(C(I), "####.00000"),
    Next I
End Sub