The objective of this paper is to describe a new method for modeling aeroelastic system in time domain based on a modification of the Laguerre Polynomials to represent complex quantities. These polynomials are used to approximate the unsteady aerodynamics forces which are defined in the frequency domain using the Doublet Lattice Method (DLM). In this approach, the size of the matrices representing the aeroelastic system remains the same as the matrices representing the structural dynamics behavior. It is an important point since classical state space aeroelastic models include lag states increasing the size of the matrices used to represent the system. The applicability of the method is demonstrated by numerical simulation performed on the benchmark wing structure. The approach offers promise mainly for complex systems such as real aircraft.ITA Aeronautical Institute of TechnologyUNESP - State University of São Paulo Fac. Eng. Ilha SolteiraUNESP - State University of São Paulo Fac. Eng. BauruUNESP - State University of São Paulo Fac. Eng. Ilha SolteiraUNESP - State University of São Paulo Fac. Eng. Bauru