@(p,q)ifft(fft([p q*0]).*fft([q p*0]))(1:end-1) ans([1 1], [1 4 6 4 1])