geomean
              Compute the geometric mean of x.
geomean(x) returns the
 geometric mean of the elements in x defined as
 $$ {\rm geomean}(x) = \left( \prod_{i=1}^N x_i \right)^\frac{1}{N}
 = exp \left({1\over N} \sum_{i=1}^N log x_i \right) $$
   where  is the length of the x vector.
 geomean(x) returns a row
 vector with the geometric mean of each columns in x.
 geomean(x)
 operates along the first nonsingleton dimension of x.
  geomean(x, "all") returns the geometric mean of all the elements
 in x.  If x contains any 0, then the returned value is 0.
 geomean(x, dim) returns the geometric mean along the
 operating dimension dim of x.  Calculating the harmonic mean of
 any subarray containing any 0 will return 0.
 geomean(x, vecdim) returns the geometric mean over the
 dimensions specified in the vector vecdim.  For example, if x is
 a 2-by-3-by-4 array, then geomean(x, [1 2]) returns a
 1-by-1-by-4 array.  Each element of the output array is the geometric mean of
 the elements on the corresponding page of x.  If vecdim indexes
 all dimensions of x, then it is equivalent to geomean (x,
 "all").  Any dimension in vecdim greater than ndims (x)
 is ignored.
 geomean(…, nanflag) specifies whether to exclude NaN
 values from the calculation, using any of the input argument combinations in
 previous syntaxes. By default, geomean includes NaN values in the calculation
 (nanflag has the value "includenan").  To exclude NaN values, set the
 value of nanflag to "omitnan".
Source Code: geomean