%-----------------------------------
%求解解析解
%-----------------------------------
syms s ;            %定義字符變量
A=[-21,19,-20;
    19,-21,20;
    40,-40,-40];
root=eig(A)         %求矩陣A的特征根，并判斷是否是病態(tài)系統(tǒng)
B=[s,0,0;
    0,s,0;
    0,0,s];
C=B-A;              %（SI-A）
D=inv(C);           %（SI-A）的逆
eAT=ilaplace(D);    %exp(AT)=L((SI-A)-1)
X=[0;0;0];
X= eAT * [1;0;-1]   %解析解
t=0:0.01:2;         %繪圖采樣點間隔和區(qū)間長度
ezplot(X(1),t)      
%ezplot(X(2),t)
%ezplot(X(3),t)