QR decomposition
// formula_24 is the precision of the QR decomposition, e.g., 1e-6.