I just finished teaching a chapter on using Laplace transform to solve constant coefficient second order linear differential equations. I touted how amazing the method was because it incorporates the initial data from the start, works for strange forcing terms, and reduces the problem of solving an ODE to computing the unit impulse response $e(t)$ (by taking the easy inverse Laplace of the reciprocal of the characteristic function) and the convolution $e*g$, where $g$ is the forcing term--notice how we do not need to compute Laplace of $g$.
I know theoretically this is significant, and that for discontinuous and non-standard forcing terms this is one of the best methods. However, I was left with the feeling that if we cannot really calculate the convolution (closed form) then this is not as impressive after all! So:
How commonly is this method actually used in practice for solving ODE -- say by engineers?
Are there ways to compute the convolution for a considerably large collection of pairs of functions?
If we cannot find convolution in closed form, is this method used to produce numerical solutions, e.g., by estimating the integral in the definition of the convolution?
References will be appreciated (over heuristics)!