code: def MakeCdf(self, steps=101):
          
          xs = [i / (steps - 1.0) for i in xrange(steps)]
          ps = [scipy.special.betainc(self.alpha, self.beta, x) for x in xs]
          cdf = Cdf(xs, ps)
          return cdf
docs: Returns the CDF of this distribution.