Three phase power flow analysis of an unbalanced power distribution system decouples voltage changes resulting from nodal admittance matrix into one contribution from a real part, conductance matrix, and other contribution from an imaginary part, and a susceptance matrix. A first voltage change and second voltage change resulting from conductance and susceptance matrices are determined respectively. The voltages of a node of the power distribution system are determined as a combination of first and second voltages.