求pi值的编程问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

以下是编写的程序

clear
local total = 0
forvalues k = 1/1000 {
    local v1 = 2 * runiform() - 1
    local v2 = 2 * runiform() - 1
    if `v1'^2 + `v2'^2 < = 1 {
        local I = 1
    }
    else {
        local I = 0
    }
    local total = `total' + `I'
}
local pai = 4 * `total'/1000
dis "`pai'"

为什么求出的pi值是3.7？ 哪里出错了？ 多谢！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：forvalues forvalue Uniform values Local

模拟1000次, 是不是太少了? 100000次呢? Stata, 没研究过, 不好意思

*可简化作：
loc N=1000
loc n=0
forv i=1/`N'{
loc n=`n'+(runiform()^2+runiform()^2<=1)
}
di "pi="4*`n'/`N'

*不指定种子，每次运行的结果一般不同。

另一种试验（Buffon's needle problem）
loc L=10
loc l=8
loc N=1000
loc n=0
forv i=1/`N'{
loc n=`n'+(`L'*runiform()<=`l'*sin(_pi*runiform()))
}
di "pi="2*`l'*`N'/(`L'*`n')

多谢楼上和楼上的楼上的回复。
改为：(`v1')^2 + (`v2')^2 < = 1 后，pi值变成3.14左右了。