function S = exm2_f(A)
[m,n] = size(A);
if m ~= n
    error('輸入矩陣應是方陣');
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;