Gpnum = [0 0 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);
sys = ss(A,B,C,D,1)
[c,t] = step(sys,0:1:10)