{

std::cout << "Examples of Arcsine distribution." << std::endl;

std::cout.precision(3); // Avoid uninformative decimal digits.

using boost::math::arcsine;

arcsine as; // Construct a default `double` standard [0, 1] arcsine distribution.

std::cout << pdf(as, 1. / 2) << std::endl; // 0.637

// pdf has a minimum at x = 0.5

std::cout << pdf(as, 1. / 4) << std::endl; // 0.735

std::cout << cdf(as, 0.05) << std::endl; // 0.144

std::cout << 2 * cdf(as, 1 - 0.975) << std::endl; // 0.202

std::cout << 2 * cdf(complement(as, 0.975)) << std::endl; // 0.202

std::cout << quantile(as, 1 - 0.2 / 2) << std::endl; // 0.976

std::cout << quantile(complement(as, 0.2 / 2)) << std::endl; // 0.976

{

using boost::math::arcsine_distribution;

arcsine_distribution<> as(2, 5); // Constructs a double arcsine distribution.

BOOST_MATH_ASSERT(as.x_min() == 2.); // as.x_min() returns 2.

BOOST_MATH_ASSERT(as.x_max() == 5.); // as.x_max() returns 5.

}

return 0;