diff --git a/matlab/cherrypick.m b/matlab/cherrypick.m
index 9e1727b88381fdb6c02d3f43822e39c2c9b83e46..834592e5ecbe91499cb7954cef4419b664758520 100644
--- a/matlab/cherrypick.m
+++ b/matlab/cherrypick.m
@@ -135,10 +135,22 @@ try
         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