d10r11.txt

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

Private Sub Command1_Click()
    'PROGRAM D10R11
    'Driver for routine ZROOTS
    M = 4
    M1 = M + 1
    Dim A(2, 5), X, ROOTS(2, 4)
    For J = 1 To M1
        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
    POLISH% = 0
    Call ZROOTS(A(), M, ROOTS(), POLISH%)
    Print Tab(5); "Unpolished roots:"
    Print Tab(3); "Root #        Real          Imag."
    For I = 1 To M
        Print Tab(5); Format$(I, "##");
        Print Tab(15); Format$(ROOTS(1, I), "#.####00");
        Print Tab(29); Format$(ROOTS(2, I), "#.####00")
    Next I
    Print Tab(5)
    Print Tab(5); "Corrupted roots:"
    For I = 1 To M
        ROOTS(1, I) = ROOTS(1, I) * (1# + 0.01 * I)
        ROOTS(2, I) = ROOTS(2, I) * (1# + 0.01 * I)
    Next
    Print Tab(3); "Roots #       Real          Imag."
    For I = 1 To M
        Print Tab(5); Format$(I, "##");
        Print Tab(15); Format$(ROOTS(1, I), "#.#####0");
        Print Tab(29); Format$(ROOTS(2, I), "#.#####0")
    Next I
    POLISH% = -1
    Call ZROOTS(A(), M, ROOTS(), POLISH%)
    Print Tab(5)
    Print Tab(5); "Polished roots:"
    Print Tab(3); "Roots #       Real          Imag."
    For I = 1 To M
        Print Tab(5); Format$(I, "##");
        Print Tab(15); Format$(ROOTS(1, I), "#.#####0");
        Print Tab(29); Format$(ROOTS(2, I), "#.#####0")
    Next I
End Sub