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