GNU Octave, version 2.9.8 (i486-pc-linux-gnu).
Copyright (C) 2006 John W. Eaton.
This is free software; see the source code for copying conditions.
There is ABSOLUTELY NO WARRANTY; not even for MERCHANTIBILITY or
FITNESS FOR A PARTICULAR PURPOSE.  For details, type `warranty'.

Additional information about Octave is available at http://www.octave.org.

Please contribute if you find this software useful.
For more information, visit http://www.octave.org/help-wanted.html

Report bugs to <bug@octave.org> (but first, please read
http://www.octave.org/bugs.html to learn how to write a helpful report).


n=5, cond(hilb(5))=476607

no piv.                   - |A*x-b| =               0 - rel.error [%] =     1.55304e-10
partial piv. (swap)       - |A*x-b| =               0 - rel.error [%] =     3.49344e-11
partial piv. (indirect)   - |A*x-b| =     1.11022e-16 - rel.error [%] =      3.4933e-11
scaled partial piv.       - |A*x-b| =               0 - rel.error [%] =     4.70864e-10

n=6, cond(hilb(6))=1.49511e+07

no piv.                   - |A*x-b| =               0 - rel.error [%] =     6.36985e-09
partial piv. (swap)       - |A*x-b| =               0 - rel.error [%] =     3.17407e-09
partial piv. (indirect)   - |A*x-b| =               0 - rel.error [%] =     3.17407e-09
scaled partial piv.       - |A*x-b| =     1.11022e-16 - rel.error [%] =     1.48292e-08

n=7, cond(hilb(7))=4.75367e+08

no piv.                   - |A*x-b| =     1.11022e-16 - rel.error [%] =     1.82942e-07
partial piv. (swap)       - |A*x-b| =               0 - rel.error [%] =     5.92343e-08
partial piv. (indirect)   - |A*x-b| =               0 - rel.error [%] =     5.92342e-08
scaled partial piv.       - |A*x-b| =     5.76889e-16 - rel.error [%] =     3.18237e-08

n=11, cond(hilb(11))=5.21823e+14

no piv.                   - |A*x-b| =     2.71948e-16 - rel.error [%] =        0.280551
partial piv. (swap)       - |A*x-b| =     2.48253e-16 - rel.error [%] =        0.085787
partial piv. (indirect)   - |A*x-b| =     6.66134e-16 - rel.error [%] =        0.085787
scaled partial piv.       - |A*x-b| =     6.18146e-16 - rel.error [%] =        0.374256

n=12, cond(hilb(12))=1.6523e+16

no piv.                   - |A*x-b| =     2.71948e-16 - rel.error [%] =         4.05758
partial piv. (swap)       - |A*x-b| =     2.48253e-16 - rel.error [%] =         0.71061
partial piv. (indirect)   - |A*x-b| =     1.11022e-16 - rel.error [%] =         0.71061
scaled partial piv.       - |A*x-b| =     4.71028e-16 - rel.error [%] =         3.27588

n=13, cond(hilb(13))=6.37569e+17

no piv.                   - |A*x-b| =     2.22045e-16 - rel.error [%] =         334.283
partial piv. (swap)       - |A*x-b| =     1.57009e-16 - rel.error [%] =         101.336
partial piv. (indirect)   - |A*x-b| =     4.71028e-16 - rel.error [%] =         101.336
scaled partial piv.       - |A*x-b| =     4.71028e-16 - rel.error [%] =         264.157