plot3d(lambda x, y: x * np.exp(-x**2 - y**2),
    ("x", -3, 3), ("y", -3, 3), use_cm=True)  # doctest: +SKIP