kld function

Kullback-Leibler Divergence