We focus on the problem of evaluating the velocity field outside a solid object moving in an incompressible Stokes flow using the boundary integral formulation. For points near the boundary, the integral is nearly singular, and accurate computation of the velocity is not routine. One way to overcome this problem is to regularize the integral kernel. The method of regularized Stokeslet (MRS) is a systematic way to regularize the kernel in this situation. For a specific blob function which is widely used, the error of the MRS is only of first order with respect to the blob parameter. We prove that this is the case for radial blob functions with decay property $ϕ(r)$=$\mathcal{O}$($r^{−3−α}$) when $r→∞$ for some constant $α>1$. We then find a class of blob functions for which the leading local error term can be removed to get second and third order errors with respect to blob parameter. Since the addition of these terms might give a flow field that is not divergence free, we introduce a modification of these terms to make the divergence of the corrected flow field close to zero while keeping the desired accuracy. Furthermore, these dominant terms are explicitly expressed in terms of blob function and so the computation time is negligible.