An improved meshfree method called the Petrov-Galerkin natural element (PG-NE) method is introduced in order to secure the numerical integration accuracy. As in the Bubnov-Galerkin natural element (BG-NE) method, we use Laplace interpolation function for the trial basis function and Delaunay triangles to define a regular integration background mesh. But, unlike the BG-NE method, the test basis function is differently chosen, based on the Petrov-Galerkin concept, such that its support coincides exactly with a regular integration region in background mesh. Illustrative numerical experiments verify that the present method successfully prevents the numerical accuracy deterioration stemming from the numerical integration error.