4.3BSD/usr/contrib/apl/lib/fn/convolve

r { a convolve b
r { b +.X ((I1)-IRb)O|((Rb),`1+Ra,b)Ra,(`1+Rb)R0