diff --git a/DynamicModel.cc b/DynamicModel.cc
index eb4e6688b7b6785907a0446e349e231a538218fb..3b274be70cc0ec487cf669ea979c9b5018d5ecd9 100644
--- a/DynamicModel.cc
+++ b/DynamicModel.cc
@@ -1195,7 +1195,7 @@ DynamicModel::writeDynamicMFile(const string &dynamic_basename) const
         cerr << "Error: Can't open file " << filename << " for writing" << endl;
         exit(EXIT_FAILURE);
       }
-    mDynamicModelFile << "function [residual, G1, g2, g3] = " << dynamic_basename << "(y, x, params, it_)" << endl
+    mDynamicModelFile << "function [residual, g1, g2, g3] = " << dynamic_basename << "(y, x, params, it_)" << endl
     << "%" << endl
     << "% Status : Computes dynamic model for Dynare" << endl
     << "%" << endl
@@ -1207,8 +1207,6 @@ DynamicModel::writeDynamicMFile(const string &dynamic_basename) const
 
     writeDynamicModel(mDynamicModelFile, false);
 
-    mDynamicModelFile << "G1 = g1;" << endl;
-
     mDynamicModelFile.close();
   }
 
@@ -1917,9 +1915,7 @@ DynamicModel::writeDynamicModel(ostream &DynamicOutput, bool use_dll) const
         << model_output.str()
         // Writing initialization instruction for matrix g1
         << "if nargout >= 2," << endl
-	<< "  if isempty(g1)" << endl// Test if g1 isn't initialized.
-        << "    g1 = zeros(" << nrows << ", " << dynJacobianColsNbr << ");" << endl
-	<< "  end" << endl
+        << "  g1 = zeros(" << nrows << ", " << dynJacobianColsNbr << ");" << endl
         << endl
         << "%" << endl
         << "% Jacobian matrix" << endl