docs: def logspline_interpolate(x1, y1, x2):
      """
      Given a function at a set of points (x1, y1), interpolate to
      evaluate it at points x2.
      """
      sp = Spline(log(x1), log(y1))
      return exp(sp(log(x2)))