One of the less-known integration methods is the weighted Newton–Cotes quadrature rule of semi-open type, which is denoted by ∫ a = x 0 b = x n + 1 = x 0 + ( n + 1 ) h f ( x ) w ( x ) d x ≃ ∑ k = 0 n w k f ( x 0 + kh ) , where w( x) is a weight function on [ a, b] and h = b - a n + 1 is a positive value. There are various cases for w( x) that one can use. Because of the special importance of the weight function of Gauss–Chebyshev quadrature rules in the numerical analysis, i.e. w ( x ) = 1 1 - x 2 , we consider this function as the main weight. Hence, in this paper, we face with the following formula: ∫ - 1 + 1 f ( x ) 1 - x 2 d x ≃ ∑ k = 0 n w k f - 1 + 2 k n + 1 , which has the precision degree is n + 1 for even n’s and n for odd n’s. In this paper, we consider bounds of above integration formula as two additional variables to reach a nonlinear system that numerically improves the precision degree up to n + 2. In this way, sevral examples are given to show the numerical superiority of our approach.