This method will decompose a matrix into an upper triangular and an orthogonal matrix.
A = [4 12 -16
12 37 -43
-16 -43 98];
R = qr(A);
This will return the upper triangular matrix while the following will return both matrices.
[Q,R] = qr(A);
The following plot will display the runtime of qr
dependent of the square root of elements of the matrix.