This document presents a novel algorithm to compute harmonic power flow in radial distribution systems accounting for unbalanced conditions and harmonic pollution. It proposes using a dynamic data structure (DDS) to model each branch and store relevant parameters. A recursive solution technique is used to calculate voltages starting from the farthest branches and moving upstream. The method is tested on IEEE 13, 37, and 123 bus test systems. Results show the total harmonic distortion at each node to evaluate the proposed approach. The algorithm aims to efficiently solve for harmonic power flow in unbalanced and non-sinusoidal radial distribution systems.