Extend medium drawing to arbitrary unwrapping numbers
Allow different periodic system sizes along each axis
Add module for distance finding in non-periodic systems
Bugfix: numpy not imported even though it is used
Replace np.zeros with np.empty for efficiency
Remove unused centre of mass finding algorithm
Change and extend docstrings to match numpy's style