已解决！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

buffon=function(n, l=0.8, a=1){
  k=0
  theta=runif(n,0,pi); x=runif(n,0,1/2)
  for(t in 1:n){if(x[t]<= l/2*sin(theta[t])) k=k+1 }
  2*l*n/(k*a)
}
buffon(100000,l=0.8,a=1)  # 结果在3.14左右

buf=function(n, l=0.8, a=1){
  k=0
  theta=runif(n,0,pi); x=runif(n,0,1/2)
  for(i in 1:n){ if(x<= 1/2*sin(theta)) k=k+1 }
  2*l*n/(k*a)
}
buf(100000,l=0.8,a=1)   # 结果却在2.51左右

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：已解决 function Theta Run For 程序

我和实验室一师兄花了接近一个小时，没找问题所在，各为也来围观一下

我找到了。

第四行，一个是l，一个是1

qoiqpwqr 发表于 2012-5-26 21:04
第四行，一个是l，一个是1

还是楼主厉害！眼睛太亮了，感谢！找了半天没看到，R里面有一点不好，就是写的程序，不能像C语言、matlab语言、SAS语言那样可以显示代码的语法结构颜色从而提高找错的效率。

楚韵荆风 发表于 2012-5-26 21:18
还是楼主厉害！眼睛太亮了，感谢！找了半天没看到，R里面有一点不好，就是写的程序，不能像C语言、matl ...

如果你用Rstudio的话，会有语法高亮的，数字和字母的颜色是不一样的。

qoiqpwqr 发表于 2012-5-27 04:53
如果你用Rstudio的话，会有语法高亮的，数字和字母的颜色是不一样的。

我发现了，以后准备用这个的，谢谢版主哈