This paper addressed the issue of the three-phase load flow program for electric power distribution systems that include distributed generations. The equipment models were selected in order to consider imbalances among phases for the load flow analysis of distribution systems. Also, power equations and measurement functions are newly derived. The load flow analysis program developed in this paper was tested for the propriety of algorithm and convergence characteristics by case studies on test systems in various scales and types.