Krylov.m

resid =

   8.7611e-01   2.3333e-01   4.1780e-02   4.6472e-03

xm =

  -3.7271e+00
   4.4943e+01
  -8.5429e+01
  -3.5649e+01
   9.8221e+01

A*xm
   9.9997e-01
   1.0005e+00
   9.9764e-01
   1.0036e+00
   9.9825e-01
bicgstab converged at iteration 9 to a solution with relative residual 7.235942e-14
xm =

   5.0000e+00
  -1.2000e+02
   6.3000e+02
  -1.1200e+03
   6.3000e+02

res =    1.9844e-13
   1.0000e+00
   1.0000e+00
   1.0000e+00
   1.0000e+00
   1.0000e+00
ans =

   5.0000e+00   5.0000e+00

V =

   7.1992e-01   9.3049e-01   1.6806e-01   1.9583e-03   5.7367e-07
   4.5509e-01   2.5489e-01  -7.9377e-01  -1.7627e-01  -1.8479e-03
   3.5097e-01   2.9936e-04  -4.5105e-01   7.3814e-01   1.6733e-01
   2.9318e-01  -1.3729e-01  -1.6384e-02   4.0248e-01  -8.6039e-01
   2.5588e-01  -2.2443e-01   3.7144e-01  -5.1194e-01   4.8139e-01

D =

Diagonal Matrix

   9.9594e+00            0            0            0            0
            0  -1.4755e+00            0            0            0
            0            0   4.4657e-02            0            0
            0            0            0  -8.5823e-05            0
            0            0            0            0   2.8800e-09

ans =

   4.4409e-15  -1.1102e-15  -1.0131e-15   1.8231e-15  -1.2496e-15
  -3.5527e-15   2.7756e-16  -1.0339e-15   2.0315e-15  -1.9802e-15
  -4.4409e-16  -5.3592e-16   0.0000e+00   9.8460e-16  -7.6647e-16
  -1.3323e-15  -3.0531e-16  -5.3950e-16   1.0823e-15  -8.9578e-16
  -1.3323e-15  -3.8858e-16  -1.6306e-16   9.5070e-16  -7.8174e-16

ans =

   4.4409e-15
  -3.5527e-15
  -4.4409e-16
  -1.3323e-15
  -1.3323e-15