chebpc.txt

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

Sub CHEBPC(C(), D(), N)
    Dim DD(50)
    For J = 1 To N
        D(J) = 0#
        DD(J) = 0#
    Next J
    D(1) = C(N)
    For J = N - 1 To 2 Step -1
        For K = N - J + 1 To 2 Step -1
            SV = D(K)
            D(K) = 2# * D(K - 1) - DD(K)
            DD(K) = SV
        Next K
        SV = D(1)
        D(1) = -DD(1) + C(J)
        DD(1) = SV
    Next J
    For J = N To 2 Step -1
        D(J) = D(J - 1) - DD(J)
    Next J
    D(1) = -DD(1) + 0.5 * C(1)
End Sub