This paper presents an iteration method of solving non linear boundary integral equations (BIE) of the plane Navier Stokes problem,which each step of the iteration is to solve a linear BIE of nonhomogeneous Stokes problem. Using the quadrature methods of [1],we give a new algorithm with a high order accuracy.The algorithm not only saves work,but also the accuracy can be improved by Richardson extrapolation.