#include #include #include int main() { using namespace boost::accumulators; accumulator_set > > acc; for(unsigned i = 125 ; i < 150; ++i) { acc(i); std::cerr << moment<4>(acc) << '\n'; } }