midinf.txt

来自「用VB实现数值积分的算法集。包括梯形求积法」· 文本 代码 · 共 26 行

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