d1r5.txt

VB常用数值算法集2

Private Sub Command1_Click()
    'program D1R5
    'Driver program for routine MPROVE
    N = 5
    Dim A(5, 5), A1(5, 5), B(5), INDX(5), B1(5)
    '输入已知的方程组的系数矩阵
    A(1, 1) = 1: A(1, 2) = 2: A(1, 3) = 3: A(1, 4) = 4: A(1, 5) = 5
    A(2, 1) = 2: A(2, 2) = 3: A(2, 3) = 4: A(2, 4) = 5: A(2, 5) = 1
    A(3, 1) = 1: A(3, 2) = 1: A(3, 3) = 1: A(3, 4) = 1: A(3, 5) = 1
    A(4, 1) = 4: A(4, 2) = 5: A(4, 3) = 1: A(4, 4) = 2: A(4, 5) = 3
    A(5, 1) = 5: A(5, 2) = 1: A(5, 3) = 2: A(5, 4) = 3: A(5, 5) = 4
    '输入已知的方程组的右端向量B
    B(1) = 1
    B(2) = 1
    B(3) = 1
    B(4) = 1
    B(5) = 1
    For I = 1 To N
        For J = 1 To N
            A1(I, J) = A(I, J)
        Next J
    Next I
    For I = 1 To N
        B1(I) = B(I)
    Next I
    Call LUDCMP(A(), N, INDX(), D)
    Call LUBKSB(A(), N, INDX(), B())
    '输出方程组的解B
    Print
    Print Tab(5); "输出方程组的解"
    Print Tab(10); Format$(B(1), "##.####")
    Print Tab(10); Format$(B(2), "##.####")
    Print Tab(10); Format$(B(3), "##.####")
    Print Tab(10); Format$(B(4), "##.####")
    Print Tab(10); Format$(B(5), "##.####")
    For I = 1 To N
        B(I) = B(I) * (1# + 0.2 * Rnd(I))
    Next I
    '输出干扰后的解B
    Print
    Print Tab(5); "输出干扰后的解"
    Print Tab(8); Format$(B(1), "##.####")
    Print Tab(8); Format$(B(2), "##.####")
    Print Tab(8); Format$(B(3), "##.####")
    Print Tab(8); Format$(B(4), "##.####")
    Print Tab(8); Format$(B(5), "##.####")
    Call MPROVE(A1(), A(), N, INDX(), B1(), B())
    '输出改善后的解B
    Print
    Print Tab(5); "输出改善后的解"
    Print Tab(10); Format$(B(1), "##.####")
    Print Tab(10); Format$(B(2), "##.####")
    Print Tab(10); Format$(B(3), "##.####")
    Print Tab(10); Format$(B(4), "##.####")
    Print Tab(10); Format$(B(5), "##.####")
End Sub