A1=A*Bet/L;
R=rcond(A1);
P=poly(A1);
w=polyroots(P);
r_w=real(w);