Function that computes the four parameters required for the GLD reconstruction problem.
Computes the centered mean, variance, skewness and kurtosis from the moments. Prints error message for insufficient moments.