This article deals with the general linearization problem of Jacobi polynomials. We provide two approaches for finding closed analytical forms of the linearization coefficients of these polynomials. The first approach is built on establishing a new formula in which the moments of the shifted Jacobi polynomials are expressed in terms of other shifted Jacobi polynomials. The derived moments formula involves a hypergeometric function of the type 4F3(1), which cannot be summed in general, but for special choices of the involved parameters, it can be summed. The reduced moments formulas lead to establishing new linearization formulas of certain parameters of Jacobi polynomials. Another approach for obtaining other linearization formulas of some Jacobi polynomials depends on making use of the connection formulas between two different Jacobi polynomials. In the two suggested approaches, we utilize some standard reduction formulas for certain hypergeometric functions of the unit argument such as Watson’s and Chu-Vandermonde identities. Furthermore, some symbolic algebraic computations such as the algorithms of Zeilberger, Petkovsek and van Hoeij may be utilized for the same purpose. As an application of some of the derived linearization formulas, we propose a numerical algorithm to solve the non-linear Riccati differential equation based on the application of the spectral tau method.