Example,
Compute the qr factorization of a matrix.
Factor the matrix `a` as `qr`, where `q` is orthonormal
(:math:`dot( q_{:,i}, q_{:,j}) = \delta_{ij}`, the Kronecker delta) and
`r` is upper-triangular.
Parameters
----------
a : array_like, shape (M, N)
Matrix to be factored.
mode : {'full', 'r', 'economic'}
Specifies the information to be returned. 'full' is the default.
mode='r' returns a "true" `r`, while 'economic' returns a "polluted"
`r` (albeit slightly faster; see Returns below).
Returns
-------
* If mode = 'full':
* q : ndarray of float or complex, shape (M, K)
* r : ndarray of float or complex, shape (K, N)
Size K = min(M, N)
* If mode = 'r':
* r : ndarray of float or complex, shape (K, N)
* If mode = 'economic':
* a2 : ndarray of float or complex, shape (M, N)
The diagonal and the upper triangle of a2 contains r,
while the rest of the matrix is undefined.