def DFT_Vandermonde(P,n,x):
    z=zero_vector(n)
    for k in range(0,n):
        z[k]=P[k]
    z2=z.column()
    return(Vandermonde(x,n)*z2)


show(DFT_Vandermonde(X^2+X+3,3,4))