qromo.txt

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

Sub QROMO(A, B, SS, PICK$)
    EPS = 0.00003
    JMAX = 14
    JMAXP = JMAX + 1
    K = 7
    Dim S(15), H(15)
    H(1) = 1#
    For J = 1 To JMAX
        If PICK$ = "MIDPNT" Then Call MIDPNT(A, B, S(J), J)
        If PICK$ = "MIDINF" Then Call MIDINF(A, B, S(J), J)
        If PICK$ = "MIDSQL" Then Call MIDSQL(A, B, S(J), J)
        If PICK$ = "MIDSQU" Then Call MIDSQU(A, B, S(J), J)
        If J > K Then
            Call POLINT(H(), S(), K, 0#, SS, DSS)
            If Abs(DSS) < EPS * Abs(SS) Then Exit Sub
        End If
        S(J + 1) = S(J)
        H(J + 1) = H(J) / 9#
    Next J
    Print "Too many steps."
End Sub