ContourPlot[x^2 + 3 x y + 3 y^2 == 1, {x, -2.5, 2.5}, {y, -2.5, 2.5}]