Abstract: The selection of numerical integration method of precise direct integration is discussed in this paper. Theoretical analysis and numerical results show that the feasibility of numerical integration method depends on the characteristics of loads. It is suggested that a numerical integration method of high algebraic accuracy be used if the load can be expressed in a polynomial form. It is shown that Cotes formula and Gauss formula are two effective integration methods for precise algorithm.