Private Sub Command1_Click()
    'program D1R11
    'Driver program for routine QRBKSB,QRDCMP
    N = 5
    Dim A(5, 5), B(5), Q(5, 5), C(5, 5), R(5), X(5)
    '輸入已知的方程組的系數矩陣
    A(1, 1) = 1.4: A(1, 2) = 2.1: A(1, 3) = 2.1: A(1, 4) = 7.4: A(1, 5) = 9.6
    A(2, 1) = 1.6: A(2, 2) = 1.5: A(2, 3) = 1.1: A(2, 4) = 0.7: A(2, 5) = 5#
    A(3, 1) = 3.8: A(3, 2) = 8#:  A(3, 3) = 9.6: A(3, 4) = 5.4: A(3, 5) = 8.8
    A(4, 1) = 4.6: A(4, 2) = 8.2: A(4, 3) = 8.4: A(4, 4) = 0.4: A(4, 5) = 8#
    A(5, 1) = 2.6: A(5, 2) = 2.9: A(5, 3) = 0.1: A(5, 4) = 9.9: A(5, 5) = 7.7
    '輸入已知的方程組的右端向量B
    B(1) = 1.1
    B(2) = 1.6
    B(3) = 4.7
    B(4) = 9.1
    B(5) = 0.1
    Print
    Print Tab(5); "已知的方程組的右端向量"
    Print Tab(14); Format$(B(1), "##.##")
    Print Tab(14); Format$(B(2), "##.##")
    Print Tab(14); Format$(B(3), "##.##")
    Print Tab(14); Format$(B(4), "##.##")
    Print Tab(14); Format$(B(5), "##.##")
    For I = 1 To N
        For J = 1 To N
           C(I, J) = A(I, J)
        Next J
    Next I
    Call QRDCMP(C(), N, N, Q())
    For I = 1 To N
        R(I) = B(I)
    Next I
    Call QRBKSB(C(), N, Q(), R(), X())
    Print
    Print Tab(5); "計算出的方程組的解"
    Print Tab(10); Format$(X(1), "#.####E+00")
    Print Tab(10); Format$(X(2), "#.####E+00")
    Print Tab(10); Format$(X(3), "#.####E+00")
    Print Tab(10); Format$(X(4), "#.####E+00")
    Print Tab(10); Format$(X(5), "#.####E+00")
    '將計算出的解乘以系數矩陣，以驗證計算結果正確
    For L = 1 To N
        B(L) = 0#
        For J = 1 To N
            B(L) = B(L) + A(L, J) * X(J)
        Next J
    Next L
    Print
    Print Tab(5); "計算出的解乘以系數矩陣的結果"
    Print Tab(14); Format$(B(1), "##.##")
    Print Tab(14); Format$(B(2), "##.##")
    Print Tab(14); Format$(B(3), "##.##")
    Print Tab(14); Format$(B(4), "##.##")
    Print Tab(14); Format$(B(5), "##.##")
End Sub