public static double SpectralRadius(MMatrix A)
{
Vector Eigenvalues = new Vector();
if (!A.IsSquared())
throw new MMatrixException("Matrix must be squared.");
Eigenvalues = A.QRIterationBasic(GlobalMath.MAXITER).DiagVector();
Eigenvalues.Sort(true);
return Math.Abs(Eigenvalues[0]);
}