*******************************************************************************
Thu Mar 28 13:18:01 2019


FIT:    data read from 'compile_spin_data.dat' u  ($4*1000):(($9-zero)/zero*1.e6  > -5. && abs($8)<0.0001 ? (($9-zero)/zero*1.e6):1/0)
        format = x:z
        #datapoints = 17
        residuals are weighted equally (unit weight)

function used for fitting: g(x)
	g(x) = -a*x**2 +b*x +c + d*x**3
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a             b             c             d            
   0 3.9755274419e+11   0.00e+00  1.01e+03    5.800000e-03   1.000000e-01   1.000000e-01   1.000000e-01
   7 6.1712841924e+10  -1.86e-02  1.01e-04    2.208610e+03  -2.296081e+04  -4.376152e+04  -4.975470e+01

After 7 iterations the fit converged.
final sum of squares of residuals : 6.17128e+10
rel. change during last iteration : -1.86022e-07

degrees of freedom    (FIT_NDF)                        : 13
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 68899.5
variance of residuals (reduced chisquare) = WSSR/ndf   : 4.74714e+09

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = 2208.61          +/- 709.3        (32.11%)
b               = -22960.8         +/- 1.133e+04    (49.33%)
c               = -43761.5         +/- 4.276e+04    (97.72%)
d               = -49.7547         +/- 12.05        (24.23%)

correlation matrix of the fit parameters:
                a      b      c      d      
a               1.000 
b              -0.953  1.000 
c              -0.704  0.838  1.000 
d              -0.981  0.881  0.608  1.000 


*******************************************************************************
Thu Mar 28 13:18:03 2019


FIT:    data read from 'compile_spin_data.dat' u ($4*1000):($26 != 0 ? (-$4**2 / (betay)**2/4. * 1.e6):1/0)
        format = x:z
        #datapoints = 35
        residuals are weighted equally (unit weight)

function used for fitting: h(x)
	h(x) = -a1*x**2
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a1           
   0 1.7898639529e-05   0.00e+00  2.30e-01    5.800000e-04
   5 1.5455587804e-32   0.00e+00  2.30e-06    5.818044e-04

After 5 iterations the fit converged.
final sum of squares of residuals : 1.54556e-32
abs. change during last iteration : 0


Hmmmm.... Sum of squared residuals is zero. Can't compute errors.

Final set of parameters 
======================= 

a1              = 0.000581804    


*******************************************************************************
Mon Apr  8 15:48:26 2019


FIT:    data read from 'compile_spin_data.dat' u  ($4*1000):(($9-zero)/zero*1.e6  > -5. && abs($8)<0.0001 ? (($9-zero)/zero*1.e6):1/0)
        format = x:z
        #datapoints = 69
        residuals are weighted equally (unit weight)

function used for fitting: g(x)
	g(x) = -a*x**2 +b*x +c + d*x**3
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a             b             c             d            
   0 7.7444415861e+11   0.00e+00  1.01e+03    5.800000e-03   1.000000e-01   1.000000e-01   1.000000e-01
   6 5.8375182715e+11  -3.36e-01  1.01e-03   -1.046642e+02   9.806040e+00  -1.448924e+04  -1.667300e-02

After 6 iterations the fit converged.
final sum of squares of residuals : 5.83752e+11
rel. change during last iteration : -3.36025e-06

degrees of freedom    (FIT_NDF)                        : 65
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 94767.1
variance of residuals (reduced chisquare) = WSSR/ndf   : 8.9808e+09

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = -104.664         +/- 24.21        (23.13%)
b               = 9.80604          +/- 1474         (1.503e+04%)
c               = -14489.2         +/- 1.366e+04    (94.31%)
d               = -0.016673        +/- 1.291        (7744%)

correlation matrix of the fit parameters:
                a      b      c      d      
a               1.000 
b              -0.000  1.000 
c               0.550  0.000  1.000 
d               0.000 -0.898 -0.000  1.000 


*******************************************************************************
Mon Apr  8 15:48:27 2019


FIT:    data read from 'compile_spin_data.dat' u ($4*1000):($26 != 0 ? (-$4**2 / (betay)**2/4. * 1.e6):1/0)
        format = x:z
        #datapoints = 117
        residuals are weighted equally (unit weight)

