In this papera developed technique based on Newton-Raphson algorithm with improved convergence characteristic is derived for the harmonic analysis of system containing nonlinear loads. This technique is able to calculate voltage spectra, power components (P, Q, S) and power factors at each load point under harmonic distortion. Additionally, the proposed technique can handle all types of nonlinear loads assuming that V-I characteristie is known for these loads.