[[58136.0, 6045.8897],[58137.0, 4221.8818],[58142.0, 4021.8887],[58143.0, 9585.4870]]
This data set (in the form [[x0, y0], [x1, y1], ... ]) causes cubicSpline $ tabulate xs ys to fail with non-invertible matrix error.
This works as expected in octave with the command interp1(xs, ys, xs, "cubic")