function [error,kappa,H,L,U]=hilb_test(n)
xexact=ones(n,1);
H=hilb(n);
b=H*xexact;
[L,U]=lu(H);
y=L\b;
x=U\y;
error=norm(x-xexact,1)/n;
kappa=cond(H,1);