de:=diff(x(t),t,t)+4*x(t)=f(t); f:=t->8*Dirac(t-2*Pi); ic:=x(0)=3,D(x)(0)=0; dsolve({de,ic},x(t)); convert(%,piecewise,t);