x = var('x')
y = function('y',x)
DE = x*diff(y,x,2) + diff(y,x)
desolve(DE, [y,x])