# snippet that can be run after doc_feffit3.lar

_ave = sigma2_eins(150, out.paramgroup.theta, path1_150)
_dlo = sigma2_eins(150, out.paramgroup.theta-out.paramgroup.theta.stderr, path1_150) - _ave
_dhi = sigma2_eins(150, out.paramgroup.theta+out.paramgroup.theta.stderr, path1_150) - _ave
print( "sigma2(T=150) = %.5f  (%+.5f, %+.5f)" % (_ave, _dlo, _dhi))
