Point orthogonal projection onto an algebraic surface is a very important topic in computer-aided geometric design and other fields. However, implementing this method is currently extremely challenging and difficult because it is difficult to achieve to desired degree of robustness. Therefore, we construct an orthogonal polynomial, which is the ninth formula, after the inner product of the eighth formula itself. Additionally, we use the Newton iterative method for the iteration. In order to ensure maximum convergence, two techniques are used before the Newton iteration: (1) Newton’s gradient descent method, which is used to make the initial iteration point fall on the algebraic surface, and (2) computation of the foot-point and moving the iterative point to the close position of the orthogonal projection point of the algebraic surface. Theoretical analysis and experimental results show that the proposed algorithm can accurately, efficiently, and robustly converge to the orthogonal projection point for test points in different spatial positions.