function y = myconv(x1,x2)
x3=fliplr(x2); len1=length(x1); len2=length(x2); len3=len1+len2-1; len4=len2-1; y=zeros(1,len3);
for q=1:1:len3 q1=q+len4; n2=zeros(1,len3); n2([q:q1])=x3; len5=length(n2); n1=zeros(1,len5); n1([len2:len3])=x1; c1=n1.*n2; c2=0; for i=1:1:len3 c2=c2+c1([i]); end sum=c2 y([q])=sum end
t=length(y); t2=(t-1)/2; v=-t2:1:t2; subplot(3,1,1),stem(x1,'g*'); subplot(3,1,2),stem(x2,'b*'); subplot(3,1,3),stem(v,y,'r*');
