diff --git a/methods/momentbased.py b/methods/momentbased.py
index a402869e7aac9b40688960ae2b472d2d83d74960..84531f98ebbcedf036b5f5d1ed86cd9fc93cedb0 100644
--- a/methods/momentbased.py
+++ b/methods/momentbased.py
@@ -104,7 +104,7 @@ class MomentBasedLbMethod(AbstractLbMethod):
                 raise NotImplementedError("Shear moments seem to be not relaxed separately - "
                                           "Can not determine their relaxation rate automatically")
 
-    def getEquilibrium(self, conservedQuantityEquations=None, includeForceTerms=True):
+    def getEquilibrium(self, conservedQuantityEquations=None, includeForceTerms=False):
         D = sp.eye(len(self.relaxationRates))
         return self._getCollisionRuleWithRelaxationMatrix(D, conservedQuantityEquations=conservedQuantityEquations,
                                                           includeForceTerms=includeForceTerms)
diff --git a/moments.py b/moments.py
index e7a0eff29075bf63933b7fd9234e285ee5fbb188..3a7c5579ad0c2705ce61352a622ef63ca7455731 100644
--- a/moments.py
+++ b/moments.py
@@ -184,8 +184,8 @@ def polynomialToExponentRepresentation(polynomial, dim=3):
 
     Example:
         >>> x , y, z = MOMENT_SYMBOLS
-        >>> polynomialToExponentRepresentation(1 + (42 * x**2 * y**2 * z) )
-        [(1, (0, 0, 0)), (42, (2, 2, 1))]
+        >>> set(polynomialToExponentRepresentation(1 + (42 * x**2 * y**2 * z) )) == {(42, (2, 2, 1)), (1, (0, 0, 0))}
+        True
     """
     assert dim <= 3
     x, y, z = MOMENT_SYMBOLS