Add qr_solve for solving linear systems

[MitchellThompkins]

qr and qr_hessenberg both return QRMatrix but there is no back-substitution step to actually solve Ax = b. This adds a qr_solve function that takes a QRMatrix and RHS vector, computes (Q^T)b, then solves via back-substitution through R.