code: def compactspace(scale, n):
      r
      logit = logistic(scale=scale).ppf
      compact_xs = np.linspace(0, 1, n + 2)[1:-1]
      return logit(compact_xs)
docs: r"""
      Returns points :math:`x` spaced in the open interval
      :math:`(-\infty, \infty)`  by linearly spacing in the compactified
      coordinate :math:`s(x) = e^{-\alpha x} / (1 + e^{-\alpha x})^2`,
      where :math:`\alpha` is a scale factor.