def CredibleInterval(pmf, percentage=<NUM_LIT>):
cdf = pmf.MakeCdf()<EOL>prob = (<NUM_LIT:1> - percentage / <NUM_LIT>) / <NUM_LIT:2><EOL>interval = cdf.Value(prob), cdf.Value(<NUM_LIT:1> - prob)<EOL>return interval<EOL>
Computes a credible interval for a given distribution.

    If percentage=90, computes the 90% CI.

    Args:
        pmf: Pmf object representing a posterior distribution
        percentage: float between 0 and 100

    Returns:
        sequence of two floats, low and high