# Function to simulate coin tosses
Coin_Toss<-function(tosses){
  a<-rbinom(tosses,1,0.5)
  b<-1:tosses
  c<-cumsum(a)
  d<-c/b
  return(d[tosses])
}

# Simulate 1000 Experiements with 500 tosses

sample<-replicate(1000,Coin_Toss(tosses=500))
summary(sample)
hist(sample,col="Orange",xlab="Probability of Head",main = "Probability of Head (1000 Experiments of Sample Size = 500)")
sd(sample)