m=100;n=1:100
ln=length(n);
for(i in 1:ln){
  bi=0;ci=0;ni=n[i]
  for(j in 1:ni){
    xj=rnorm(m)
    yj=shapiro.test(xj)$p.value
    if(yj<0.05){
      bi=bi+1
    }
    else{
      ci=ci+1
    }
  }
  print(bi/ni)
}