E.g. find critical points on a surface:
e = (x-y+2)*(x**2+y**2-9) solve([e.diff(x), e.diff(y)], [x, y])