e13p12.m

% Example 13.12   PH3

format compact

num = [0 0 0 4];
den = [1 3 2 0];
[A,B,C,D] = tf2ss(num,den);
[Ad,Bd] = c2d(A,B,0.1);
[Gnum,Gden] = ss2tf(Ad,Bd,C,0);
w1 = i*1.704;
z1 = (1 + (0.1/2)*w1)/