Numerical Analysis, Centered Finite Difference in Python, Quadratic Rate of Convergence