[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.
I'm assuming that the covariance matrix is just a diagonal matrix, i.e., sigma(0) = Sigma(0,0), sigma(1) = Sigma(1,1), etc...