This paper introduces a modified edition of classical Cespedes' load flow method to radial distribution system analysis. In the developed approach, a distribution network is modeled in different complex reference systems and reduced to a set of connected equivalent subnetworks, each without resistance, while graph topology and node voltage solution are preserved. Active power losses are then not dissipated in the modeled subnetworks and active power flows can be obtained as a consequence of radiality. Thus, the proposed method preprocesses a series of variable transformations concomitant to an iterative algorithm using a forward-backward sweep to arrive at the load flow solution. The proposed approach has been tested using literature and actual distribution networks, and efficiency improvements are verified in comparison to Cespedes' load flow method. [1,2] provides the steady-state condition for power systems and is one of the most important numerical tools to system planning and designing. In order to design solutions to power distribution system real-time operation [3], modeling and analysis might take into account unbalanced operation and detailed features of each component connected to the networks. On the other hand, in long-term planning stages, some hypotheses and simplifications can be assumed (e.g. balanced operation, constant power load) due to, for instance, uncertainties regarding future load profiles and distributed generation productions. This is the case in adequacy evaluations [4,5], where hundreds of thousands of load flow analysis may be executed, while performance indices are estimated by modeling failure/repair cycles of components and load/generation profiles as stochastic processes.Detailed distribution system modeling and analysis is a well established topic, covered in books such as [6]. Simplified/single-phase load flow modeling and analysis dates mostly to the 80's and 90's and can be divided into two groups: the first group comprises Newtonian based methods adjusted to distribution system analysis [7,8] and the second group includes methods based on iterative forward-backward sweep processes [9,10,11,12,13,14]. The sweep based techniques are known by their efficiency and take advantage of the fact that distribution networks usually are radially operated.Forward-backward sweep methods employ the following general procedures: (a) assuming a flat start or an approximate solution, current or power downstream each node is estimated in a backward sweep (from end nodes towards the substation node); (b) using estimates from the previous step, node voltages are updated in a forward sweep (from the substation node towards end nodes); (c) these two steps are repeated up to the convergence of node voltages.Among the several variations of sweep techniques, the one proposed by Baran and Wu [11,12] stands out by employing a set of equations, known as Distflow branch equations, that recursively relates active power, reactive power and voltage magnitudes. Cespedes [9] steps ahead by addressin...