code: def _B(self, x, a, b):
          
          return special.betainc(a, b, x) * special.beta(a, b)
docs: incomplete Beta function as described in Mamon&Lokas A13
  
          :param x:
          :param a:
          :param b:
          :return: