Added the first implementation of QR factorization.

An extension of previous factorization (triangular with householder transformation).
There seems to a large round error.
Requires performance evaluation.