The Uzawa method is an iterative approach to find approximated solutions to the Stokes equations. This method solves velocity variables involving augmented Lagrangian operator and then updates pressure variable by Richardson update. In this paper, we construct a new version of the Uzawa method to find optimal numerical solutions of the Stokes equations including corner singularities. The proposed method is based on the dual singular function method which was developed for elliptic boundary value problems. We estimate the solvability of the proposed formulation and special orthogonality form for two singular functions. Numerical convergence tests are presented to verify our assertion.