Files
SCR/SCR3.1/TP4/Exo2/montecarlo.c
2025-12-05 13:48:10 +01:00

19 lines
317 B
C

uint64_t shots_in=0,shots=0;
double x,y;
int main(int argc, char **argv)
{
shots_in = 0;
shots = 0;
for (;;){
x = rand()/(RAND_MAX*1.0);
y = rand()/(RAND_MAX*1.0);
shots ++;
if ((x*x+y*y)<=1){
shots_in ++;
}
}
/* la probabilité vaut tirsIn/tirs,
Elle converge lentement vers pi/4*/
}