d10r10.txt

VB常用数值算法集2

Private Sub Command1_Click()
    'PROGRAM D10R10
    'Driver for routine LAGUER
    M = 4
    MP1 = M + 1
    NTRY = 21
    EPS = 0.000001
    Dim A(2, 5), Y(2, 21), X(2)
    For J = 1 To MP1
        For I = 1 To 2
            A(I, J) = 0
        Next I
    Next J
    A(2, 1) = 2#
    A(1, 3) = -1#
    A(2, 3) = -2#
    A(1, 5) = 1#
    Print
    Print Tab(5); "Roots of polynomial x^4-(1+2i)*x^2+2i"
    Print
    Print Tab(5); "             Real         Complex"
    Print
    N = 0
    POLISH% = 0
    For I = 1 To NTRY
        X(1) = (I - 11#) / 10#
        X(2) = (I - 11#) / 10#
        Call LAGUER(A(), M, X(), EPS, POLISH%)
        If N = 0 Then
            N = 1
            Y(1, 1) = X(1)
            Y(2, 1) = X(2)
            Print Tab(5); Format$(N, "#");
            Print Tab(15); Format$(X(1), "#.#####0");
            Print Tab(29); Format$(X(2), "#.#####0")
        Else
            IFLAG = 0
            For J = 1 To N
                AAA = CABS(X(1) - Y(1, J), X(2) - Y(2, J))
                BBB = EPS * CABS(X(1), X(2))
                If AAA <= BBB Then IFLAG = 1
            Next J
            If IFLAG = 0 Then
                N = N + 1
                Y(1, N) = X(1)
                Y(2, N) = X(2)
                Print Tab(5); Format$(N, "#");
                Print Tab(15); Format$(X(1), "#.#####0");
                Print Tab(29); Format$(X(2), "#.#####0")
            End If
        End If
    Next I
End Sub