#include <cmath>

#include <iostream>

using std::cout;

using std::endl;

using std::log;

#include <boost/units/quantity.hpp>

#include <boost/units/io.hpp>

#include <boost/units/conversion.hpp>

namespace bu = boost::units;

using bu::quantity;

using bu::conversion_factor;

#include <boost/units/systems/si/prefixes.hpp>

namespace si = boost::units::si;

#include <boost/units/systems/information.hpp>

using namespace bu::information;

template <typename Sys>

quantity<bu::unit<bu::information_dimension, Sys> >

bernoulli_entropy(double p, const bu::unit<bu::information_dimension, Sys>&) {

typedef bu::unit<bu::information_dimension, Sys> requested_unit;

return quantity<requested_unit>((-(p*log(p) + (1-p)*log(1-p)))*nats);

}

int main(int argc, char** argv) {

}