Compute discrete KL divergence and compare numerical versus closed-form KL for univariate Gaussian distributions.
