x = var('x')   
y = function('y',x)   
sol=desolve (diff(y,x) + 6*x^2*y - (exp(-2*x^3))/x^2, [y,x], ics=(1,0))
sol

---- /!\ '''Edit conflict - other version:''' ----

/!\ Edit conflict - your version:

/!\ End of edit conflict