/* 利用蒙特卡洛算法近似求圆周率Pi*/
/*程序使用:VC++6.0 */
#include
#include
#include
#define COUNT 800
/*循环取样次数,每次取样范围依次变大*/
void main()
{ double x,y; int num=0; int i;
for(i=0;i
/*RAND_MAX=32767,包含在
y=rand()*1.0/RAND_MAX;
if((x*x+y*y)<=1) num++;
/*统计落在四分之一圆之内的点数*/
}
printf("Pi值等于:%f\n",num*4.0/COUNT);
}
望采纳
程序:(该算法可以修改后用Mathematica计算或者Matlab)
/* 利用蒙特卡洛算法近似求圆周率Pi*/
/*程序使用:VC++6.0 */
#include
#include
#include
#define COUNT 800
/*循环取样次数,每次取样范围依次变大*/
void main()
{ double x,y; int num=0; int i;
for(i=0;i
/*RAND_MAX=32767,包含在
y=rand()*1.0/RAND_MAX;
if((x*x+y*y)<=1) num++;
/*统计落在四分之一圆之内的点数*/
}
printf("Pi值等于:%f\n",num*4.0/COUNT);
}