Interest in meshless methods has grown rapidly in recent years in solving boundary value problems (BVPs) arising in science and engineering. In this paper, we present the moving least square radial reproducing polynomial (MLSRRP) meshless method as a generalization of the moving least square reproducing kernel particle method (MLSRKPM). The proposed method is established upon the extension of the MLSRKPM basis by using the radial basis functions. Some important properties of the shape functions are discussed. An interpolation error estimate is given to assess the convergence rate of the approximation. Also, for some class of time-dependent partial differential equations, the error estimate is acquired. The efficiency of the present method is examined by several test problems. The studied method is applied to the parabolic two-dimensional transient heat conduction equation and the hyperbolic two-dimensional sine-Gordon equation which are discretized by the aid of the meshless local Petrov–Galerkin (MLPG) method.