ContourPlot[ y^2 + x^2 == 1 - 1337*x^2*y^2, {x, -1.3, 1.3}, {y, -1.3, 1.3}]