“半个数学系 + 一个计算机科学系 = Deep Learning初级班”spa
simulation = function(sampleSize){ c = rep(0,sampleSize) // <-- 分配了空间 countIn = 0 for(i in 1:sampleSize){ x = runif(1,-1,1) y = runif(1,-1,1) if(sqrt(x*x + y*y) <= 1){ countIn = countIn + 1 } piHat = (countIn / i) * 4 c[i] = piHat } return(c) } sampleSize = 1000 res = simulation(sampleSize)
/* res: 存储着全部的点 */
plot(res[1:sampleSize], type = 'l') // 画点连成的线 lines(rep(pi, sampleSize)[1:sampleSize], col = 'red') // 画PI,参考线