z = 2; t1 = 1; t2 = 2; ts = 1;
G = tf(z, conv([t1 1],[t2 1]));
zinteg = tf([1 0],[1 -1], ts);
