Mathematical modeling of scientific and engineering processes often yield Boundary Value Problems (BVPs). One of the broad categories of numerical methods for solving BVPs is the finite difference methods, in which the differential equation is replaced by a set of difference equations which are solved by direct or iterative methods. In this paper, we use some properties of matrices to analyze the stability and convergence of the prominent finite difference methods - two-step Obrechkoff method - for solving the boundary value problem $u^{\prime \prime} = f(t,u)$, $a < x < b$, $u(a) = \eta_1$, $u(b) = \eta_2$. Conditions for the stability and convergence of the two-step Obrechkoff method method were established.