In recent years, mesh free methods have attracted more and more attention and regarded as a potentialnumerical method in computational mechanics. However, the integrands emerged in these methods are rational functions and ordinary gauss quadrature does not lead to suitable results as in FEM. But using a polynomialapproximation for integrands over small local quadrature domain, a novel point integration technique is presented for the LRPIM mesh free method. Radial basis functions (RBFs) which possess the Delta function property are used toconstruct shape functions. The accuracy, stability and efficiency of the present method are proved in some numerical examples