n=[-10:10];
Dn=0.504./(1+j*4*n);
subplot(211);
stem(n,abs(Dn));
xlabel('n');
ylabel('|D_n|');
subplot(212);
stem(n,angle(Dn));
xlabel('n');
ylabel('\angle D_n');