function [arith, geom, med] = meansAndMedian(v)
% [arith, geom, median] = meansAndMedian(v)

n = length(v);
arith = sum(v) / n;
geom  = prod(v).^(1/n);
mid = ceil(n/2);
if mod(n,2) == 0
    med = (v(mid) + v(mid+1)) /2;
else
    med = v(mid);
end