Le C en 7 jours - L'algorithme de Monte-Carlo

13

L'algorithme de Monte-Carlo

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
int main( ) {
const long NTHROWS = 1000000; // how many times do we throw the dart
unsigned long inside = 0;
unsigned seed = time( 0 );
double x, y;
unsigned long i;
srand( seed );
for ( i = 0; i < NTHROWS; i++ ) {
x = (double)rand()/RAND_MAX;
y = (double)rand()/RAND_MAX;
if ( x*x + y*y < 1.0 )
inside++;
}
printf("PI=%f\n", 4.0*inside/NTHROWS);
}

$ gcc -o mc mc.c
$ ./mc
PI=3.139084

Réglez NTHROWS à 10000000 pour obtenir une valeur plus proche (et tester la vitesse de votre CPU) :

$ gcc -o mc mc.c
$ ./mc
PI=3.141202

Commentaires