diff --git a/src/api/python.rs b/src/api/python.rs
index 3ddb78e3..ce287ed7 100644
--- a/src/api/python.rs
+++ b/src/api/python.rs
@@ -8743,7 +8743,7 @@ impl PythonMatrix {
         })
     }
 
-    /// Create a new row vector from a list of scalars.
+    /// Create a new column vector from a list of scalars.
     #[classmethod]
     pub fn vec(
         _cls: &PyType,
diff --git a/src/tensors/matrix.rs b/src/tensors/matrix.rs
index 91ec7f40..ab9bbfed 100644
--- a/src/tensors/matrix.rs
+++ b/src/tensors/matrix.rs
@@ -59,7 +59,7 @@ impl<F: Ring> Matrix<F> {
         m
     }
 
-    /// Create a new row vector from a list of scalars.
+    /// Create a new column vector from a list of scalars.
     pub fn new_vec(data: Vec<F::Element>, field: F) -> Matrix<F> {
         Matrix {
             nrows: data.len() as u32,
diff --git a/symbolica.pyi b/symbolica.pyi
index 987a20e7..8090fc3e 100644
--- a/symbolica.pyi
+++ b/symbolica.pyi
@@ -2755,7 +2755,7 @@ class Matrix:
 
     @classmethod
     def vec(cls, entries: Sequence[RationalPolynomial | Polynomial | Expression | int]) -> Matrix:
-        """Create a new row vector from a list of scalars."""
+        """Create a new column vector from a list of scalars."""
 
     @classmethod
     def from_linear(cls, nrows: int, ncols: int, entries: Sequence[RationalPolynomial | Polynomial | Expression | int]) -> Matrix: