This paper presents a new simulation method to analyze the transient stability of a power system including three phase unbalanced impedances. The method is based on the phase coordinate method, because it is easy to analyze the power system which has elements of unbalanced three phase impedances by the phase coordinate method. This paper consists of two main parts: unbalanced power flow calculation to determine the initial conditions, and the transient stability simulation of power system following faults.