Based on the basic idea of iterative refinement and the special structure of coefficient matrices,a method called the refined iterative method for solving special block tridiagonal Toeplitz linear systems of equations is proposed,which could greatly improve the accuracy of exact solution.The method is characterized by its simple scheme,numerical stability,high precision and fast convergence.Finally,we apply this method to surface-fitting of cubic uniform B-spline sur...