Added uncertainty bands to fit; redo fit up to MaxAttempts times to ensure congergence; print boundaries; print edge position on -log(L) plot; added load and unload functions for samples to reduce memory footprint (still experimental, currently commented out);