code: def logit(x, a=0., b=1.):
          r
          return numpy.log(x-a) - numpy.log(b-x)
docs: r"""Computes the logit function with domain :math:`x \in (a, b)`.
  
          This is given by:
  
          .. math::
  
              \mathrm{logit}(x; a, b) = \log\left(\frac{x-a}{b-x}\right).
  
          Note that this is also the inverse of the logistic function with range
          :math:`(a, b)`.
  
          Parameters
          ----------
          x : float
              The value to evaluate.
          a : float, optional
              The minimum bound of the domain of x. Default is 0.
          b : float, optional
              The maximum bound of the domain of x. Default is 1.
  
          Returns
          -------
          float
              The logit of x.