"""

=========================

voxels([x, y, z], filled)

=========================

See `~mpl_toolkits.mplot3d.axes3d.Axes3D.voxels`.

"""

import matplotlib.pyplot as plt

import numpy as np

plt.style.use('_mpl-gallery')

# Prepare some coordinates

x, y, z = np.indices((8, 8, 8))

# Draw cuboids in the top left and bottom right corners

cube1 = (x < 3) & (y < 3) & (z < 3)

cube2 = (x >= 5) & (y >= 5) & (z >= 5)

# Combine the objects into a single boolean array

voxelarray = cube1 | cube2

# Plot

fig, ax = plt.subplots(subplot_kw={"projection": "3d"})

ax.voxels(voxelarray, edgecolor='k')

ax.set(xticklabels=[],

yticklabels=[],

zticklabels=[])

plt.show()