The main purpose of this paper is to apply the boundary integral equation (BIE) method to the analysis of spoof localized surface plasmons (spoof LSPs) excited in a perfectly conducting cylinder with longitudinal corrugations. Frequency domain BIE schemes based on electric field integral equation (EFIE), magnetic field integral equation (MFIE) and combined field integral equation (CFIE) formulations are used to solve two-dimensional electromagnetic (EM) problems of scattering from the cylinder illuminated by a transverse electric plane wave. In this approach effects of spoof LSPs are included in the secondary surface current and charge densities resulting from the interaction between the plane wave and the cylinder. Numerical results obtained with the BIE schemes are validated by comparison with that of a recently proposed modal solution based on the metamaterial approximation.