| 1 | initial version |
You don't consider the right equation to solve:
y=function('y')(x)
dy = diff(y)
derivative = solve(diff(y==exp(x*y)), dy)
derivative[0].subs(exp(x*y)==y)
This yields
diff(y(x), x) == -y(x)^2/(x*y(x) - 1)
that is, $$y'=\frac{y^2}{1-xy},$$ in agreement with your hand calculation.
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.