R で書いたモンテカルロ法

#
# この中身を I:￥rdata￥monsim.r に保存してください
#
mon.sim <- function(ite){
  sc <- numeric(ite)
  for(i in 1:ite){
    x<- runif(200)
    y<-ifelse(x<0.4,1,ifelse(x<0.7,2,
       ifelse(x<0.9,3,4)))  # 本当は１行
    z1 <- table(y)
    z2 <- c(40,30,20,10)*2
    z3 <- (z1-z2)^2 / z2
    sc[i] <- sum(z3)
  }
  return(sc)
}