n>;
            ind=0;
            <span class="keyword">for</span> j=1:size(Y,1),
                ind = max(min(find(Y(j,:)&gt;=0.5)),ind);
            <span class="keyword">end</span>
            tdom   = max(eig(inv(GQ)*C));
            elmore = max(sum((inv(GQ)*C)'));
            tthres = T(ind);
            plot( tdom   * [1;1], [0;1], <span class="string">'--'</span>, <span class="keyword">...</span>
                  elmore * [1;1], [0;1], <span class="string">'--'</span>, <span class="keyword">...</