R Code (Yoye)

function (m,n,p,alpha,m2,n2,round1)
{
    set.seed(5)
    x=rbinom(m,n,p)
    z=length(x)

        for(i in 1:z){
            p[i]=(x[i]/m)
        }
        pbar=(sum(p)/n)
        cat('\n','M =',p,' ','\n')
        Mbar=pbar
cat ('\n','M bar =',Mbar,'\n')

####################### Calculate UCL,CL,LCL ##############################
######## Calculate CAL(L) ##########
    CLl = (Mbar*alpha)
    LCLl= CLl-(3*(sqrt(((CLl)*(1-CLl))/n)))
    UCLl= CLl+(3*(sqrt(((CLl)*(1-CLl))/n)))

cat('\n','--------------- Real UCL,CL,LCL for (L) ------------------','\n')
cat('\n','UCL(L) =',UCLl,'CL(L) =',CLl,'LCL(L) =',LCLl,'\n')


######## Calculate CAL(R) ##########
    CLr = 1-((1-((Mbar*alpha)*alpha)))
    LCLr= CLr-(3*(sqrt(((CLr)*(1-CLr))/n)))
    UCLr= CLr+(3*(sqrt(((CLr)*(1-CLr))/n)))

cat('\n','--------------- Real UCL,CL,LCL for (R) ------------------','\n')
cat('\n','UCL(R) =',UCLr,'CL(R) =',CLr,'LCL(R) =',LCLr,'\n')

######## Check Condition for UCL,LCL (L) ##########
        if(LCLl <= 0){
            LCLl=0
        } else LCLl=LCLl
        if(UCLl > 1){
            UCLl=1
        } else UCLl=UCLl

######## Check Condition for UCL,LCL (R) ##########
        if(LCLr <= 0){
            LCLr=0
        }else LCLr=LCLr

        if(UCLr > 1){
            UCLr=1
        } else UCLr=UCLr

cat('\n','--------------- Calculate UCL,CL,LCL for (L) ------------------','\n')
cat('\n','UCL(L) =',UCLl,'CL(L) =',CLl,'LCL(L) =',LCLl,'\n')
cat('\n','--------------- Calculate UCL,CL,LCL for (R) ------------------','\n')
cat('\n','UCL(R) =',UCLr,'CL(R) =',CLr,'LCL(R) =',LCLr,'\n')


######### Check Condition Process (Not work!) #########
######## Calculate Step 2 ##############

    for (i in 1:round1){
        x2=rbinom(m2,n2,Mbar)
        z1=length(x2)
        pnew=rep(0,m2)
        rl=rep(0,round1)

    for (v in 1:z1){
        pnew[v]=(x2[v]/m2)
    }
#cat('\n',pnew,'\n')
#cat('\n',length(pnew),'\n')

        in.control=TRUE
            j=0
            k=0

        while(in.control)
            {
                j=j+1
                k=k+1
            if((pnew[k]>UCLl)||(pnew[k]<LCLl))
                {
                in.control=FALSE
                }
            }
            rl[i]=j

}
#cat('\n',rl,'\n')
}