inverse laplace transform(1/[(p+1)(p-1)^3] + (p²-3p+1)/(p-1)^3)*heaviside(t)