q=var('x')
w=function('y',x)
deq1=diff(y,x)+6*x^2*y-(e^-(2*x^3))/x^2==0
r=desolve(deq1,[y,x],ics=(1,0))
show(r.simplify_full())