d9r3.txt

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

Private Sub Command1_Click()
    'PROGRAM D9R3
    'Driver for routine COVSRT
    MA = 10
    MFIT = 5
    Dim COVAR(10, 10), LISTA(5)
    For I = 1 To MA
        For J = 1 To MA
            COVAR(I, J) = 0#
            If I <= 5 And J <= 5 Then
                COVAR(I, J) = I + J - 1
            End If
        Next J
    Next I
    Print Tab(5); "Original matrix"
    For I = 1 To MA
        For J = 1 To MA
            Print Tab(J * 6); Format$(COVAR(I, J), "#.0");
        Next J
    Next I
    'Test 1 - spread by 2
    Print Tab(5); "Test #1 - Spread by two"
    For I = 1 To MFIT
        LISTA(I) = 2 * I
    Next I
    Call COVSRT(COVAR(), MA, MA, LISTA(), MFIT)
    For I = 1 To MA
        For J = 1 To MA
            Print Tab(J * 6); Format$(COVAR(I, J), "#.0");
        Next J
    Next I
    'Test 2 - reverse
    Print Tab(5); "Test #2 - Reverse"
    For I = 1 To MA
        For J = 1 To MA
            COVAR(I, J) = 0#
            If I <= 5 And J <= 5 Then
                COVAR(I, J) = I + J - 1
            End If
        Next J
    Next I
    For I = 1 To MFIT
        LISTA(I) = MFIT + 1 - I
    Next I
    Call COVSRT(COVAR(), MA, MA, LISTA(), MFIT)
    For I = 1 To MA
        For J = 1 To MA
            Print Tab(J * 6); Format$(COVAR(I, J), "#.0");
        Next J
    Next I
    'Test 3 - spread and reverse
    Print Tab(5); "Test #3 - Spread and reverse"
    For I = 1 To MA
        For J = 1 To MA
            COVAR(I, J) = 0#
            If I <= 5 And J <= 5 Then
                COVAR(I, J) = I + J - 1
            End If
        Next J
    Next I
    For I = 1 To MFIT
        LISTA(I) = MA + 2 - 2 * I
    Next I
    Call COVSRT(COVAR(), MA, MA, LISTA(), MFIT)
    For I = 1 To MA
        For J = 1 To MA
            Print Tab(J * 6); Format$(COVAR(I, J), "#.0");
        Next J
    Next I
End Sub