v1.0.1 - accepted paper version

Mostly changes in the plotting routine from v1.0.0.
Improved convergence of numerical solution method when combining uncertainty, with better error message upon failure.