Numerical methods were employed to find periodic orbits and the construction of invariant manifolds of a saddle periodic orbit. A plane endomorphism which was a model of coupled oscillators was studied. The appearance and destruction of invariant curves were determined by numerical simulation methods. Parameter values associated with the rising of periodic orbits were obtained.