function y=convolution1(x,h)

M=length(x);
N=length(h);
temp=zeros(1,M+N-1);
for n=1:1:M+N-1
    for m=1:1:M
        k=n-m+1;
        if(k>=1&&k<=N)
            temp(n)=temp(n)+x(m)*h(n-m+1);
        end
    end
end

y=temp;