Transfer function method was employed to study the flutter of highaspectratio flexible aircraft wings carrying external stores. Firstly, the flexible wing was modeled as a curve beam, and the flutter differential equations of a curvebeam element were established. Further, by using the transfer function method, the differential equations were formulated in a statespace form. Moreover, considering the deformation compatibility and internal force balance at the position of the external store, the governing equations of the whole wing were obtained by element assembling. ,Both the flutter velocity and flutter frequency were obtained by solving complex eigenvalue equations. The results are good agreement with finite element analysis, which indicates that the present method is accurate and efficient for flutter analysis of the flexible wing. Finally, the effects of the mass, moment inertia, number and position of the external store on the aeroelastic characteristics of the flexible wing were investigated.