No longer uses the initial fit parameters as starting points for the Monte Carlo error estimation, since it sometimes gets stuck in a local minimum or something and gives unreasonably small uncertainties. Also made the main part of the script something that could be called from ipython, for easier debugging.