A = hilb(15) B = ones(1,15) linsolve(A,B') Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 1.543404e-18. ans = 1.0e+08 * 0.0000 -0.0000 0.0007 -0.0097 0.0738 -0.3122 0.7102 -0.5709 -1.0953 3.3002 -2.9619 -0.2353 2.4214 -1.7331 0.4121 cond(A) ans = 2.4469e+17 rcond(A) ans = 1.5434e-18 eig(A) ans = -0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0004 0.0056 0.0572 0.4266 1.8459