calculates the sum of squares.
[y] = nan_sumsq(x [, DIM])
dimension
first DIMENSION, with more than 1 element
estimated standard deviation
features: - can deal with NaN's (missing values) - dimension argument also in Octave - compatible to Matlab and Octave