docs: public double getSkewness() {
  		double skewness = NaN;
  		if (_n >= 3L) {
  			final double var = _m2.value()/(_n - 1.0);
  			if (var < 10E-20) {
  				skewness = 0.0d;
  			} else {
  				skewness = (_n*_m3.value())/
  						((_n - 1.0)*(_n - 2.0)*sqrt(var)*var);
  			}
  		}
  
  		return skewness;
  	}