x=var('x')#6.  
y=function('y',x)
DE=diff(y,x,x)-4*diff(y,x)+5*y-exp(2*x)*csc(x)
desolve(DE,[y,x])