In this paper we analyze the numerical aspects of the various methods that have been utilized to analyze thin wire antennas. First, we derive the properties of the operators for Pocklington's and Hallen's integral equation. On the basis of these properties, we discuss the various iterative methods used to find current distribution on thin wire structures. An attempt has been made to resolve the question of numerical stability associated with various entire domain and subdomain expansion functions in Galerkin's method. It has been shown that the sequence of solutions generated by the iterative methods monotonically approaches the exact solution provided the excitations chosen for these problems are in the range of the operator. Such a statement may not hold for Galerkin's methods if the inverse operator is unbounded. Moreover, if the excitation function is not in the range of the operator, then the sequence of solutions forms an asymptotic series. Examples have been presented to illustrate this point.