region plot x^2 + 3 y^2 = 21