in fatal/math/statistical_moments.h [79:99]
void add(value_type const &sample) {
auto const n = samples_++;
auto const delta = sample - moment_1_;
auto const normalized_delta = delta / samples_;
auto const normalized_delta_squared = normalized_delta * normalized_delta;
auto const term1 = delta * normalized_delta * n;
moment_4_ += term1 * normalized_delta_squared
* (samples_ * samples_ - 3 * samples_ + 3)
+ 6 * normalized_delta_squared * moment_2_
- 4 * normalized_delta * moment_3_;
moment_3_ += term1 * normalized_delta * (n - 1)
- 3 * normalized_delta * moment_2_;
moment_2_ += term1;
moment_1_ += normalized_delta;
}