function used for fitting: h(x)
	h(x) = -a1*x**2
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a1           
   0 1.5041560833e-03   0.00e+00  2.51e-01    5.800000e-04
   5 5.2987946931e-32   0.00e+00  2.51e-06    5.717292e-04

After 5 iterations the fit converged.
final sum of squares of residuals : 5.29879e-32
abs. change during last iteration : 0


Hmmmm.... Sum of squared residuals is zero. Can't compute errors.

Final set of parameters 
======================= 

a1              = 0.000571729    


*******************************************************************************
Mon Apr  8 15:48:31 2019


FIT:    data read from 'compile_spin_data.dat' u ($8*100):($11/$13*$8*f)
        format = x:z
        x range restricted to [-0.350000 : 0.350000]
        #datapoints = 121
        residuals are weighted equally (unit weight)

function used for fitting: h(x)
	h(x) = a*x**2+b*x+c + a2*cos(b2*x) +d*x**3  +e*x**4
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a             b             c             d             e            
   0 1.3647302411e+02   0.00e+00  1.43e-02    1.000000e+00   1.000000e-02   1.000000e-02   1.000000e-02   1.000000e-02
   6 2.0825727240e-03  -2.18e-01  1.43e-08    4.481877e+00  -9.015166e-03  -9.988530e-01   3.279350e+00  -1.153779e+02

After 6 iterations the fit converged.
final sum of squares of residuals : 0.00208257
rel. change during last iteration : -2.18481e-06

degrees of freedom    (FIT_NDF)                        : 116
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.00423712
variance of residuals (reduced chisquare) = WSSR/ndf   : 1.79532e-05

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = 4.48188          +/- 0.05705      (1.273%)
b               = -0.00901517      +/- 0.0116       (128.7%)
c               = -0.998853        +/- 0.0004749    (0.04755%)
d               = 3.27935          +/- 0.2761       (8.42%)
e               = -115.378         +/- 1.076        (0.9323%)

correlation matrix of the fit parameters:
                a      b      c      d      e      
a               1.000 
b              -0.547  1.000 
c              -0.393  0.069  1.000 
d               0.625 -0.942 -0.074  1.000 
e              -0.921  0.711  0.217 -0.829  1.000 


*******************************************************************************
Mon Apr  8 15:48:31 2019


FIT:    data read from 'compile_spin_data.dat' u ($8*100):(f($8))*1.e6
        format = x:z
        x range restricted to [-0.350000 : 0.350000]
        #datapoints = 121
        residuals are weighted equally (unit weight)

function used for fitting: n(x)
	n(x) =a3*cos(b3*x)-a3
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a3            b3           
   0 5.6044233832e-03   0.00e+00  1.23e+00    2.279070e+02   4.753540e-01
3923 3.7105132196e-06  -5.78e-01  1.23e-03    4.507938e+02   3.364000e-01

After 3923 iterations the fit converged.
final sum of squares of residuals : 3.71051e-06
rel. change during last iteration : -5.77589e-06

degrees of freedom    (FIT_NDF)                        : 119
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.000176581
variance of residuals (reduced chisquare) = WSSR/ndf   : 3.11808e-08

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a3              = 450.794          +/- 41.46        (9.198%)
b3              = 0.3364           +/- 0.01547      (4.6%)

correlation matrix of the fit parameters:
                a3     b3     
a3              1.000 
b3             -1.000  1.000 


*******************************************************************************
Mon Apr  8 15:49:02 2019


FIT:    data read from 'compile_spin_data.dat' u ($8*100):($11/$13*$8*f)
        format = x:z
        x range restricted to [-0.350000 : 0.350000]
        #datapoints = 121
        residuals are weighted equally (unit weight)

function used for fitting: h(x)
	h(x) = a*x**2+b*x+c + a2*cos(b2*x) +d*x**3  +e*x**4
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a             b             c             d             e            
   0 1.3647302411e+02   0.00e+00  1.43e-02    1.000000e+00   1.000000e-02   1.000000e-02   1.000000e-02   1.000000e-02
   6 2.0825727240e-03  -2.18e-01  1.43e-08    4.481877e+00  -9.015166e-03  -9.988530e-01   3.279350e+00  -1.153779e+02

