function [val1,val2]=spect_rad(A)
L=extractL(A);
D=diag(diag(A));
I=eye(size(A));
val1=max(abs(eig(I-inv(D)*A)))
val2=max(abs(eig(I-inv(L)*A)))
