x=var('x')#3.  
y=function('y',x)
DE=diff(y,x,x)+2*diff(y,x)+y-exp(-x)*cos(x)
desolve(DE,[y,x])