After 6 iterations the fit converged.
final sum of squares of residuals : 0.00208257
rel. change during last iteration : -2.18481e-06

degrees of freedom    (FIT_NDF)                        : 116
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.00423712
variance of residuals (reduced chisquare) = WSSR/ndf   : 1.79532e-05

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = 4.48188          +/- 0.05705      (1.273%)
b               = -0.00901517      +/- 0.0116       (128.7%)
c               = -0.998853        +/- 0.0004749    (0.04755%)
d               = 3.27935          +/- 0.2761       (8.42%)
e               = -115.378         +/- 1.076        (0.9323%)

correlation matrix of the fit parameters:
                a      b      c      d      e      
a               1.000 
b              -0.547  1.000 
c              -0.393  0.069  1.000 
d               0.625 -0.942 -0.074  1.000 
e              -0.921  0.711  0.217 -0.829  1.000 


*******************************************************************************
Mon Apr  8 15:49:02 2019


FIT:    data read from 'compile_spin_data.dat' u ($8*100):(f($8))*1.e6
        format = x:z
        x range restricted to [-0.350000 : 0.350000]
        #datapoints = 121
        residuals are weighted equally (unit weight)

function used for fitting: n(x)
	n(x) =a3*cos(b3*x)-a3
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a3            b3           
   0 5.6044233832e-03   0.00e+00  1.23e+00    2.279070e+02   4.753540e-01
3923 3.7105132196e-06  -5.78e-01  1.23e-03    4.507938e+02   3.364000e-01

After 3923 iterations the fit converged.
final sum of squares of residuals : 3.71051e-06
rel. change during last iteration : -5.77589e-06

degrees of freedom    (FIT_NDF)                        : 119
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.000176581
variance of residuals (reduced chisquare) = WSSR/ndf   : 3.11808e-08

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a3              = 450.794          +/- 41.46        (9.198%)
b3              = 0.3364           +/- 0.01547      (4.6%)

correlation matrix of the fit parameters:
                a3     b3     
a3              1.000 
b3             -1.000  1.000 


*******************************************************************************
Mon Apr  8 18:46:05 2019


FIT:    data read from 'compile_spin_data.dat' u  ($4*1000):(($9-zero)/zero*1.e6  > -5. && abs($8)<0.0001 ? (($9-zero)/zero*1.e6):1/0)
        format = x:z
        #datapoints = 69
        residuals are weighted equally (unit weight)

function used for fitting: g(x)
	g(x) = -a*x**2 +b*x +c + d*x**3
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a             b             c             d            
   0 7.7444415861e+11   0.00e+00  1.01e+03    5.800000e-03   1.000000e-01   1.000000e-01   1.000000e-01
   6 5.8375182715e+11  -3.36e-01  1.01e-03   -1.046642e+02   9.806040e+00  -1.448924e+04  -1.667300e-02

After 6 iterations the fit converged.
final sum of squares of residuals : 5.83752e+11
rel. change during last iteration : -3.36025e-06

degrees of freedom    (FIT_NDF)                        : 65
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 94767.1
variance of residuals (reduced chisquare) = WSSR/ndf   : 8.9808e+09

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = -104.664         +/- 24.21        (23.13%)
b               = 9.80604          +/- 1474         (1.503e+04%)
c               = -14489.2         +/- 1.366e+04    (94.31%)
d               = -0.016673        +/- 1.291        (7744%)

correlation matrix of the fit parameters:
                a      b      c      d      
a               1.000 
b              -0.000  1.000 
c               0.550  0.000  1.000 
d               0.000 -0.898 -0.000  1.000 


*******************************************************************************
Mon Apr  8 18:46:06 2019


FIT:    data read from 'compile_spin_data.dat' u ($4*1000):($26 != 0 ? (-$4**2 / (betay)**2/4. * 1.e6):1/0)
        format = x:z
        #datapoints = 117
        residuals are weighted equally (unit weight)

function used for fitting: h(x)
	h(x) = -a1*x**2
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a1           
   0 1.5041560833e-03   0.00e+00  2.51e-01    5.800000e-04
   5 5.2987946931e-32   0.00e+00  2.51e-06    5.717292e-04

