laplace inverse[k/(p^2+k^2),p,t]*heaviside(t)