Two semi-analytic solutions for a perturbed two-body problem known as Lagrange planetary equations (LPE) were compared to a numerical integration of the equation of motion with same perturbation force. To avoid the critical conditions inherited from the configuration of LPE, non-singular orbital elements (EOE) had been introduced. In this study, two types of orbital elements, classical Keplerian orbital elements (COE) and EOE were used for the solution of the LPE. The effectiveness of EOE and the discrepancy between EOE and COE were investigated by using several near critical conditions. The near one revolution, one day, and seven days evolutions of each orbital element described in LPE with COE and EOE were analyzed by comparing it with the directly converted orbital elements from the numerically integrated state vector in Cartesian coordinate. As a result, LPE with EOE has an advantage in long term calculation over LPE with COE in case of relatively small eccentricity.