#include <stdio.h> #include <stdlib.h> #include <time.h> #include <unistd.h> int main() { srand48(time(NULL) ^ getpid()); for (double i = 0.0; i <= 1.0001; i += 0.05) { unsigned c[2] = { 0, 0 }; for (unsigned j = 0; j < 10000000; ++j) { for (double x = i;;) { double d = x < 0.5 ? x : 1.0 - x; x = drand48() < 0.5 ? x - d : x + d; if (x <= 0.0 || x >= 1.0) { c[x <= 0.0 ? 0 : 1]++; break; } } } printf("x=%.3f p=%.4f\n", i, (double)c[1] / (double)(c[0] + c[1])); } }