Nonlinear mode coupling of the (shear) Alfvén eigenmodes (AEs) is studied both analytically and numerically to understand the effect of the magnetohydrodynamic nonlinearity on the saturation mechanism. The energy transfer via the mode coupling can restrict the amplitude to a level that is estimated by the frequency mismatch and the coupling coefficient among the coupled modes. For the cases of cylindrical and toroidal geometries, new numerical codes are developed to calculate the coupling coefficient directly according to the Lagrangian theory of weakly nonlinear perturbations. It is shown that the shear AEs can couple with each other through the small compressional components of their eigenfunctions. By assuming low-β plasma, an analytic estimation of the coupling coefficient (∝ √ β) is derived in the cylindrical case. A coupling of reversed shear AEs in toroidal geometry indicates a more effective energy transfer at a lower amplitude level owing to the enhanced compressional perturbations in the poloidal sidebands. * a ,ξ b ,ξ c ) [9]. This formalism will be