Sympy's solve method doesn't return all solutions on this relatively simple problem: import sympy as sp y1, y2, x = sp.symbols('y1 y2 x') y1 = x**2 y2 = x sp.solve(y1, y2, dict = True) This returns [{x: 0}] but there should be an x=1 solution as well. Am I using this correctly? Continue reading...