After 5 iterations the fit converged.
final sum of squares of residuals : 5.29879e-32
abs. change during last iteration : 0


Hmmmm.... Sum of squared residuals is zero. Can't compute errors.

Final set of parameters 
======================= 

a1              = 0.000571729    


*******************************************************************************
Mon Apr  8 18:46:37 2019


FIT:    data read from 'compile_spin_data.dat' u  ($4*1000):(($9-zero)/zero*1.e6  > -5. && abs($8)<0.0001 ? (($9-zero)/zero*1.e6):1/0)
        format = x:z
        #datapoints = 69
        residuals are weighted equally (unit weight)

function used for fitting: g(x)
	g(x) = -a*x**2 +b*x +c + d*x**3
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a             b             c             d            
   0 7.7444415861e+11   0.00e+00  1.01e+03    5.800000e-03   1.000000e-01   1.000000e-01   1.000000e-01
   6 5.8375182715e+11  -3.36e-01  1.01e-03   -1.046642e+02   9.806040e+00  -1.448924e+04  -1.667300e-02

After 6 iterations the fit converged.
final sum of squares of residuals : 5.83752e+11
rel. change during last iteration : -3.36025e-06

degrees of freedom    (FIT_NDF)                        : 65
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 94767.1
variance of residuals (reduced chisquare) = WSSR/ndf   : 8.9808e+09

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = -104.664         +/- 24.21        (23.13%)
b               = 9.80604          +/- 1474         (1.503e+04%)
c               = -14489.2         +/- 1.366e+04    (94.31%)
d               = -0.016673        +/- 1.291        (7744%)

correlation matrix of the fit parameters:
                a      b      c      d      
a               1.000 
b              -0.000  1.000 
c               0.550  0.000  1.000 
d               0.000 -0.898 -0.000  1.000 


*******************************************************************************
Mon Apr  8 18:46:40 2019


FIT:    data read from 'compile_spin_data.dat' u ($4*1000):($26 != 0 ? (-$4**2 / (betay)**2/4. * 1.e6):1/0)
        format = x:z
        #datapoints = 117
        residuals are weighted equally (unit weight)

function used for fitting: h(x)
	h(x) = -a1*x**2
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a1           
   0 1.5041560833e-03   0.00e+00  2.51e-01    5.800000e-04
   5 5.2987946931e-32   0.00e+00  2.51e-06    5.717292e-04

After 5 iterations the fit converged.
final sum of squares of residuals : 5.29879e-32
abs. change during last iteration : 0


Hmmmm.... Sum of squared residuals is zero. Can't compute errors.

Final set of parameters 
======================= 

a1              = 0.000571729    


*******************************************************************************
Mon Apr  8 18:51:48 2019


FIT:    data read from 'compile_spin_data.dat' u  ($4*1000):(($9-zero)/zero*1.e6  > -5. && abs($8)<0.0001 ? (($9-zero)/zero*1.e6):1/0)
        format = x:z
        #datapoints = 69
        residuals are weighted equally (unit weight)

function used for fitting: g(x)
	g(x) = -a*x**2 +b*x +c + d*x**3
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a             b             c             d            
   0 7.7444415861e+11   0.00e+00  1.01e+03    5.800000e-03   1.000000e-01   1.000000e-01   1.000000e-01
   6 5.8375182715e+11  -3.36e-01  1.01e-03   -1.046642e+02   9.806040e+00  -1.448924e+04  -1.667300e-02

After 6 iterations the fit converged.
final sum of squares of residuals : 5.83752e+11
rel. change during last iteration : -3.36025e-06

degrees of freedom    (FIT_NDF)                        : 65
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 94767.1
variance of residuals (reduced chisquare) = WSSR/ndf   : 8.9808e+09

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = -104.664         +/- 24.21        (23.13%)
b               = 9.80604          +/- 1474         (1.503e+04%)
c               = -14489.2         +/- 1.366e+04    (94.31%)
d               = -0.016673        +/- 1.291        (7744%)

correlation matrix of the fit parameters:
                a      b      c      d      
