%Este programa calcula algunas matrices de Hilbert, su número 
% de condición y el producto del det de ella con el de de su inversa.
clear
for n=2:14
   for i=1:n
      for j=1:n
         a(i,j)=1/(i+j-1);
      end
   end
   c=cond(a);
   d=det(a)*det(a^(-1));
   fprintf('n=%3.0f cond(a)=%e det*det=%e\n', n,c,d)
   end