cum3(a, b = a, c = a, unbiased = TRUE)a.
If supplied then it must be the same length as a.a.
If supplied then it must be the same length as a.n/((n-1)*(n-2)) where n is
the sample size, if unbiased is FALSE then a multiplier of 1/n is used.
This is multiplied by sum((a-mean(a))*(b-mean(b))*(c-mean(c))) to give the
required estimate.