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

    end
end
[U,S,V]=svd(A);
b=U'*y;

plot(b)