diff --git a/equationcollection/equationcollection.py b/equationcollection/equationcollection.py
index 05787b5f3e6766b8415136fe463118d352d654ee..86fb5700cab04af7e1cf7a2c022dff05be2160d3 100644
--- a/equationcollection/equationcollection.py
+++ b/equationcollection/equationcollection.py
@@ -66,8 +66,8 @@ class EquationCollection(object):
         newEquations = [fastSubs(eq, substitutionDict) for eq in self.mainEquations]
         if addSubstitutionsAsSubexpressions:
             newSubexpressions = [sp.Eq(b, a) for a, b in substitutionDict.items()] + newSubexpressions
-
-        return self.copy(newEquations, sortEquationsTopologically(newSubexpressions))
+            newSubexpressions = sortEquationsTopologically(newSubexpressions)
+        return self.copy(newEquations, newSubexpressions)
 
     def addSimplificationHint(self, key, value):
         """