Givens QR factorization:

matrix A:
    0.683     0.713     0.548 
    0.672     0.953     0.936 
    0.225     0.599     0.396 
    0.645     0.175     0.398 

matrix Q^T*Q:
    1.000     0.000     0.000 
    0.000     1.000     0.000 
    0.000     0.000     1.000 

matrix R:
    1.177     1.169     1.146 
    0.000     0.664     0.376 
    0.000     0.000     0.188 

matrix Q*R:
    0.683     0.713     0.548 
    0.672     0.953     0.936 
    0.225     0.599     0.396 
    0.645     0.175     0.398 
b=        1        2        3
Givens QR solution of a square system:
matrix A:
    0.057     0.717     0.584 
    0.721     0.913     0.744 
    0.119     0.286     0.763 
          vector b:        1        2        3
solution x to Ax=b:     1.12    -2.53     4.71
         vector Ax:        1        2        3
A^(-1):
   -1.965     1.542     0.002 
    1.876     0.106    -1.540 
   -0.398    -0.280     1.890 
A A^(-1):
    1.000     0.000     0.000 
    0.000     1.000     0.000 
    0.000     0.000     1.000 
A^(-1) A:
    1.000     0.000     0.000 
    0.000     1.000     0.000 
    0.000     0.000     1.000