solve y'+2y==DiracDelta[t] + HeavisideTheta[t]