Plotting equation (y-x**2)**2 = x**5 with:
from sympy import init_printing, plot_implicit, symbols, Eq
init_printing()
x, y = symbols('x y', real=True)
eq = Eq((y-x**2)**2, x**5)
plot_implicit(eq, (x,0,1), (y,0,.2))
Gives:
The curve should be of constant width. What is the reason for the varying width? how can it be improved?
I tried using adaptive=False, the width is more constant at the cost of a terrible aliasing, and no longer passing through the origin:

