e13p1.m

% Example 13.1  discrete system simulation   PH3

format compact

Gpnum = [1];  Gpden = [1 1 0];
[Ac,Bc,Cc,Dc] = tf2ss(Gpnum,Gpden);
[A,B] = c2d(Ac,Bc,1);
[Gznum,Gzden] = ss2tf(A,B,Cc,Dc)
pause
Tznum = Gznum,  Tzden = Gznum + Gzden
pause
[A,B,C,D] = tf2ss(Tznum,Tzden);
  x = [0;0];
  iter = 12
    itsol