Geometric mean and variance
[m,v] = geostat(p)
[m,v] = geostat(p) returns
the mean m and variance v of
a geometric distribution with corresponding probability parameters
in p. p can be a vector, a matrix,
or a multidimensional array. The parameters in p must
lie in the interval [0,1].