Added new module diffusion.py

So far, contains only one routine to compute the diffusion coefficient of a
collection of particles in 3D.