code: def Prob(self, x):
          
          if x < self.xs[0]: return 0.0
          index = bisect.bisect(self.xs, x)
          p = self.ps[index - 1]
          return p
docs: Returns CDF(x), the probability that corresponds to value x.
  
          Args:
              x: number
  
          Returns:
              float probability