var('x')
y=function('y',x)
print desolve(x^2*diff(y,x)+2*y*x-sinh(5*x),[y,x])