72 long double factorial(
long double);
73 long double power_of_4(
long double);
80 long double n = 0.0, i=0.0;
81 long double prevsum = 0.0, sum = 0.0, prevsumdiff = 0.0, sumdiff = 0.0, term1 = 0.0;
82 for(n = 0.0; n <= 30000.0; n++)
84 term1 = factorial(2*n) / power_of_4(n);
85 for (i=1.0; i <= n ; i++)
86 sum = sum + (1.0 / (factorial(n+i) * factorial(n-i)));
88 cout <<
"Probability of good choice for population of " << 2*n <<
"=" << sum*100.0 <<endl;
89 sumdiff = sum - prevsum;
90 cout <<
"prob - prevprob = " << sumdiff << endl;
91 cout <<
"Convergence test: (sum - prevsum)/prevsum = " << sumdiff/prevsum << endl;
93 prevsumdiff = sumdiff;
100 long double factorial(
long double n)
105 return (
long double) n*factorial(n-1);
108 long double power_of_4(
long double n)
110 long double power = 1.0 ;
112 for (i=n;i > 0.0;i--)