% Define the polynomial degree and end points of the interval --
degree = 10;
n = degree+1;
a = -5;
b = 5;
h = (b-a)/(n-1);
% Construct the abscissas --
for i =1:n;
    x(i) = a + (i-1)*h;
end
% Evaluate the function at the abscissas --
y = 1./(1+x.^2);
y=y';
% Set up the Vandermonde system --
A = zeros(n,n);
A(:,1) = 1;
for j=2:n
    for i=1:n
        A(i,j) = x(i)^(j-1);
    end
end
a = A\y;




