Since I is the convolution of f(x) = x^2 and g(x) = cos(3x), you probably figured out that
L[I](s) = L[f](s) L[g](s) = 2/(s^2(s^2+9)). N
ow we want to take the inverse laplace transform, but first we should change this function into a slightly nicer form. The common way to do this is using partial fractions:
we know that 2/(s^2(s^2+9)) = A/s+B/s^2+(Cs+D)/(s^2+9).
It's not to bad to apply the inverse laplace transform to each of these terms, and L^-1 is linear, so once we solve for A,B,C,D, we'll pretty much be done. In the end you should get the same answer as Yefim got above using integration by parts :)