蒙特卡罗算法计算圆周率π

(1)均匀撒点:在C语言中可以使用随机函数来实现,产生【0,1】之间的随机坐标【x,y】。
(2)区域判断:图中阴影部分的特点是距离坐标原点的距离小于等于1,这样,可以通过计算判断x2+y2<=1来实现。

#include<stdio.h>
#include<time.h>
#include<stdlib.h>
double MontePI(int n)
{
	double PI;
	double x,y;
	int i,sum;
	sum=0;
	srand(time(NULL));
	for(i=1;i<n;i++)
	{
		x=(double)rand()/RAND_MAX;//产生0~1之间的一个随机数
		y=(double)rand()/RAND_MAX;//产生0~1之间的一个随机数
		if(x*x+y*y<=1)//若在阴影区域
			sum++;//计数
	}
	PI=4.0*sum/n;//计算PI
}
void main()
{
	int n;
	double PI;
	printf("蒙特卡罗概率算法计算π:\n");
	printf("输入点的数量:\n");
	scanf("%d",&n);//输入撒点个数
	PI=MontePI(n);//计算PI
	printf("PI=%f\n",PI);//输出结果
}

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值