%M551_2_16_12 - see text, sec 2.8, pp. 63-64 format long A=[.780 .563; .913 .659] A = 0.780000000000000 0.563000000000000 0.913000000000000 0.659000000000000 b=[.217;.254] b = 0.217000000000000 0.254000000000000 x=A\b x = 0.999999999981709 -0.999999999974659 b-A*x ans = 1.0e-016 * -0.832667268468867 0 xex=[1;-1] xex = 1 -1 b-A*xex ans = 1.0e-016 * -0.832667268468867 0 xstar=[-0.443;1] xstar = -0.443000000000000 1.000000000000000 b-A*xstar ans = 1.0e-003 * -0.459999999999933 -0.541000000000014 svd(A) ans = 1.480952058643204 0.000000675241305 cond(A) ans = 2.193218999650770e+006 diary off