In this work, a nodal and iterative technique to evaluate the effective multiplication factor as well as the neutron flux, in multigroup diffusion problems, is presented. An iterative scheme, similar to the source iteration method, is implemented to decouple the system of differential equations which is the fundamental mathematical model. Then, analytical solutions are derived for the one-dimensional transverse integrated equations, of each energy group, resulting from a nodal approach. Constant approximations are assumed for the unknown transverse leakage terms in the contours of the nodes. In addition, constant and linear representations are investigated to express the fluxes in the source term to be updated in the iterative process. Numerical results for the effective multiplication factor were obtained for a series of two-dimensional multigroup problems with up-scattering and down-scattering. The procedure is simple, fast, the analysis of the results indicated a satisfactory agreement with results available in the literature and the use of different approximations to the source term seems to be a good alternative, instead of using higher-order approximations on the contour of the nodes, to improve accuracy.