x=var('x')
y=function('y',x)
deq1=x*diff(y,x)+4*y-8*x^4==0
r=desolve(deq1,[y,x],ics=(1,2))
show(r.simplify_full())