clc;
s = tf('s');
[G,Gd,A,B,C,D,E]=computeTransferFunctions;

%[z11, p11, k11] = zpkdata(G(1,1), 'v')
%[z12, p12, k12] = zpkdata(G(1,2), 'v')
[z21, p21, k21] = zpkdata(G(2,1), 'v')
%[z22, p22, k22] = zpkdata(G(2,2), 'v')
%[zd11, pd11, kd11] = zpkdata(Gd(1,1), 'v')
%[zd12, pd12, kd12] = zpkdata(Gd(1,2), 'v')
%[zd21, pd21, kd21] = zpkdata(Gd(2,1), 'v')
[zd22, pd22, kd22] = zpkdata(Gd(2,2), 'v')


%G21 = (0.2375*s - 0.059375)/(s*(s+9.75));
%G21=minreal(G21);
%G21=set(zpk(G21), 'DisplayFormat','time constant');
%[z, p, k] = zpkdata(G21, 'v')

