This work presents three efficient higher order analytical nodal methods for the numerical solution of a two-dimensional multigroup neutron diffusion equation in Cartesian geometry based on the use of the successive polynomial-weighted transverse integrations technique to convert a one-group diffusion equation to a system of coupled one-dimensional ordinary differential equations. These equations are then solved analytically over each homogenized cell after adequate approximations of the resulting effective sources after transversal integrations. Coupling between the approximate transverse flux-moments is achieved by imposing uniqueness constraint on their moments values. Adjacent elements are coupled by enforcing continuity conditions on the flux and current moments at interfaces cells. The weighted cell-balance equations and current-continuity conditions are then used to derive the discrete equations. These methods are applied for solving numerically various 2D benchmark problems and theirs performances discussed. Numerical results demonstrates more efficiency for the third higher order analytical nodal method for which the alone unknowns considered are the transverse flux moments on the interfaces of the homogenized elements.