md_trajectory = nlread('alumina_trajectory.hdf5')[-1]

# Calculate the voidsize distribution of the trajectory.
vd = VoidsizeDistribution(md_trajectory)

# Get the time and autocorrelation values.
voidsizes = vd.voidsizes().inUnitsOf(Ang)
histogram = vd.data()

# Plot the data using pylab.
import pylab

pylab.plot(
    voidsizes,
    histogram,
    label='Voidsize distribution'
)
pylab.xlabel('Voidsizes (Ang)')
pylab.ylabel('Distribution')
pylab.legend()

pylab.show()