%%% Kullback-Leibler divergence between two Bernouilli   %%%
%%% distributions with respective probabilities a and b. %%%

function res=Kullback_Bernoulli(a,b)
res = a*log(a/b) + (1-a)*log((1-a)/(1-b));
end
