function d = KL(nu,mu)
mu = mu+realmin;
d = entropy(nu./mu,mu);
return