{

cout << "Example 1 using Inverse Gamma function. " << endl;

#ifdef BOOST_NO_CXX11_NUMERIC_LIMITS

int max_digits10 = 2 + (boost::math::policies::digits<double, boost::math::policies::policy<> >() * 30103UL) / 100000UL;

cout << "BOOST_NO_CXX11_NUMERIC_LIMITS is defined" << endl;

#else

int max_digits10 = std::numeric_limits<double>::max_digits10;

#endif

cout << "Show all potentially significant decimal digits std::numeric_limits<double>::max_digits10 = "

<< max_digits10 << endl;

cout.precision(max_digits10); //

double x = 1.;

double a = 10;

double r = boost::math::gamma_q_inv(a ,x);

cout << " x = " << x << ", = gamma_q_inv(a,x)" << r << endl; //

return 0;