x = var('x')
y = function('y',x)
DE = diff(y, x) - x*e^(-x^2)
desolve(DE, y,x)