code: def MaxLikeInterval(self, percentage=90):
          
          interval = []
          total = 0
  
          t = [(prob, val) for val, prob in self.Items()]
          t.sort(reverse=True)
  
          for prob, val in t:
              interval.append(val)
              total += prob
              if total >= percentage / 100.0:
                  break
  
          return interval
docs: Returns the maximum-likelihood credible interval.
  
          If percentage=90, computes a 90% CI containing the values
          with the highest likelihoods.
  
          percentage: float between 0 and 100
  
          Returns: list of values from the suite