kl.dist function

Kullback-Leibler distance