diff --git a/libs/geom/include/psemek/geom/gauss.hpp b/libs/geom/include/psemek/geom/gauss.hpp
index 479d5de2..417cdced 100644
--- a/libs/geom/include/psemek/geom/gauss.hpp
+++ b/libs/geom/include/psemek/geom/gauss.hpp
@@ -197,6 +197,14 @@ namespace psemek::geom
 		return true;
 	}
 
+	template <typename T, std::size_t N>
+	std::optional<vector<T, N>> solve(matrix<T, N, N> const & m, vector<T, N> v)
+	{
+		if (!gauss(m, v))
+			return std::nullopt;
+		return v;
+	}
+
 	template <typename T, std::size_t N>
 	std::optional<matrix<T, N, N>> inverse(matrix<T, N, N> m)
 	{