gibbsNorm function

Fit a univariate normal (Gaussian) distribution to the observed data.