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

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