clc
clear

global numc denc h numd dend

a1 = 1;
a2 = 2;
a3 = 3;
h = 0.2;

numc = [1,a1];
numc = (a2*a3/a1)*[1,a1];
denc = conv([1,a2],[1,a3]);

G = tf(numc,denc)

Gz = c2d(G,h)


%numd = [0.4929, -0.4459];

%dend = [1, -1.56,0.6065];

[numd dend] = tfdata(Gz,'v');
