n=0:20;
x=sin(0.2*pi*n);
b=[1 0 -0.25];
a=[1 -0.5];
y=filter(b,a,x);
stem(n,y);
xlabel('n')
ylabel('y(n)');
title ('output using filter()');
grid on;
disp(y);
disp(x);