function [s] = sst(x)
m = mean(x);
xx = x - m;
s = sum(xx .* xx);