setup_and_solve_Vandermonde
xpts = linspace(-5,5,200)
for i = 1:200
    z1(i) = interp_poly_val(a,xpts(i));
end
z2 = 1./(1+xpts.^2);
plot(xpts,z1,xpts,z2)   
    