rand('state',1),A=rand(100,100)-0.5;
t0=clock;[V,D]=eig(A);T_full=etime(clock,t0)%指令eig的運作時間。
options.tol=1e-8;				%為eigs設定計算精度。
options.disp=0;					%使中間迭代結果不顯示。
t0=clock;[v,d]=eigs(A,1,'lr',options);%計算最大實部特征值和特征向量。
T_part=etime(clock,t0)			%指令eigs的運作時間。
[Dmr,k]=max(real(diag(D)));	%在eig求得的全部特征值中找最大實部的那個。
d,D(1,1)

vk1=V(:,k);							%與d相同的特征向量應是V的第k列。
vk1=vk1/norm(vk1);v=v/norm(v);	%向量長度歸一。
V_err=acos(norm(vk1'*v))*180/pi	%求復數向量之間的夾角（度）。
D_err=abs(D(k,k)-d)/abs(d)	%求兩個特征值間的相對誤差。