a               1.000 
b              -0.000  1.000 
c               0.550  0.000  1.000 
d               0.000 -0.898 -0.000  1.000 


*******************************************************************************
Mon Apr  8 18:51:50 2019


FIT:    data read from 'compile_spin_data.dat' u ($4*1000):($26 != 0 ? (-$4**2 / (betay)**2/4. * 1.e6):1/0)
        format = x:z
        #datapoints = 117
        residuals are weighted equally (unit weight)

function used for fitting: h(x)
	h(x) = -a1*x**2
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a1           
   0 1.5041560833e-03   0.00e+00  2.51e-01    5.800000e-04
   5 5.2987946931e-32   0.00e+00  2.51e-06    5.717292e-04

After 5 iterations the fit converged.
final sum of squares of residuals : 5.29879e-32
abs. change during last iteration : 0


Hmmmm.... Sum of squared residuals is zero. Can't compute errors.

Final set of parameters 
======================= 

a1              = 0.000571729    


*******************************************************************************
Mon Apr  8 18:52:50 2019


FIT:    data read from 'compile_spin_data.dat' u  ($4*1000):(($9-zero)/zero*1.e6  > -5. && abs($8)<0.0001 ? (($9-zero)/zero*1.e6):1/0)
        format = x:z
        #datapoints = 69
        residuals are weighted equally (unit weight)

function used for fitting: g(x)
	g(x) = -a*x**2 +b*x +c + d*x**3
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a             b             c             d            
   0 7.7444415861e+11   0.00e+00  1.01e+03    5.800000e-03   1.000000e-01   1.000000e-01   1.000000e-01
   6 5.8375182715e+11  -3.36e-01  1.01e-03   -1.046642e+02   9.806040e+00  -1.448924e+04  -1.667300e-02

After 6 iterations the fit converged.
final sum of squares of residuals : 5.83752e+11
rel. change during last iteration : -3.36025e-06

degrees of freedom    (FIT_NDF)                        : 65
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 94767.1
variance of residuals (reduced chisquare) = WSSR/ndf   : 8.9808e+09

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = -104.664         +/- 24.21        (23.13%)
b               = 9.80604          +/- 1474         (1.503e+04%)
c               = -14489.2         +/- 1.366e+04    (94.31%)
d               = -0.016673        +/- 1.291        (7744%)

correlation matrix of the fit parameters:
                a      b      c      d      
a               1.000 
b              -0.000  1.000 
c               0.550  0.000  1.000 
d               0.000 -0.898 -0.000  1.000 


*******************************************************************************
Mon Apr  8 18:52:51 2019


FIT:    data read from 'compile_spin_data.dat' u ($4*1000):($26 != 0 ? (-$4**2 / (betay)**2/4. * 1.e6):1/0)
        format = x:z
        #datapoints = 117
        residuals are weighted equally (unit weight)

function used for fitting: h(x)
	h(x) = -a1*x**2
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a1           
   0 1.5041560833e-03   0.00e+00  2.51e-01    5.800000e-04
   5 5.2987946931e-32   0.00e+00  2.51e-06    5.717292e-04

After 5 iterations the fit converged.
final sum of squares of residuals : 5.29879e-32
abs. change during last iteration : 0


Hmmmm.... Sum of squared residuals is zero. Can't compute errors.

Final set of parameters 
======================= 

a1              = 0.000571729    


*******************************************************************************
Mon Apr  8 18:53:38 2019


FIT:    data read from 'compile_spin_data.dat' u  ($4*1000):(($9-zero)/zero*1.e6  > -50. &&($9-zero)/zero*1.e6 < 50. && abs($8)<0.0001 ? (($9-zero)/zero*1.e6):1/0)
        format = x:z
        #datapoints = 67
        residuals are weighted equally (unit weight)

function used for fitting: g(x)
	g(x) = -a*x**2 +b*x +c + d*x**3
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a             b             c             d            
   0 1.9780899808e+08   0.00e+00  8.58e+02    5.800000e-03   1.000000e-01   1.000000e-01   1.000000e-01
   7 8.5640716555e-02  -1.25e-06  8.58e-05    6.566589e-04   4.478418e-10  -2.294944e-03  -3.826486e-13

