function exercise2c
  a = [2 3 5 10 15 25 50 100]';

  %n=30
  %a = ones(n,1)
  %a(1) = 2
  %for i=2:n
  %  a(i) = a(i) + a(i-1)
  %end
  
  
  for i=1:size(a,1)
    
    [H,h] = Hilbert(a(i));
    
    x = CG(H,h);
    cor = ones(size(H,1),1);
    n = norm(cor-x,'inf');    

    fprintf('%3d -> %1.20f\n',a(i),n)
  end  