pi的模拟

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

一个比较啰嗦的π的模拟。运用投飞镖法计算圆周率。貌似程序可以写的更简洁些。 require(plotrix)
require(grid)

simp <- function(n){
  x <- runif(n)
  y <- runif(n)
  xy <- x^2 + y^2
  ind <- xy <=1
  ind2 <- ind+1
  up <- length(xy[which(xy <= 1)])
  pai <- up /n * 4
  plot(c(0, 1), c(0, 1), type = 'n', asp = 1,
    main = c('pi = ', paste(pai)))
  rect(0, 1, 1, 0)
  draw.circle(0, 0, 1)
  points(x, y, pch = 20, cex = .5,
    col = c('red', 'blue')[ind2])
}
simp(10)
simp(100)
simp(500)
simp(1000)
simp(2000)
simp(10000)

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：require Equi grid plot REQ function require points 圆周率 程序

图画得很好啊

您好，很想问一下这个图是怎么画的，期待回复：）

mengtong22 发表于 2015-12-29 21:07
您好，很想问一下这个图是怎么画的，期待回复：）

不用回复啦==是我太愚笨，没仔细看：（