def get_coverage(p05, p95, truth):
    """How often is the truth between p05 and p95."""
    ncovered = 0
    ntot = 0
    for yval, lower, upper in zip(truth, p05, p95):
        ntot += 1
        if lower <= yval <= upper:
            ncovered += 1
    return ncovered/ntot