{

cout << "Example using Inverse Gamma distribution. " << endl;

// TODO - awaiting a real example using Bayesian statistics.

#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 shape = 1.;

double scale = 1.;

double x = 0.5;

// Construction using default RealType double, and default shape and scale..

inverse_gamma_distribution<> my_inverse_gamma(shape, scale); // (alpha, beta)

cout << "my_inverse_gamma.shape() = " << my_inverse_gamma.shape()

<< ", scale = "<< my_inverse_gamma.scale() << endl;

cout << "x = " << x << ", pdf = " << pdf(my_inverse_gamma, x)

<< ", cdf = " << cdf(my_inverse_gamma, x) << endl;

// Construct using typedef and default shape and scale parameters.

inverse_gamma my_ig;

inverse_gamma my_ig23(2, 3);

cout << "my_inverse_gamma.shape() = " << my_ig23.shape()

<< ", scale = "<< my_ig23.scale() << endl;

cout << "x = " << x << ", pdf = " << pdf(my_ig23, x)

<< ", cdf = " << cdf(my_ig23, x) << endl;

// Example of providing an 'out of domain' or 'bad' parameter,

// here a shape < 1, for which mean is not defined.

// Try block is essential to catch the exception message.

// (Uses the default policy which is to throw on all errors).

try

{

inverse_gamma if051(0.5, 1);

//inverse_gamma if051(0.5, 1);

cout << "mean(if051) = " << mean(if051) << endl;

}

catch(const std::exception& e)

{ // Always useful to include try & catch blocks because default policies

// are to throw exceptions on arguments that cause errors like underflow, overflow.

// Lacking try & catch blocks, the program will abort without a message below,

// which may give some helpful clues as to the cause of the exception.

std::cout <<

"\n""Message from thrown exception was:\n " << e.what() << std::endl;

}

return 0;