x = var('x')
y = function('y',x)

DE = x*diff(y,x) - x - y

desolve(DE, [y,x])