BF: correction to fit cumulative normal function 

Fitting a cumulative normal function gave incorrect results on recovering parameters from simulated data. I looked up the function on http://en.wikipedia.org/wiki/Normal_distribution and there seemed to be some differences. I updated to code according to this definition and did so for the inverse function and recovering the parameters from simulated data worked fine.