function f = variance(inverse_uncertainty_distribution)
e = quadl(inverse_uncertainty_distribution,realmin,1-realmin);
b = 0; t = 1; small=0.00001;
while t-b>small
    c = (b+t)./2;
    if inverse_uncertainty_distribution(c) < e
        b=c;
    else t=c;
    end
end
phie = (b+t)./2;
f = 0; alpha = small;
while alpha<phie
    alpha = alpha+small;
    f = f+(e-inverse_uncertainty_distribution(alpha)).*alpha.*(inverse_uncertainty_distribution(alpha)...
        -inverse_uncertainty_distribution(alpha-small))
end
alpha = phie;
while alpha<1-small
    alpha = alpha+small;
    f = f+(inverse_uncertainty_distribution(alpha)-e).*(1-alpha).*(inverse_uncertainty_distribution(alpha)...
        -inverse_uncertainty_distribution(alpha-small))
end
f = 2.*f;
end