Cube lattice

using AbstractPlotting

r = range(-3, stop = 3, length = 100)
me = [((1 ./ x).^2 + (1 ./ y).^2 + (1 ./ z).^2) for x=r, y=r, z=r]
me2 = me .* (abs.(me) .> 1.5)
contour(me2, colormap = :Set2)