In this work, by using Orlicz Spaces and variational methods, we can find enough conditions for the existence of least energy solution for a class of quasilinear systems of the formwhere V 1 , V 2 : R 2 → R are positive continuous potentials and h, g : R 2 × R × R → R are continuous functions with critical exponential growth in the sense of Trudinger-Moser inequality.