code: def laplace_like(x, mu, tau):
      R
  
      return flib.gamma(np.abs(np.array(x) - mu), 1, tau) - \
          np.size(x) * np.log(2)
docs: R"""
      Laplace (double exponential) log-likelihood.
  
      The Laplace (or double exponential) distribution describes the
      difference between two independent, identically distributed exponential
      events. It is often used as a heavier-tailed alternative to the normal.
  
      .. math::
          f(x \mid \mu, \tau) = \frac{\tau}{2}e^{-\tau |x-\mu|}
  
      :Parameters:
        - `x` : :math:`-\infty < x < \infty`
        - `mu` : Location parameter :math:`-\infty < mu < \infty`
        - `tau` : Scale parameter :math:`\tau > 0`
  
      .. note::
        - :math:`E(X) = \mu`
        - :math:`Var(X) = \frac{2}{\tau^2}`