t(x,topcs);
0052     x=nt_normcol(x);
0053 <span class="keyword">end</span>
0054 
0055 <span class="comment">% form cross products</span>
0056 xx=nt_xprod(x);
0057 
0058 <span class="comment">% DSS finds the most reproducible quadratic form</span>
0059 keep1=[]; keep2=0;<span class="comment">%10^-12;</span>
0060 [toquads,pwr0,pwr1]=nt_dss1(xx,[],keep1,keep2);
0061 
0062 <span class="comment">% illustrate</span>
0063 figure(1); clf
0064 plot(pwr1./pwr0,<span class="string">'.-'</span>)
0065 z=nt_mmat(xx,toqu