Implemented a change in the functions so that they return more information than before. simsum() still needs to be updated again to fully make more use of this new information.

Also, I fixed a pretty silly bug in my code where the mean of a posterior distribution (normal) was being computed as mu/(mu + sigma^2)