Bro the discrepancy (2.8820544184511335e-19) is simply a result of floating point arithmetic, and not because of np.linspace or np.savez or even solve_ivp. Even when you use identical expressions in different parts of your code, the limited precision of floating point representation (and the way operations like division and multiplication are performed) can introduce such minuscule differences.