#include <boost/math/tools/roots.hpp> // root finding.

#include <boost/math/distributions/normal.hpp> // for normal_distribution

using boost::math::normal; // typedef provides default type is double.

#include <iostream>

using std::cout; using std::endl; using std::left; using std::showpoint; using std::noshowpoint;

#include <iomanip>

using std::setw; using std::setprecision;

#include <limits>

using std::numeric_limits;

#include <stdexcept>

int main()

} // int main()