A=[-10 -35 -50 -24;1 0 0 0;0 1 0 0;0 0 1 0];
B=[1;0;0;0];C=[1 7 24 24];D=[0];
[z,p,k]=ss2zp(A,B,C,D,1);                       %求零極點增益模型
Flagz=0;                                     %設標志
n=length(A);
for i=1:n 
    if real(p(i))>0
        Flagz=1;
    end
end
disp('系統的零極點模型為');z,p,k
if Flagz==1
    disp('系統不穩定');
else
    disp('系統是穩定的');
end
step(A,B,C,D);                                %求系統狀態方程的階躍響應