answered 4 years ago
The code
x = var('x') y = function("y")(x) eq = x*e**y == x - y dy = diff(y,x) sol = solve(diff(eq), dy) show(dy.subs(sol))
yields −ey(x)−1xey(x)+1 as expected