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]); }