An efficient and stable point collocation scheme based on a meshfree method is studied for the stationary incompressible Navier-Stokes equations. We describe the diffuse derivatives associated with the moving least square method. Using these diffuse derivatives, we propose a point collocation method to fit in solving the Navier-Stokes equations which improves the stability of the direct point collocation scheme. The convergence of the numerical solution is investigated from numerical examples. The driven cavity ow and the backward facing step ow are implemented for the reliability of the scheme. Also, the viscous ow on complicated geometry is successfully calculated such as the ow past a circular cylinder in duct.