public static Complex CentralMoments(double[] angles, int order)
{
double theta = Mean(angles);
double v = angles.Length / theta;
double[] alpha = new double[angles.Length];
for (int i = 0; i < alpha.Length; i++)
alpha[i] = Distance(angles[i], v);
return NoncentralMoments(angles, order);
}