x=var('x')
y=function('y',x)
de=diff(y,x,2)+2*diff(y,x)+5*y-1.25*exp(0.5*x)-40*cos(4*x)-55*sin(4*x)
desolve(de,[y,x],ics=[0,0.2,60.1])