public static Complex NoncentralMoments(double[] angles, int order)
{
double cbar = 0;
double sbar = 0;
for (int i = 0; i < angles.Length; i++)
{
cbar += Math.Cos(angles[i] * order);
sbar += Math.Sin(angles[i] * order);
}
return new Complex(cbar, sbar);
}
}