calculates the sum of squares.
[y] = moc_sumsq(x [, DIM])
dimension
first DIMENSION, with more than 1 element
estimated standard deviation
Sum of squares of elements along dimension DIM. If DIM is omitted, it defaults to the first non-singleton dimension.