from sympy import *
from spb import *
var("x:z")
plot3d_implicit(
    x*y*z < 1, (x, -5, 5), (y, -5, 5), (z, -5, 5),
    n=150, backend=KB
)