function y = gmean(x);
y = exp(mean(log(x)));
