%===================================================
% program 7_1
% Solving a constant coefficient difference equation
% The input sequence must be a closed representation
%===================================================
N=input('Desired impulse response length = ');
n=0:N;
p=input('Type in the vector p = ');
d=input('Type in the vector d = ');
x=input('Type in the input sequence x(n) = ');%[1 zeros(1,N-1)];
y=filter(p,d,x);
stem(n,y)
title('The output sequence y(n)')
xlabel('Time index n');ylabel('Amplitude')