pkg load symbolic

n = 0:1:127 % sets the limits of the functions

% define functions
x= (3*sin(n*pi/8)+3).*(heaviside(n) - heaviside(n-63))
z= 7*cos(n*pi/16).*cos(n*pi).*(heaviside(n) - heaviside(n-63))
w=x+z
h=(1/2)*sinc(0.5*(n-31)).*(heaviside(n) - heaviside(n-63))

y=conv(x,h,'same')



% plot the functions
figure
plot (n,x)

figure
plot (n,z)

figure
plot (n,w)

figure
plot (n,h)

figure
plot (n,y)