Earthquake-induced slope instability represents a significant risk, frequently resulting in considerable damage to transportation infrastructure and cities. This research enhances Sarma's method for determining the seismic slope safety factor through limit equilibrium theory, utilizing a Python program to improve computational efficiency. The modified method was implemented in the Sichuan earthquake zone in China, achieving significant insights for the assessment and mitigation of seismic slope hazards. The findings indicate that the seismic slope safety factor exhibits a sinusoidal pattern over time, with variations determined by slope height, angle, wave velocity, and frequency. A progressive decrease in the safety factor was noted with increasing slope height, whereas smaller slope angles diminished the extent of these variations. The interaction between seismic wave frequency and slope height is critical, as specific frequencies resonate with the slope's natural frequency, thereby diminishing stability. The improved Sarma method, along with the computational efficiency of Python, presents a viable approach for analyzing seismic slope stability. The precise evaluation of seismic slope stability is crucial for identifying high-risk zones and executing prompt protective strategies to reduce disaster risks in regions susceptible to earthquakes.