diff --git a/spec/API_specification/linear_algebra_functions.md b/spec/API_specification/linear_algebra_functions.md
index 209a2f787..109d23666 100644
--- a/spec/API_specification/linear_algebra_functions.md
+++ b/spec/API_specification/linear_algebra_functions.md
@@ -259,10 +259,36 @@ Computes the outer product of two vectors `x1` and `x2`.
 
 TODO
 
-(function-qr)=
-### qr()
+(function-linalg-qr)=
+### linalg.qr(x, /, *, mode='reduced')
 
-TODO
+Computes the qr factorization of a matrix (or a stack of matrices), where `q` is an orthonormal matrix (or a stack of matrices) and `r` is an upper-triangular matrix (or a stack of matrices).
+
+#### Parameters
+
+-   **x**: _<array>_
+
+    -   input array having shape `(..., M, N)` and whose innermost two dimensions form `MxN` matrices. Should have a floating-point data type.
+
+-   **mode**: _str_
+
+    -   factorization mode. Should be one of the following modes:
+
+        -   `'reduced'`: compute only the leading `K` columns of `q`, such that `q` and `r` have dimensions `(..., M, K)` and `(..., K, N)`, respectively, and where `K = min(M, N)`.
+        -   `'complete'`: compute `q` and `r` with dimensions `(..., M, M)` and `(..., M, N)`, respectively.
+
+        Default: `'reduced'`.
+
+#### Returns
+
+-   **out**: _Tuple\[ <array>, <array> ]_
+
+    -   a namedtuple `(q, r)` whose
+
+        -   first element must be an array whose shape depends on the value of `mode` and contain orthonormal matrices. If `mode` is `'complete'`, the array must have shape `(..., M, M)`. If `mode` is `'reduced'`, the array must have shape `(..., M, K)`, where `K = min(M, N)`. The first `x.ndim-2` dimensions must have the same size as those of the input `x`.
+        -   second element must be an array whose shape depends on the value of `mode` and contain upper-triangular matrices. If `mode` is `'complete'`, the array must have shape `(..., M, M)`. If `mode` is `'reduced'`, the array must have shape `(..., K, N)`, where `K = min(M, N)`. The first `x.ndim-2` dimensions must have the same size as those of the input `x`.
+
+        Each returned array must have a floating-point data type determined by {ref}`type-promotion`.
 
 (function-slogdet)=
 ### slogdet()