simulator1.cpp

来自「圆周率统计计算方法仿真（即蒙特卡罗Monte Carlo仿真） 文件 应用」· C++ 代码 · 共 37 行

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

#define PI 3.141592653589793

void main( )
{

	 int total_num = 0;
	 int success_num = 0;
	 
	 double precision = 0.01 ;
	 double value1 ;
	 double value2 ;
	 double my_pi = 0;

	 srand( (unsigned)time( NULL ) ) ;

	 do
	 {
		     total_num++ ;
			 value1 = (double) rand()/RAND_MAX ;
			 value2 = (double) rand()/RAND_MAX ;
			 if( (value1 * value1 + value2 * value2) < 1 )
				 success_num++ ;			
					 
			 my_pi = 4*(double)success_num/(double)total_num ;

	 }while( (my_pi > ( PI + precision)) || (my_pi < (PI - precision)) );
	

	 printf("圆周率的值%f    落入1/4圆中的点数%d  单位正方形中的总点数%d\n",PI,success_num,total_num);
	 printf("圆周率的估计值为：%f\n所用点数：%d\n精度要求:%f\n",my_pi,total_num,precision);  

}