When I tried to fit against the original (non-log) data, the only one that gave good results was fit = fitspline(x,y).
It’s just that spline doesn’t really tell me anything useful.
-------- Spline fit ---------------------------
x spline: x = [0.001, 0.011...
y spline: y = [1005.2007999999998, 90.34990000000002...
-----------------------------------------------