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);