function Vk=makeV(x,k)
% function Vk=makeV(x,k)
% input x: sequence of n values as column vector
% k: degree of polynomial
% output Vk: 

[m,n]=size(x);
if m<n
x=x';
m=n;
end

Vk=zeros(m,k+1);

for col=1:1:(k+1)
   Vk(:,col)=x.^(col-1);
end


   

