function [u]=coj(s,r)
m=length(s);
n=length(r);
p=m;
if n > m,
p=n;
end
u=zeros(p,1);
for i=1:m,
for j=1:n,
k=i-j;
num=abs(k)+1;
u(num)=u(num)+s(i)*r(j);
end;
end;
end
