function[R]=Hilbert(N)
A=hilb(N);
B=inv(A);
mu=norm(A,2)*norm(B,2);
mu;
H=A;
n=size(H,1);
b=zeros(n,1);
sum=0;
for i=1:n
    for j=1:n
        sum=sum+H(i,j);
    end
    b(i)=sum;
    sum=0;
end
x1=ones(N,1)
xstar=H\b
err=x1-xstar
R=norm(err,inf);