直接上代码:
#include <iostream>
#include <ctime>
using namespace std;
void main()
{
srand(unsigned(time(NULL)));
int total = 10000000;
float k, j;
int num = 0;
for (int i = 0; i < total; i++)
{
k = 1.0*rand() / RAND_MAX;
j = 1.0*rand() / RAND_MAX;
if (k*k + j*j <= 1)
num++;
}
cout << "Pi = " <<4.0* num / total << endl;
system("pause");
}
蒙特卡洛探讨待续: