Comentarios en: Hallando Pi gracias a la probabilidad http://2geek2curious.com/2008/07/hallando-pi-gracias-a-la-probabilidad/ Éste es un blog creado por unos cuantos amigos un tanto frikis (con la acepción equivalente a geek) para hablar de temas varios que nos llaman la atención, como pueden ser la electrónica, tecnología en general, cocina, deportes, viajes o cualquier cosa que se nos ocurra. Sun, 29 Jan 2012 21:41:07 +0000 hourly 1 http://wordpress.org/?v=3.3.1 Por: Oscar http://2geek2curious.com/2008/07/hallando-pi-gracias-a-la-probabilidad/comment-page-1/#comment-30 Oscar Wed, 09 Jul 2008 09:44:24 +0000 http://2geek2curious.com/?p=50#comment-30 en la Escuela lo hacemos con un excel y sale muy bien en la Escuela lo hacemos con un excel y sale muy bien

]]> Por: luis http://2geek2curious.com/2008/07/hallando-pi-gracias-a-la-probabilidad/comment-page-1/#comment-27 luis Mon, 07 Jul 2008 09:36:29 +0000 http://2geek2curious.com/?p=50#comment-27 Me ha gustado, pero como es simulación va muy lenta he probado a hacerlo en c: <tt> #include <stdio.h> #include <math.h> #include <stdlib.h> #include <time.h> int main(void) {   double x, y, dist, pi;   unsigned long n_in = 0, n_out = 0, n = 0;   int i;   srand(time(NULL));   while(1) {     for(i=0; i<1000000; ++i) {       x = (double)(rand())/RAND_MAX - 0.5;       y = (double)(rand())/RAND_MAX - 0.5;       dist = sqrt(x*x+y*y);       if ( dist <= 0.5) ++n_in;       else ++n_out;       ++n;     }     pi = 4 * (double)n_in / n;     printf("Total: %6u\t\tIn: %6u\t\tOut: %6u\t\tPi: %1.6f\n", n, n_in, n_out, pi);   } } </tt> Después de 248 millones de puntos, la estimación de Pi seguía variando mucho, pero su útimo valor era de 3.141498 Me ha gustado, pero como es simulación va muy lenta he probado a hacerlo en c:


#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <time.h>

int main(void)
{
  double x, y, dist, pi;
  unsigned long n_in = 0, n_out = 0, n = 0;
  int i;

  srand(time(NULL));

  while(1) {
    for(i=0; i<1000000; ++i) {
      x = (double)(rand())/RAND_MAX - 0.5;
      y = (double)(rand())/RAND_MAX - 0.5;
      dist = sqrt(x*x+y*y);

      if ( dist < = 0.5) ++n_in;
      else ++n_out;
      ++n;
    }
    pi = 4 * (double)n_in / n;
    printf("Total: %6u\t\tIn: %6u\t\tOut: %6u\t\tPi: %1.6f\n", n, n_in, n_out, pi);
  }
}

Después de 248 millones de puntos, la estimación de Pi seguía variando mucho, pero su útimo valor era de 3.141498

]]>