Resonance effect on the orbital elements of geosynchronous artificial satellite due to the non-zonal geopotential has been calculated. For the perturbation of a artificial satellite, perturbation effects due to the non-zonal geopotential is less than due to the $J_2$ or Luni-Solar perturbation, but non-zonal harmonics resonance exist. So, we calculate the perturbation of geosynchronous artificial satellite orbit due to the non-zonal harmonics resonance. The effect on the orbit eccentricity of non-zonal harmonics resonance is represented by a phase plane plot of ec. The effect on the orbit eccentricity of non-zonal harmonics resonance is represented by a phase plane plot of $e_c$ verse $e_s$. The evolution of mean longitude and semi-major axis are obtained.