*! version 1.0.0 STB-34 sts12
program define bartcdf
	version 4.0

	local a `1'
	local eps=.00001

	if `a' < 0.3 {
		local prob=0.
	}
	else {
		local prob=1.
		local i=1

		while `i' <= 100 {

			local del=2*(-1.)^`i'*exp(-2*`a'^2*`i'^2)
*			di `i' " " `prob' " " `del'

			if abs(`del') < `eps'*`prob' { 
				local i=101
			}
			else {
				local prob=`prob'+`del'
				local i=`i'+1
			}
		}

	}

	global S_1=`prob'
*	di `prob'
end