%tiempo continuo
%num=input('num=');
%den=input('den=');
%H1=tf(num,den)
%figure(1),impulse(H1)
%segundo orden
%num=[1 0];
%den=[1 0 16];
%H2=tf(num,den)
%figure(2),impulse(H2)
%tiempo discreto
Ts=1;
numz=input('numz=');
denz=input('denz=');
Hz=tf(numz,denz,Ts)
hk=dimpulse(numz,denz,41); 
k=0:40;
figure(3),stem(k,hk);
yk=dstep(numz,denz,13);
l=0:12;
figure(4),stem(l,yk)