diff --git a/lbmpy/methods/creationfunctions.py b/lbmpy/methods/creationfunctions.py
index 93174eb7ed3bbfee0e4b92b048b680ed1da65e46..131ad9b493b1ca209557e6de67a4a5b06d5a3274 100644
--- a/lbmpy/methods/creationfunctions.py
+++ b/lbmpy/methods/creationfunctions.py
@@ -396,7 +396,7 @@ def create_mrt_orthogonal(stencil, relaxation_rate_getter=None, maxwellian_momen
         nested_moments = [
             [one, x, y, z],  # [0, 3, 5, 7]
             [sq - 1],  # [1]
-            [3 * sq ** 2 - 6 * sq + 1],  # [2]
+            [3 * sq ** 2 - 9 * sq + 4],  # [2]
             [(3 * sq - 5) * x, (3 * sq - 5) * y, (3 * sq - 5) * z],  # [4, 6, 8]
             [3 * x ** 2 - sq, y ** 2 - z ** 2, x * y, y * z, x * z],  # [9, 10, 11, 12, 13]
             [x * y * z]
diff --git a/lbmpy_tests/test_momentbased_methods_equilibrium.py b/lbmpy_tests/test_momentbased_methods_equilibrium.py
index f2698f3c8910f50493fba4c70cd23d549d85c024..5f3381d37f3ac88bacf0c4c9a5515ec959255081 100644
--- a/lbmpy_tests/test_momentbased_methods_equilibrium.py
+++ b/lbmpy_tests/test_momentbased_methods_equilibrium.py
@@ -72,6 +72,9 @@ def test_mrt_orthogonal():
     m = create_mrt_orthogonal(get_stencil("D2Q9"), maxwellian_moments=True)
     assert m.is_orthogonal
 
+    m = create_mrt_orthogonal(get_stencil("D3Q15"), maxwellian_moments=True)
+    assert m.is_weighted_orthogonal
+
     m = create_mrt_orthogonal(get_stencil("D3Q19"), maxwellian_moments=True)
     assert m.is_weighted_orthogonal