13

L'algorithme de Monte-Carlo

  1. #include <stdio.h>
  2.  
  3. #include <stdlib.h>
  4. #include <time.h>
  5.  
  6. int main( ) {
  7.     const long NTHROWS = 1000000;   // how many times do we throw the dart
  8.  
  9.     unsigned long inside = 0;
  10.     unsigned seed = time( 0 );
  11.  
  12.     double x, y;
  13.     unsigned long i;
  14.    
  15.     srand( seed );
  16.  
  17.     for ( i = 0; i < NTHROWS; i++ ) {
  18.         x = (double)rand()/RAND_MAX;
  19.         y = (double)rand()/RAND_MAX;
  20.         if ( x*x + y*y < 1.0 )
  21.             inside++;
  22.     }
  23.     printf("PI=%f\n", 4.0*inside/NTHROWS);
  24. }
$ 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

Votre commentaire :
[p] [b] [i] [u] [s] [quote] [pre] [br] [code] [url] [email] strip aide 2000

Entrez un maximum de 2000 caractères.
Améliorez la présentation de votre texte avec les balises de formatage suivantes :
[p]paragraphe[/p], [b]gras[/b], [i]italique[/i], [u]souligné[/u], [s]barré[/s], [quote]citation[/quote], [pre]tel quel[/pre], [br]à la ligne,
[url]http://www.izend.org[/url], [url=http://www.izend.org]site[/url], [email]izend@izend.org[/email], [email=izend@izend.org]izend[/email],
[code]commande[/code], [code=langage]code source en c, java, php, html, javascript, xml, css, sql, bash, dos, make, etc.[/code].