qtrap.txt

来自「本书共有数值计算中常用的VisualBasic子过程近200个」· 文本 代码 · 共 11 行

Sub QTRAP(A, B, S)
    EPS = 0.000001
    JMAX = 20
    OLDS = -1E+30
    For J = 1 To JMAX
        Call TRAPZD(A, B, S, J)
        If Abs(S - OLDS) < EPS * Abs(OLDS) Then Exit Sub
        OLDS = S
    Next J
    Print "Too many steps."
End Sub