In this paper, a methodology for accurate determination of multiple harmonic sources in a power system is presented. In the proposed methodology, assuming that no prior information regarding the possible locations of the harmonic sources in the system is available, the problem of determination of location of harmonic sources is solved in two stages. In the first stage, weighted least square estimation technique is used to find possible candidate buses (where the harmonic sources may be connected) among the system buses. In the second stage, using the concept of Euclidean norm, the exact buses (where the harmonic sources are present) among the possible candidate buses are obtained. To demonstrate the effectiveness of the proposed technique, the IEEE-14 and IEEE-30 bus test systems are used. (C) 2004 Elsevier Ltd. All rights reserved.