trapzd.txt

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

Sub TRAPZD(A, B, S, N)
    If N = 1 Then
        S = 0.5 * (B - A) * (FUNC(A) + FUNC(B))
        IT = 1
    Else
        IT = 2 ^ (N - 2)
        TNM = IT
        DEL = (B - A) / TNM
        X = A + 0.5 * DEL
        Sum = 0#
        For J = 1 To IT
            Sum = Sum + FUNC(X)
            X = X + DEL
        Next J
        S = 0.5 * (S + (B - A) * Sum / TNM)
    End If
End Sub