#include #include #include #include /* krävs för log()-funktionen */ /* Programmet simulerar tio poissonfördelade utfall där förväntade antal händelser är LAMBDA. */ #define LAMBDA ((float)1/6) main() { int i,antal; double intervall, tid; randomize(); /* initiera slumptalsgeneratorn */ for(i=0;i<10;i++){ tid=0; antal=0; while(tid < 1){ /* upprepa tills det gått en timme */ intervall = -log((double) rand()/((double)RAND_MAX))/LAMBDA; /* beräkna tiden till nästa händelse */ tid += intervall; antal++; /* antal händelser */ } printf("%d\n", antal-1); } }