The complex aeroelastic equation of an aircraft wing with freeplay nonlinearity was built and applied for aeroelastic response analysis. This took into account multiple lumped mass points and directions of the control surface using the Roger rational function approximation. The aerodynamic analysis was based on the subsonic doublet lattice method and the external program of time history analysis method was employed. The aeroelastic response results show that amplitudes of generalized displacement considering the freeplay effect is decreased compared with those not considering the freeplay effect, and influenced by freeplay stiffness coefficients. Constant and variable step length calculations of the fourth order Runge-Kutta method can influence the results and efficiency.