A Computational Procedure for Obtaining the Poles of a Spherical Harmonic of order N; Application to the Multipole Expansion of Electrostatic Interaction Abstract: The solution for the poles of a general spherical surface harmonic is analyzed with application to the multipole calculation for the interaction of two charge distributions. The theoretical advantages are discussed and have been verified in test calculation using a FORTRAN program available on request.