InverseLaplaceTransform[1/(s^2+s+1),s,t]*heaviside(t)