I am new to sympy but I already get a nice output when I plot the implicit function (actually the formula for Cassini's ovals) using sympy: from sympy import plot_implicit, symbols, Eq, solve x, y = symbols('x y') k=2.7 a=3 eq = Eq((x**2 + y**2)**2-2*a**2*(x**2-y**2), k**4-a**4) plot_implicit(eq) Now is it actually possible to somehow get the x and y values corresponding to the plot? or alternatively solve the implicit equation without plotting at all? thanks! Continue reading...