(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);//输出结果
}
2197

被折叠的 条评论
为什么被折叠?



