在一个2*2的正方形内投掷飞镖，根据飞镖落在单位圆中的数量，计算PI(π)的值

编程实现如下模拟：

飞镖被随机地投掷到以（1，1）和（-1，-1）为对角的正方形内。若飞镖落在单位圆内【即以（0，0）为圆心，1为半径的圆】，算命中。否则未命中。

运行这个模拟并用它求PI的近似值。

 

模型如下图所示：

程序如下：

// Circle_PI.cpp : 定义控制台应用程序的入口点。//#include "stdafx.h"#include <iostream>#include <cstdlib>#include <ctime>#include <iomanip>using namespace std;const int MAX_DAST=100000;//投掷总数double rand_double(double a,double b);void rand_seed();int main(void){int PI_DAST=0;//投掷在圆里面的总数rand_seed();for(int i=0;i<MAX_DAST;i++){double x=rand_double(-1,1);double y=rand_double(-1,1);if( x*x+y*y<=1)PI_DAST++;}double PI=PI_DAST*4.0/MAX_DAST;cout<<fixed;cout<<"PI="<<PI<<endl;system("pause");return 0;}/**生成某一区域内的随机浮点数@param a 该区域的下边界@param b 该区域的上边界@param return 随机浮点数x，x∈[a，b]*/double rand_double(double a,double b){return a+(b-a)*rand()*1.0/RAND_MAX;}/**设置随机数生成器的种子*/void rand_seed(){int seed=static_cast<int>(time(0));srand(seed);}

结果截图：