function V  =  VandermondeMatrix(x)

n  =  size(x,1)-1;
V  =  zeros(n+1);

for i=1:n+1
    for j=1:n+1
        V(i,j)   =  x(i)^(j-1)   ;
    end
end