Sub PCSHFT(A, B, D(), N)
    CONST1 = 2# / (B - A)
    FAC = CONST1
    For J = 2 To N
        D(J) = D(J) * FAC
        FAC = FAC * CONST1
    Next J
    CONST1 = 0.5 * (A + B)
    For J = 1 To N - 1
        For K = N - 1 To J Step -1
            D(K) = D(K) - CONST1 * D(K + 1)
        Next K
    Next J
End Sub