"should be around +1 " ... WRONG
"should be around -1 " ... WRONG
Here is the verdict, both Mathematica and Mathcad do likely under evaluate by some unknown value. Thiele is much closer, but that does not prove your spline is correct. If I would be asked to declare, then it would be Thiele. Not all 15 decimals were entered in Mathematica.
jmG