% Problem 8.22
disp('F-15 Aircraft model for subsonic conditions:')
A1=[-0.00819 -25.70839 0       -32.17095
   -0.00019 -1.27626  1         0
    0.00069  1.02176 -2.40523   0
    0        0        1         0]
B1=[-6.80939 -0.14968 -14.06111 0]'
C1=eye(4)
D1=[0 0 0 0]'

disp('F-15 Aircraft model for suoersonic conditions:')
A2=[-0.01172 -95.91071  0       -32.11294
    -0.00011  -1.87942  1         0
     0.00056  -3.61627 -3.44478   0
     0         0        1         0]
B2=[-25.40405 -0.22042 -53.4246 0]'
C2=eye(4)
D2=[0 0 0 0]'

[num1,den1]=ss2tf(A1,B1,C1,D1,1)
[num2,den2]=ss2tf(A2,B2,C2,D2,1)
for i=1: 1: 4,
  for j=1: 1: 5,
    if abs(num1(i,j)) < 1.e-10
       num1(i,j)=0;
    end
    if abs(num2(i,j)) < 1.e-10
       num2(i,j)=0;
    end
  end
end