diff --git a/matlab/compute_Pinf_Pstar.m b/matlab/compute_Pinf_Pstar.m
index aa56b39e69dd1a4fd29e47676d061c5724187c94..150acac1877f43df6aae9ea0e2b1f8fafcdf3796 100644
--- a/matlab/compute_Pinf_Pstar.m
+++ b/matlab/compute_Pinf_Pstar.m
@@ -169,6 +169,7 @@ if np0
     QTinf(1:np0,np0+1:np0+nk) = STinf0;
     QTinf([indx0(:); indx(:)],:) = QTinf;
     STinf1 = [zeros(np0+np,np0) [STinf0; eye(nk); zeros(np-nk,nk)] zeros(np0+np,np-nk)];
+    mf = ismember([indx0(:); indx(:)],mf);
     for k = 1:nk
         if norm(QTinf(mf,:)*ST([indx0(:); indx(:)],k+np0)) < 1e-8
             Pinf(k+np0,k+np0) = 0;