x=var('x')#1
y=function('y',x)
DE=diff(y,x,x)-4*diff(y,x)+4*y-exp(2*x)/x
desolve(DE,[y,x])