var('y')
implicit_plot(x^3 + y^3 == 16, (x, -4, 4), (y, -4, 4))
A = plot3d(exp(x*y), (x, -2.1, 2.1), (y, -2.1, 2.1), aspect_ratio = [20, 20, 1])
B = parametric_plot3d([t, (16 - t^3)^(1/3), exp(t*(16 - t^3)^(1/3))], (t, -2, 2.5), color = 'black')
show(A + B)
C = contour_plot(exp(x*y), (x, -4, 4), (y, -4, 4), contours = [20, 40, 60, 80, 100, 120], fill = False)
D = implicit_plot(x^3 + y^3 == 16, (x, -4, 4), (y, -4, 4), color = 'green')
show(C + D)