Implementing the online algorithm used here:https://www.johndcook.com/blog/2013/06/11/computing-skewness-and-kurtosis-in-one-pass/ to compute skewness, kurtosis and bimodality of a dataset represented here by an array.

Bimodality Coefficient Difference calculator