convolution(y_,x_):=inverselaplace(laplace(x)*laplace(y));
convolution(a_+b_,x_):=convolution(a,x)+convolution(b,x);

convolution(a_):=int(replace(a,x,y)*(x-y),y,0,x);