var('y') plot3d(x^2 + y^2 - (1/2)*x^2*y^2, (x, -2, 2), (y, -2, 2), mesh = True)
contour_plot(x^2 + y^2 - (1/2)*x^2*y^2, (x, -2, 2), (y, -2, 2), colorbar = True)
plot3d(x^2 + y^2 - (1/2)*x^2*y, (x, -4, 4), (y, -2, 4), mesh = False, aspect_ratio = [1, 1, 0.2], adaptive = True)
contour_plot(x^2 + y^2 - (1/2)*x^2*y, (x, -3, 3), (y, -3, 3), fill = False, contours = [1, 2, 3, 4, 5, 6, 7, 8], labels = True)