10

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

Comments

Your comment:
[p] [b] [i] [u] [s] [quote] [pre] [br] [code] [url] [email] strip help 2000

Enter a maximum of 2000 characters.
Improve the presentation of your text with the following formatting tags:
[p]paragraph[/p], [b]bold[/b], [i]italics[/i], [u]underline[/u], [s]strike[/s], [quote]citation[/quote], [pre]as is[/pre], [br]line break,
[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]command[/code], [code=language]source code in c, java, php, html, javascript, xml, css, sql, bash, dos, make, etc.[/code].