compare ccc op by ton

function (m,n,p0,alpha)
{
arl.ccc = 0
arl.op = 0
rl.ccc = rep(0,m)
rl.op = rep(0,m)
######### compute ############
	ucl.ccc=(log(alpha/2)/log(1-p0))
	lcl.ccc=(log(1-(alpha/2))/log(1-p0))
	gam.al=(log(log(1-alpha/2)/log(alpha/2))/log((alpha/2)/(1-alpha/2)))
	ucl.op=(gam.al*log(alpha/2)/log(1-p0))
	lcl.op=(gam.al*log(1-(alpha/2))/log(1-p0))

cat('\n',ucl.ccc,lcl.ccc,ucl.op,lcl.op,'\n')

########### compare ucl.ccc #############
	for(i in 1:m)
	{
		x=rep(0,n)
      	x=rgeom(n,p0)
		in.control=TRUE
		j=0
		k=0

	while(in.control)
		{
			j=j+1
			k=k+1
		if((x[k]>ucl.ccc)||(x[k]<lcl.ccc))
			{
			in.control=FALSE
			}
		}
		rl.ccc[i]=j
	}
########### compare ucl.op #############
for(i in 1:m)
	{
		x=rep(0,n)
      	x=rgeom(n,p0)
		in.control=TRUE
		j=0
		k=0

	while(in.control)
		{
			j=j+1
			k=k+1
		if((x[k]>ucl.op)||(x[k]<lcl.op))
			{
			in.control=FALSE
			}
		}
		rl.op[i]=j
	}

######### print ###########################
arl.ccc=mean(rl.ccc)
arl.op=mean(rl.op)
cat('\n','arl.ccc =',arl.ccc,'arl.op =',arl.op,'\n')
}