diff --git a/matlab/get_companion_matrix.m b/matlab/get_companion_matrix.m
index 7ae48f0f3..2a689faa5 100644
--- a/matlab/get_companion_matrix.m
+++ b/matlab/get_companion_matrix.m
@@ -76,8 +76,6 @@ switch auxiliary_model_type
     % Get the indices of trend and EC equations in the auxiliary model.
     target_eqnums_in_auxiliary_model = find(M_.trend_component.(auxiliary_model_name).targets);
     ecm_eqnums_in_auxiliary_model = find(~M_.trend_component.(auxiliary_model_name).targets);
-    % Get the indices of trend equations in model.
-    target_eqnums = M_.trend_component.(auxiliary_model_name).target_eqn;
     % REMARK It is assumed that the non trend equations are the error correction
     %        equations. We assume that the model can be cast in the following form:
     %
diff --git a/matlab/get_companion_matrix_legacy.m b/matlab/get_companion_matrix_legacy.m
index c1893c756..a7d6345e7 100644
--- a/matlab/get_companion_matrix_legacy.m
+++ b/matlab/get_companion_matrix_legacy.m
@@ -73,8 +73,6 @@ switch auxiliary_model_type
     % Get the indices of trend and EC equations in the auxiliary model.
     target_eqnums_in_auxiliary_model = find(M_.trend_component.(auxiliary_model_name).targets);
     ecm_eqnums_in_auxiliary_model = find(~M_.trend_component.(auxiliary_model_name).targets);
-    % Get the indices of trend equations in model.
-    target_eqnums = M_.trend_component.(auxiliary_model_name).target_eqn;
     % REMARK It is assumed that the non trend equations are the error correction
     %        equations. We assume that the model can be cast in the following form:
     %