%rutina para graficar en tiempo discreto una 
%funcion coseno con periodicidad N0=24
k=[0:60];
A=1;
alpha=0.9;
tetha0=pi/12;
fk=A*cos(tetha0*k);
figure(4),stem(k,fk,'r')
hold on
hk=alpha.^k;
gk=hk.*fk;
figure(4),stem(k,hk,'b')
hold on
figure(4),stem(k,gk,'k')