143       double probability = pdf(flip, successes);  in main()  local
145         << probability << " or 1 in " << 1. / probability  in main()
146         << ", or " << probability * 100. << "%" << endl;  in main()
155       double probability = cdf(flip, successes); // P(X <= heads)  in main()  local
157         << probability << " or 1 in " << 1. / probability << ", or "  in main()
158         << probability * 100. << "%"<< endl;  in main()
163     double probability = 0.3;  in main()  local
164     double q = quantile(flip, probability);  in main()
165 …std::cout << "Quantile (flip, " << probability << ") = " << q << std::endl; // Quantile (flip, 0.3…  in main()
166     probability = 0.6;  in main()
167     q = quantile(flip, probability);  in main()
168 …std::cout << "Quantile (flip, " << probability << ") = " << q << std::endl; // Quantile (flip, 0.6…  in main()