68     binomial flip(flips, success_fraction);  in main()  local
75       << mean(flip) << " heads (or tails)." << endl;  in main()
76     cout << "Mode is " << mode(flip) << endl;  in main()
77     cout << "Standard deviation is " << standard_deviation(flip) << endl;  in main()
79       <<  ceil(mean(flip) - standard_deviation(flip))  << " and "  in main()
80       << floor(mean(flip) + standard_deviation(flip)) << " correct." << endl;  in main()
81     cout << "Skewness is " << skewness(flip) << endl;  in main()
85     cout << "Skewness if success_fraction is " << flip.success_fraction()  in main()
86       << " is " << skewness(flip) << endl << endl; // Expect zero for a 'fair' coin.  in main()
90     cout << "For " << flip.trials() << " coin flips: " << endl;  in main()
91     cout << "Probability of getting no heads is " << pdf(flip, 0) << endl;  in main()
92     cout << "Probability of getting at least one head is " << 1. - pdf(flip, 0) << endl;  in main()
97       << pdf(flip, 0) + pdf(flip, 1) << endl; // sum of exactly == probabilities  in main()
101     cout << "Probability of getting 0 or 1 (<= 1) heads is " << cdf(flip, 1) << endl;  in main()
102     cout << "Probability of getting 9 or 10 heads is " << pdf(flip, 9) + pdf(flip, 10) << endl;  in main()
106     cout << "Probability of getting 9 or 10 heads is " << 1. - cdf(flip, 8) << endl;  in main()
110     cout << "Probability of getting 9 or 10 heads is " << cdf(complement(flip, 8)) << endl;  in main()
121       << pdf(flip, 4) + pdf(flip, 5) + pdf(flip, 6) << endl;  in main()
127       << cdf(flip, 6) - cdf(flip, 3) << endl;  in main()
133       << cdf(flip, 7) - cdf(flip, 2) << endl;  in main()
143       double probability = pdf(flip, successes);  in main()
155       double probability = cdf(flip, successes); // P(X <= heads)  in main()
164     double q = quantile(flip, probability);  in main()
167     q = quantile(flip, probability);  in main()