SUMMARYThis paper discusses the calculation of potential distribution of impressed cathodic protection (CP) models with non-linear polarization curves. We propose a relaxation iterative algorithm for the nonlinear problem and prove both theoretically and numerically that this iterative sequence is convergent for any physical polarization curves. This feature is of signiÿcant importance in developing a computer code for the design of CP systems.