I need help making a program in python for Natural Cubic Spline, (the picture above shows the algorithm from Numercial Analysis Textbook)

Can anyone help create the code in python using the same exact method shown below? I tried following the steps but it's just not working

The output should be something like this:

(a2-a1)-ho (al-ao) x=1c1 and (an-1-an-2) Cn (an-an-1)- 0 The matrix A is strictly diagonally dominant, so it satisfies the hypotheses of Theore in Section 6.6. Therefore, the linear system has a unique solution for Co. c1, .. , cn. The solution to the cubic spline problem with the boundary conditions S"(n)- S(xn) = 0 can be obtained by applying Algorithm 3.4. Natural Cubic Spline To construct the cubic spline interpolant S for the function f, defined at the numbers ALGORITHM 3.4 to