regionplot(x^2 + y^2 - 1)^3 - x^2 y^3 <= 0 red