lass="comment">% +1 and one -1. RHO is computed from the weights W as follows:</span>
<span class="comment">%    RHO = max(abs(eig( eye(n,n) - (1/n)*ones(n,n) - A*W*A' ))).</span>
<span class="comment">%</span>
<span class="comment">% The M.-H. weight on an edge is one over the maximum of the degrees of the</span>
<span class="comment">% adjacent nodes.</span>
<span class="comment">%</span>
<span class="comment">% For more details, see the references:</span>
<span class="comment">% "Fast linear iterations f