diff --git a/matrix.cpp b/matrix.cpp
index dcdb151d..28644528 100644
--- a/matrix.cpp
+++ b/matrix.cpp
@@ -130,7 +130,22 @@ util::matrix<S,T>::inverse (void) const
 }
 
 
-//-----------------------------------------------------------------------------
+///////////////////////////////////////////////////////////////////////////////
+template <size_t S, typename T>
+matrix<S,T>
+util::transposed (const matrix<S,T> &m)
+{
+    util::matrix<S,T> res;
+
+    for (size_t y = 0; y < S; ++y)
+        for (size_t x = 0; x < S; ++x)
+            res[y][x] = m[x][y];
+
+    return res;
+}
+
+
+///////////////////////////////////////////////////////////////////////////////
 template <size_t S, typename T>
 matrix<S,T>
 matrix<S,T>::operator* (const matrix<S,T> &rhs) const
diff --git a/matrix.hpp b/matrix.hpp
index 70fb5987..2f1903a6 100644
--- a/matrix.hpp
+++ b/matrix.hpp
@@ -137,6 +137,12 @@ namespace util {
     matrix<S,T>
     inverse (const matrix<S,T>&);
 
+    template <size_t S, typename T>
+    matrix<S,T>
+    transposed (const matrix<S,T>&);
+
+    ///////////////////////////////////////////////////////////////////////////
+
     template <size_t S, typename T>
     matrix<S,T>
     abs (const matrix<S,T>&);