In order to improve the approximation of spatial derivatives without meshes,a set of mesh free methods for spatial derivatives is developed,which is compatible with the coordinates of Cartesian,cylindrical and spherical.Based on the comparisons between numerical and theoretical solutions,errors and convergences are assessed by a posteriori methods,which show the approximations for functions and derivatives are of the second accuracy order,and the scale of the support domain has some influences on numerical errors but not on accuracy orders.With a discrete scale h=0.01,the relative errors of the numerical simulation for the selected functions and their derivatives are within 0.65%.