A =
     1.000     1.000     1.000
     1.000     3.000     1.000
     1.000     1.000     5.000
 b =
     1.000
     4.000
     1.000
 ---------------------
 Q =
     0.577    -0.408    -0.707
     0.577     0.816     0.000
     0.577    -0.408     0.707
 R =
     1.732     2.887     4.041
     0.000     1.633    -1.633
     0.000     0.000     2.828
 QR =
     1.000     1.000     1.000
     1.000     3.000     1.000
     1.000     1.000     5.000
 Q^T*Q =
     1.000     0.000     0.000
     0.000     1.000     0.000
     0.000     0.000     1.000
 Ax=b => x =
    -4.919
     2.803
     0.354
 ---------------------
abs(det(A)) =    8.000
 ---------------------
 A^-1 =
     1.750    -0.500    -0.250
    -0.500     0.500     0.000
    -0.250     0.000     0.250
 A^-1*A =
     1.000     0.000    -0.000
    -0.000     1.000     0.000
    -0.000     0.000     1.000
 Unit(4) =
     1.000     0.000     0.000     0.000
     0.000     1.000     0.000     0.000
     0.000     0.000     1.000     0.000
     0.000     0.000     0.000     1.000