蒙特卡洛法计算圆周率近似值(C语言)

C语言实现蒙特卡洛法计算圆周率近似值

随机生成ALL个范围0-1的数来模拟1×1的正方形，计算位于四分之一圆（X×X+Y×Y<1）内的个数inround,PI＝4×ALL/inround。其中宏的参数(ALL和RAND_MAX)对于精确度非常关键。

#include<time.h>#include<stdio.h>#include<stdlib.h>#define ALL 100000int main (){srand(time(0));double X[ALL],Y[ALL];int inround=0,index;for(index=0;index<ALL;index++){X[index]=(double)rand()/(double)RAND_MAX;Y[index]=(double)rand()/(double)RAND_MAX;if(X[index]*X[index]+Y[index]*Y[index]<=1)inround++;/*printf("%lf %lf\n",X[index],Y[index]); */}//计算平均数一定程度上检验生成的随机数是否合理// for(index=1;index<ALL;index++){X[0]+=X[index];Y[0]+=Y[index];}printf("\n%lf%lf",X[0]/ALL,Y[0]/ALL);//越接近0.5越合理// printf("\n圆周率近似= %lf\n",4*(double)inround/(double)ALL); return 0;}