%%rutina para obtener graficas de tiempo discreto a partir de una funcion
%%en k
D=5;
k=[0:60];
alpha=1.011;
hk=alpha.^k;
figure(1),stem(k,hk);
numz=[1 0];
denz=[1 -alpha];
Hz=tf(numz,denz,1);
figure(2),dimpulse(numz,denz,61)