% 《工程線性代數（MATLAB版）》第六章例題6.2.1程序ea621a
%
A=[3,2,4;2,0,2;4,2,3];      % 輸入系數矩陣
f=poly(A),r=roots(f),       % 求特征多項式并求其根（特征值）
B1=real([A-r(1)*eye(3)]);   % 把三個特征值代入方程，去掉產生的虛數誤差 
p1=null(B1,'r')               % 求其基礎解系
B2=real([A-r(2)*eye(3)]);   % 把三個特征值代入方程，去掉產生的虛數誤差 
p2=null(B2,'r')               % 求其基礎解系
B3=real([A-r(3)*eye(3)]);   % 把三個特征值代入方程，去掉產生的虛數誤差 
p3=null(B3,'r')               % 求其基礎解系