In this paper, a meshless method is introduced for NDT computation. Compared with the conventional FEM, it can avoid the onerous mesh generation and updating, only a distribution of points and the description of the boundaries are needed. The mathematical background for moving least square approximation employed in the method is given, and the numerical implementation is discussed. Application of the method for MFL computation and comparison with the results from FEM are also presented.