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
for j=1:20
    plot(t,A(:,j))
    hold on
end