Wyznaczanie wielkości liczby PI metodą Monte Carlo

 #include <iostream>

#include <cstdlib>

#include <cmath>

#include <ctime>

using namespace std;

int main()

{

double a,b; //zmienne pomocnicze - współrzędne losowanego punktu

long  T = 0;//ilość wylosowanych punktów, które znajdują się w kole

long R;     //ilość wszystkich wylosowanych punktów znajdujących się w kwadracie 

long ile;

srand(time(NULL));

cout<<"Na podstawie ilu punktow mam znalezc wartosc liczby pi: ";

cin>>ile;

for (R=0; R<ile; R++)

{

a=(double)rand()/(RAND_MAX);

b=(double)rand()/(RAND_MAX);

if (a*a*b*b<=1)

T++;

}

cout<<"liczba pi ma wartosc "<<(double)(4*T)/R<<endl;

return 0;

}