[bdg-utils-13] Added code to implement a poisson mixture model.

Adds code to fit a mixture of poisson random variables. Additionally, adds a few mathematical utils. Closes #13.
This should probably be made explicit then.  Also, why limit yourself to a diagonal matrix?