added n data points variable
added printing chi2
added calculation of actual number of points used in the fit
1sigma and 2sigma contours
added one more pull plot