diff --git a/matlab/mr_gstep.m b/matlab/mr_gstep.m
index db25d43a4..381b21fa0 100644
--- a/matlab/mr_gstep.m
+++ b/matlab/mr_gstep.m
@@ -57,7 +57,7 @@ while i<n
         gg(i)=(f1(i)'-f_1(i)')./(2.*h1(i));
         hh(i) = 1/max(1.e-9,abs( (f1(i)+f_1(i)-2*f0)./(h1(i)*h1(i)) ));
         if gg(i)*(hh(i)*gg(i))/2 > htol
-            [f0 x fc retcode] = csminit(func0,x,f0,gg,0,diag(hh),DynareDataset,DynareOptions,Model,EstimatedParameters,BayesInfo,DynareResults);
+            [f0 x fc retcode] = csminit1(func0,x,f0,gg,0,diag(hh),DynareDataset,DynareOptions,Model,EstimatedParameters,BayesInfo,DynareResults);
             ig(i)=1;
             fprintf(['Done for param %s = %8.4f\n'],BayesInfo.name{i},x(i))
         end