diff --git a/matlab/cherrypick.m b/matlab/cherrypick.m
index 85bfe6af4e01de4f3b182c7e617c732b69b50c71..3224d4f45fd222737db649d1dde23a8ac9c65814 100644
--- a/matlab/cherrypick.m
+++ b/matlab/cherrypick.m
@@ -134,10 +134,22 @@ for i=1:length(eqtags)
     end
     if ~isempty(ispac)
         [rhs, growthneutralitycorrection] = write_expectations(eqtags{i}, ispac.name, 'pac');
-        lhs = sprintf('%s_PE', eqtags{i});
-        RHS = strrep(RHS, sprintf('pac_expectation(model_name = %s)', ispac.name), lhs);
-        if ~isempty(growthneutralitycorrection)
-            RHS = sprintf('%s + %s', RHS, growthneutralitycorrection);
+        if ~isempty(rhs)
+            lhs = sprintf('%s_PE', eqtags{i});
+            if isempty(growthneutralitycorrection)
+                RHS = strrep(RHS, sprintf('pac_expectation(model_name = %s)', ispac.name), lhs);
+            else
+                RHS = strrep(RHS, sprintf('pac_expectation(model_name = %s)', ispac.name), sprintf('%s+%s', lhs, growthneutralitycorrection));
+            end
+        else
+            % MCE version of the PAC equation.
+            [rhs, growthneutralitycorrection] = write_pac_mce_expectations(eqtags{i}, ispac.name);
+            lhs = sprintf('%s_Z', eqtags{i});
+            if isempty(growthneutralitycorrection)
+                RHS = strrep(RHS, sprintf('pac_expectation(model_name = %s)', ispac.name), lhs);
+            else
+                RHS = strrep(RHS, sprintf('pac_expectation(model_name = %s)', ispac.name), sprintf('%s+%s', lhs, growthneutralitycorrection));
+            end
         end
     end
     % Print equation for unrolled PAC/VAR-expectation and update