The mesh free collocation method has a long history of development and is simple to implement and computationally efficient. However, it is often found unstable and less accurate, especially for problems governed by partial differential equations with Neumann (derivative) boundary conditions, such as solid mechanics problems with stress (natural) boundary conditions.