This paper describes the comparison of Newton’s method and Steepest Descent method for determining the coordinates of earthquake source. Here, the Steepest Descent method is used because it is a Newton-based method. The earthquake case used in this research is Italian earthquake August 24, 2016, which has a seismic phase of Pg. The calculation was supported by Azimuth Coordinate equations to find the coordinates and Haversine formula to find the distance between five earthquake stations to the earthquake source. The final result of calculations was path’s graph from the iteration of Steepest Descent method. Moreover, the results will be compared with the results of Newton’s method that has been successfully approaching the point of earthquake source in the same case study of previous research. The result shows the number of final iterations of two methods using tolerance number 0.01, minimum velocity number 3093 m/s and three cases of initial guess in the form of city coordinate of Rome, Milan, and Palermo. Newton’s method generates 12 iterations in every case, Steepest Descent method generate 7, 6, 5 iteration respectively. However, the final numerical errors for Rome, Milan and Palermo initial guess are 0.1598 by Newton’s method, while Steepest Descent method are 0.1566, 0.1567 and 0.1567.