`
yangyi
  • 浏览: 113211 次
  • 性别: Icon_minigender_1
  • 来自: 北京
社区版块
存档分类
最新评论

用蒙特卡罗法求圆周率

阅读更多
用蒙特卡罗法求解圆周率,很简单,很慢。问题是怎么才能知道已经达到了所需要的精度,还是需要利用概率论对数据进行加工

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <time.h>

int main(int argc, char** argv){
	double insideCount=0, loopCount=300000000, i;
	srand((int)time(0));

	for(i=0;i<loopCount;i++){
		if(sqrt(pow(rand(), 2) + pow(rand(), 2)) <= RAND_MAX){
			insideCount++;
		}
	}

	printf("%f", 4*(insideCount/loopCount);
}
分享到:
评论

相关推荐

Global site tag (gtag.js) - Google Analytics