chint.txt

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

Sub CHINT(A, B, C(), CINT1(), N)
    CON = 0.25 * (B - A)
    Sum = 0#
    FAC = 1#
    For J = 2 To N - 1
        CINT1(J) = CON * (C(J - 1) - C(J + 1)) / (J - 1)
        Sum = Sum + FAC * CINT1(J)
        FAC = -FAC
    Next J
    CINT1(N) = CON * C(N - 1) / (N - 1)
    Sum = Sum + FAC * CINT1(N)
    CINT1(1) = 2# * Sum
End Sub