An effective homogenization of control rods is studied by considering neutron leakage as a fixed boundary source. Good numerical convergence is seen in the iterative calculation of effective homogenized cross sections that preserve integrated reaction rates. Usefulness of the present method is demonstrated by treating off-center control rods in a one-dimensional slab model of an FBR.