function S=exm2_f(A)
[m,n]=size(A);
if m~=n
	error('輸入矩陣應(yīng)是方陣！');
end;
e=eig(A);
%檢查輸入矩陣的特征值是否各異
same=0;
for i=1:m-1
   for j=(i+1):m
      if e(j)==e(i)
         same=1;
      end
   end
end
% A可以對角化的條件是A具有互異特征值或者A為埃爾米特矩陣。
if any(any((A'-A)))&(same==1)
   error('矩陣無法對角化！');
end
[v,d]=eig(A);
S=v;