Updated to use c++ sampling functions to fix problem with nose-dividng probabilities
Switched over to using a discrete distribution sampler in c++ from rmultinom
added more graph output including one to plot the log uniform draws
updated to use c++ log uniform draw