% 《工程線性代數（MATLAB版）》第六章例題6.4程序ea604
%
A=[3,2,4;2,0,2;4,2,3];      % 輸入系數矩陣
f=poly(A),r=roots(f),       % 求特征多項式并求其根（特征值）
for k=1:3
B(k)=real([A-r(k)*eye(3)]);   % 把三個特征值代入方程，去掉產生的虛數誤差 
% B=rref(B,1e-12),    % 若因計算誤差造成B非奇異,無基礎解，加這一句
p=null(B(k),'r')               % 求其基礎解系
end