% Produce entropy of uncertain variable from inverse uncertainty
% distribution
function f = entropy(inverse_uncertainty_distribution)
   y=@(alpha) inverse_uncertainty_distribution(alpha).*log(alpha./(1-alpha));
   f = quadl(y,0.0000001,0.9999999);
end