4

Monte-Carlo

Rand  |  Prime

mc.cpp
  1. #include <iostream>
  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.     double x, y;
  12.  
  13.     srand( seed );
  14.  
  15.     for ( unsigned long i = 0; i < NTHROWS; i++ ) {
  16.         x = float(rand())/RAND_MAX;
  17.         y = float(rand())/RAND_MAX;
  18.         if ( x*x + y*y < 1 )
  19.             inside++;
  20.     }
  21.     std::cout << "PI=" << 4.0*inside/NTHROWS << std::endl;
  22. }

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].