Your Code Goes Here
x=var('x')
y=function('y',x)
deq1=x^2*diff(y,x)+2*x*y-sinh(5*x)==0
r=desolve(deq1,[y,x],y)
show(r.simplify_full())