integrate exp(-p*t)*((2* (t-3)^2+3) theta(t-3)) for t=0 to inf