
function demo_hankel

ns=1:25;

for n=ns
    A=hankel_matrix(n);
    xsol=ones(size(A,1),1);
    b=A*xsol;
    x_LU=solve_linear_LU(A,b);
    x_LUP=solve_linear_LUP(A,b);
    
    e_LU=norm(xsol-x_LU,inf);
    e_LUP=norm(xsol-x_LUP,inf);
    
    fprintf('\n $%2.0f$ & $%1.2e$ & $%1.2e$',n,e_LU,e_LUP)
end