The authors have developed 2D, 3D and axisymmetric boundary-element-method BEM and BEM system based on these programs, in order to quantitatively estimate cathodic protection and macro-cell corrosion such as galvanic corrosion and differential-aeration-cell-corrosion. Because the knowledge of physical quantities potential and current density on the surface of corroding materials is of prime importance in corrosion problems, the authors have developed the BEM, which does not require internal elements discretization with internal element. Whereas, conventional analysis techniques, e.g. finite-difference-method FDM and finite-element-method FEM require internal element discretization which brings difficulties on the analysis of complicated 3D regions. In this paper, the history of development of corrosion analysis are reviewed briefly, Next, the BEM and BEM system are explained, and validity and usefulness of the BEM system are demonstrated by showing verification tests and application examples for cathodic protection.