In this paper, the state space method is used to analyze the static and dynamic behaviors of laminated magneto-electro-elastic rectangular plates with simply supported boundary conditions. Multilayers are considered as well as viscoelastic interfaces. The Kelvin-Voigt model is used to take into accounted the viscoelastic interface effects. The mathematical formulation is elaborated in a general form and an arbitrary number of layers as well as the orthotropic behavior can be considered. The coupling time and thickness variables formalism is presented in an elegant matrix way allowing a versatility and capability of the formalism to handle more general multilayered plates. The solution procedure is based on the transfer relationships coupled with the Lagrange polynomials in discretised time intervals. The static behavior of multilayered magneto-electro-elastic rectangular plates is analyzed and the main results are the viscoelastic interface effect on the dynamic behavior of multifunctional structures. New results have been obtained by this coupling state space-Lagrange polynomial methodology. The mathematical procedure, elaborated in this paper, is original and allows the numerical study of the effects of the viscoelastic interfaces. The effectiveness of the proposed methods has been demonstrated by performing variant numerical tests.