Radial distribution systems (RDS) require special load flow methods to solve power flow equations owing to their high R/X ratio. Increasing use of power electronic devices and effect of magnetic saturation cause harmonics in RDS. This paper proposes a novel algorithm to compute the power flow solution of a RDS accounting for all the harmonic components. It uses a recursive solution technique. The proposed method uses a novel dynamic data structure reported in the paper. The proposed method is tested upon a 33-bus RDS and the results are reported.