argout==0;
0127     
0128     <span class="comment">% no output, just plot</span>
0129     disp([<span class="string">'cluster1: '</span>,num2str(100*numel(find(A==1))/numel(A)), <span class="string">'%'</span>]);
0130 
0131     figure(100); clf
0132     subplot 311;
0133     plot(x); hold on
0134     xx=x;
0135     xx(IDX{1},:)=nan;
0136     plot(xx,<span class="string">'k'</span>);
0137     axis tight
0138     title(<span class="string">'black: cluster B'</span>);
0139     
0140     z1=<a href="nt_mmat.ht