=[1:k-1,k+1:n];
0028     [topcs,eigenvalues]=<a href="nt_pcarot.html" class="code" title="function [topcs,eigenvalues]=nt_pcarot(cov,N)">nt_pcarot</a>(cc(idx,idx)); <span class="comment">% PCA to orthogonalize the other channels</span>
0029     topcs=topcs(:,find(eigenvalues/max(eigenvalues)&gt;thresh));
0030     b=x(:,idx)*topcs(:,:);
0031     b=<a href="nt_normcol.html" class="code" title="function [y,norm]=nt_normcol(x,w)">nt_normcol</a>(b); <span class="comment">% this could be optimized</span>
0032    