diff --git a/matlab/perfect-foresight-models/sim1.m b/matlab/perfect-foresight-models/sim1.m
index 7f9fc84c6e4ac1d2dad90b9daafcb18f4f989837..8f0f65a90ed8aa338664a3d0818e4c3d8c24bce6 100644
--- a/matlab/perfect-foresight-models/sim1.m
+++ b/matlab/perfect-foresight-models/sim1.m
@@ -68,8 +68,9 @@ for iter = 1:options.simul.maxit
     % A is the stacked Jacobian with period x equations alongs the rows and
     % periods times variables (in declaration order) along the columns
     if options.debug && iter==1
-        row=find(all(A==0,2));
-        column=find(all(A==0,1));
+        [row,col]=find(A);
+        row=setdiff(1:periods*ny,row);
+        column=setdiff(1:periods*ny,col);
         if ~isempty(row) || ~isempty(column)
             fprintf('The stacked Jacobian is singular. The problem derives from:\n')
             if ~isempty(row)