Abstract: A numerical approach is presented for simulating the stability of harmonic resonance overvoltage in power systems and the relationship is discussed between the stability of harmonic solutions and the change of parameters of system by using the methods of Galerkin and Runge Kutta.It can be used for describing further bifurcation phenomenon in power systems.