  A new weak Galerkin finite element method, called generalized weak Galerkin
method ({g}WG), is introduced for Stokes equations in this paper by using a new
definition of the weak gradient. Error estimates in energy norm and $L^2$ norm
for the velocity and $L^2$ norm for the pressure are derived for elements with
arbitrary combination of polynomials. Some numerical examples are presented to
verify the effectiveness, theoretical convergence orders, and robustness of the
proposed scheme.
