#include <iostream>
#include "math.h"
#include <TRandom3.h>
#include <iomanip>
#define PI 3.1415926535897932384626433832795
using namespace std;

int main() {
  TRandom3  gen(0);
  double random;
  double x;
  double LSinTh=0;
  double theta;
  double pi;
   int N;
  int Ns=0;
    double d;
  double l;
  double Prob;
  double SigmaProb;
  double SigmaPi;

  
  cout <<" Inserisci la larghezza d delle strisce (in cm)" << endl;
  cin >> d;
  
   cout <<" Inserisci la lunghezza dell'ago  (in cm)" << endl;
  cin >> l;
  if(l> d){
    cout<< "Ago piu' lungo della larghezza della striscia! Deve essere piu' corto!! " << endl;
    return 1;
    }
  cout<< "Inserisci il numero di punti Monte Carlo = " << endl;
  cin >> N;
  
 for (int i=0; i<N; i++){
    
    x=d*gen.Rndm();
	theta = PI*gen.Rndm();
	LSinTh = l*sin(theta);
     if(x < d && LSinTh > x) Ns++;
	 
 }

 Prob = (double)Ns/(double)N;
 pi = 2*l/(d*Prob);
 SigmaProb = sqrt(Ns -Ns*Ns/N)/N;
 SigmaPi = 2*l*SigmaProb/(d*Prob*Prob);
 
 cout << "La stima MC di pi greco   con "<< N << " punti  e' " << setprecision(6) << fixed<< pi << " +- " << SigmaPi << endl;
  return 0;
}