Monte Carlo Simulations

[alcarney]

Once #12 has been completed it would be good as a comparison to the polynomial chaos method to look at doing a Monte-Carlo simulation on the previous problems.

This means:

• Sampling various values of the parameters a, b and c from a distribution and solving the resulting system.
• Computing things like the mean and standard deviation of the solution

This will give a good comparison to the polynomial chaos approach and should be relatively easy to do with the existing codebase.

[alcarney]

Would've been nice but at this stage no time :(