When the matrices $A$ and $Q$ have special structure, the structure-preserving algorithm was used to compute the stabilizing solution of the complex matrix equation $X+A^TX^{-1}A=Q.$ In this paper, we study the numerical methods to solve the complex symmetric stabilizing solution of the general matrix equation $X+A^TX^{-1}A=Q.$ We not only establish the global convergence for the methods under an assumption, but also show the feasibility and effectiveness of them by numerical experiments.