Sub DDPOLY(C(), NC, X, PD(), ND)
    PD(1) = C(NC)
    For J = 2 To ND
        PD(J) = 0#
    Next J
    For I = NC - 1 To 1 Step -1
        If NC + 1 - I < ND Then
           NND = NC + 1 - I
        Else
           NND = ND
        End If
        For J = NND To 2 Step -1
            PD(J) = PD(J) * X + PD(J - 1)
        Next J
        PD(1) = PD(1) * X + C(I)
    Next I
    CONST1 = 2#
    For I = 3 To ND
        PD(I) = CONST1 * PD(I)
        CONST1 = CONST1 * I
    Next I
End Sub