The primary focus of this article is on applying specific generalized Jacobi polynomials (GJPs) as basis functions to obtain the solution of linear and non-linear even-order two-point BVPs. These GJPs are orthogonal polynomials that are expressed as Legendre polynomial combinations. The linear even-order BVPs are treated using the Petrov–Galerkin method. In addition, a formula for the first-order derivative of these polynomials is expressed in terms of their original ones. This relation is the key to constructing an operational matrix of the GJPs that can be used to treat the non-linear two-point BVPs. In fact, a numerical approach is proposed using this operational matrix of derivatives to convert the non-linear differential equations into effectively solvable non-linear systems of equations. The convergence of the proposed generalized Jacobi expansion is investigated. To show the precision and viability of our suggested algorithms, some examples are given.