The vandermonde command solves a transposed vandermonde system. If an optional third argument P is provided, the computation is done modulo P.
A = [1,2,3]; B = [4,5,6]; vandermonde(A,B); [5/2, 2, -1/2]