baci<-function(p,N)
{
brp=0
for(i in 1:N)
{
x=runif(1)
if (x<p) brp=brp+1
}
return(brp/N-p)
}

freq=c()
n=100
p=1/2
for(i in 1:n)
{
freq=c(freq,baci(p,100))
}


igra<-function(n=40,a=20,p=1/2)
{
trajektorija=c(a)
trenutno=a
for(i in 1:n)
{
x=runif(1)
if(x<p) trenutno=trenutno+1
else trenutno=trenutno-1
trajektorija=c(trajektorija,trenutno)
}
return(trajektorija)
}


oceniNula<-function(n=40,a=0,b=0,p=1/2,N)
{
brojac=0
for(i in 1:N)
{
trajektorija=igra(n,a,p)
if (trajektorija[(n+1)]==b) brojac= brojac+1
}
return(brojac/N)
}

brDece<-function(Nmax=10,p=1/2)
{
OK=TRUE
brojac=0
while(OK&&(brojac<(Nmax+1)))
{
x=runif(1)
if(x<p) OK=FALSE
brojac=brojac+1
}
return(brojac)
}

ocenaBRdece<-function(N,p=1/2,Nmax=10)
{
br=0
for(i in 1:N)
{
br=br+brDece(Nmax,p)
}
return(br/N)
}

niz=c()
p=seq(from=0.45,to=0.55,length=100)
for(i in 1:100)
niz=c(niz,ocenaBRdece(10000,p[i],10))



brDece1<-function(Nmax=10,p=1/2)
{
OKm=TRUE
OKf=TRUE
brojac=0
while((OKm||OKf)&&(brojac<(Nmax+1)))
{
x=runif(1)
if(x<p) {OKm=FALSE}
else {OKf=FALSE}
brojac=brojac+1
}
return(brojac)
}

ocenaBRdece1<-function(N,p=1/2,Nmax=10)
{
br=0
for(i in 1:N)
{
br=br+brDece1(Nmax,p)
}
return(br/N)
}

niz1=c()
p=seq(from=0.45,to=0.55,length=100)
for(i in 1:100)
niz1=c(niz1,ocenaBRdece1(10000,p[i],10))