diff --git a/matlab/perfect-foresight-models/sim1_linear.m b/matlab/perfect-foresight-models/sim1_linear.m
index 080a60a5dd75ce70d7422d022a2537e010b3c7bb..72a44ff9dd0b498aafa47a51bceceb12f4d7cc31 100644
--- a/matlab/perfect-foresight-models/sim1_linear.m
+++ b/matlab/perfect-foresight-models/sim1_linear.m
@@ -108,11 +108,12 @@ if verbose
 end
 
 dynamicmodel = str2func([M.fname,'_dynamic']);
+
 z = steadystate_y([ip; ic; in]);
+x = repmat(transpose(steadystate_x), 1+M.maximum_exo_lag+M.maximum_exo_lead, 1);
 
 % Evaluate the Jacobian of the dynamic model at the deterministic steady state.
-[d1,jacobian] = dynamicmodel(z, repmat(transpose(steadystate_x),options.periods+M.maximum_lag+M.maximum_lead,1),...
-    params, steadystate_y, M.maximum_lag+1);
+[d1, jacobian] = dynamicmodel(z, x, params, steadystate_y, M.maximum_lag+1);
 
 % Check that the dynamic model was evaluated at the steady state.
 if max(abs(d1))>1e-12