m=101;
n=20;
t=linspace(-1,1,m);
y=rand(1,m);
c=polyfit(t,y,n-1);
A=zeros(m,n);
y=y';
for i=1:n
    for j=1:m
          A(j,i)=y(j)^(i-1);

    end
end
[U,S,V]=svd(A);
[Q,R]=qr(A);
plot(U(:,2))
hold on
plot(Q(:,2))