Sympyで♥

from sympy import *

x = Symbol("x")
y = Symbol("y")

exp = x*x + ((y - ((x*x) ** (1/3))) ** 2) - 1
s = solve(exp, y)
p = plot(s[0], s[1], (x,-1.7, 1.7), show=False)
p[0].line_color = p[1].line_color = "#FF14A0"
p.show()