After 7 iterations the fit converged.
final sum of squares of residuals : 0.0856407
rel. change during last iteration : -1.24991e-11

degrees of freedom    (FIT_NDF)                        : 63
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.0368697
variance of residuals (reduced chisquare) = WSSR/ndf   : 0.00135938

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = 0.000656659      +/- 1.069e-05    (1.627%)
b               = 4.47842e-10      +/- 0.0005958    (1.33e+08%)
c               = -0.00229494      +/- 0.005362     (233.7%)
d               = -3.82649e-13     +/- 5.726e-07    (1.496e+08%)

correlation matrix of the fit parameters:
                a      b      c      d      
a               1.000 
b              -0.000  1.000 
c               0.543 -0.000  1.000 
d               0.000 -0.889  0.000  1.000 


*******************************************************************************
Mon Apr  8 18:53:39 2019


FIT:    data read from 'compile_spin_data.dat' u ($4*1000):($26 != 0 ? (-$4**2 / (betay)**2/4. * 1.e6):1/0)
        format = x:z
        #datapoints = 117
        residuals are weighted equally (unit weight)

function used for fitting: h(x)
	h(x) = -a1*x**2
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a1           
   0 1.5041560833e-03   0.00e+00  2.51e-01    5.800000e-04
   5 5.2987946931e-32   0.00e+00  2.51e-06    5.717292e-04

After 5 iterations the fit converged.
final sum of squares of residuals : 5.29879e-32
abs. change during last iteration : 0


Hmmmm.... Sum of squared residuals is zero. Can't compute errors.

Final set of parameters 
======================= 

a1              = 0.000571729    


*******************************************************************************
Mon Apr  8 18:54:17 2019


FIT:    data read from 'compile_spin_data.dat' u  ($4*1000):(($9-zero)/zero*1.e6  > -50. &&($9-zero)/zero*1.e6 < 50. && abs($8)<0.0001 ? (($9-zero)/zero*1.e6):1/0)
        format = x:z
        #datapoints = 67
        residuals are weighted equally (unit weight)

function used for fitting: g(x)
	g(x) = -a*x**2 +b*x +c + d*x**3
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a             b             c             d            
   0 1.9780899808e+08   0.00e+00  8.58e+02    5.800000e-03   1.000000e-01   1.000000e-01   1.000000e-01
   7 8.5640716555e-02  -1.25e-06  8.58e-05    6.566589e-04   4.478418e-10  -2.294944e-03  -3.826486e-13

After 7 iterations the fit converged.
final sum of squares of residuals : 0.0856407
rel. change during last iteration : -1.24991e-11

degrees of freedom    (FIT_NDF)                        : 63
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.0368697
variance of residuals (reduced chisquare) = WSSR/ndf   : 0.00135938

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = 0.000656659      +/- 1.069e-05    (1.627%)
b               = 4.47842e-10      +/- 0.0005958    (1.33e+08%)
c               = -0.00229494      +/- 0.005362     (233.7%)
d               = -3.82649e-13     +/- 5.726e-07    (1.496e+08%)

correlation matrix of the fit parameters:
                a      b      c      d      
a               1.000 
b              -0.000  1.000 
c               0.543 -0.000  1.000 
d               0.000 -0.889  0.000  1.000 


*******************************************************************************
Mon Apr  8 18:54:18 2019


FIT:    data read from 'compile_spin_data.dat' u ($4*1000):($26 != 0 ? (-$4**2 / (betay)**2/4. * 1.e6):1/0)
        format = x:z
        #datapoints = 117
        residuals are weighted equally (unit weight)

function used for fitting: h(x)
	h(x) = -a1*x**2
fitted parameters initialized with current variable values

iter      chisq       delta/lim  lambda   a1           
   0 1.5041560833e-03   0.00e+00  2.51e-01    5.800000e-04
   5 5.2987946931e-32   0.00e+00  2.51e-06    5.717292e-04

After 5 iterations the fit converged.
final sum of squares of residuals : 5.29879e-32
abs. change during last iteration : 0


Hmmmm.... Sum of squared residuals is zero. Can't compute errors.

Final set of parameters 
======================= 

a1              = 0.000571729