%exm09634_0.m		例 9.6.3.4_1 的主運行程序

%求坐標原點處的線性化模型
[A,B,C,D]=linmod('exm09634_1');A  

%求[ 1 , 0.5]坐標處的線性化模型
[A1,B1,C1,D1]=linmod('exm09634_1',[1,0.5]);A1  

%由特征根看模型穩定性
eA=eig(A)';eA1=eig(A1)';
disp('原點處線性化模型的特征根,表明“穩定”。')
disp(eA)
disp(' ')
disp('[ 1 , 0.5]坐標處線性化模型的特征根，表明“不穩定”。')
disp(eA1)