KL_div function

Computes the Kullback-Leibler divergence based on kernel density estimates