function lam2 = lambda2(A) ae = sort(abs(eig(A))); lam2 = ae(end-1); return