The electrostatics of two cylinders charged to the symmetrical or anti-symmetrical potential is investigated by using the null-field boundary integral equation (BIE) in conjunction with the degenerate kernel of the bipolar coordinates. The undetermined coefficient is obtained according to the Fredholm alternative theorem. The uniqueness of solution, infinite solution, and no solution are examined therein. A single cylinder (circle or ellipse) is also provided for comparison. The link to the general solution space is also done. The condition at infinity is also correspondingly examined. The flux equilibrium along circular boundaries and the infinite boundary is also checked as well as the contribution of the boundary integral (single and double layer potential) at infinity in the BIE is addressed. Ordinary and degenerate scales in the BIE are both discussed. Furthermore, the solution space represented by the BIE is explained after comparing it with the general solution. The present finding is compared to those of Darevski[2] and Lekner[4] for identity.