L = 87;
hh = dtmfdesign([697 770 852 941 1209 1336 1477 1633], L, 8000);
ww = linspace(0,pi,500);
FRESP = freqz(hh(:, 6), 1, ww);
plot(ww, abs(FRESP)), grid
colors = ['b', 'g', 'r', 'c', 'm', 'r', 'b', 'g'];
hold on
for j = 1 : 8
    FRESP = freqz(hh(:, j), 1, ww);
    plot(ww, abs(FRESP), colors(j)), grid
    stem(2*pi*[697 770 852 941 1209 1336 1477 1633]/8000, ones(1,8), 'r')
xlabel('Scaled \omega=0 to \pi','FontSize',16,'FontWeight','bold')
ylabel('Magnitude','FontSize',16,'FontWeight','bold')
title(['FREQUENCY RESPONSES  (L=' num2str(L) ')'],'FontSize',16,'FontWeight','bold')
end



%hold off

%dtmfscore(hh(:,2),hh(:,5))