Cross sections are homogenized over an entire node in nodal model implementation. The presence of a control rod (CR) partially inserted in the node has occasioned axial heterogeneity and generates a homogenization problem. If the homogenization process is only the volume-weighted average for nuclear parameters, the calculation of the multiplication factor and power in steady-state problems may mean relevant errors and for time-dependent problems may have caused the well-known cusping problem, which arises in three-dimensional transient simulations with CR motions. The major purpose of this technical note is to introduce an alternative method, based on the nodal expansion method, to deal with partially inserted CRs in nodes. One-dimensional equations, acquired through transverse integration of the neutron diffusion equation, have been modified to formulate the alternative method, which was evaluated in a transient problem. Furthermore, the alternative method gives satisfactory results and corrects the cusping effect in the case analyzed in this technical note.