heaviside(-t)*inverselaplacetransform(1/((1+p^2)