A local symmetric weak form (LSWF) for linear potential problems is developed, and a truly meshless method, based on the LSWF and the moving least squares approximation, is presented for solving potential problems with high accuracy. The essential boundary conditions in the present formulation are imposed by a penalty method. The present method does not need a "finite element mesh", either for purposes of interpolation of the solution variables, or for the integration of the "energy". All integrals can be easily evaluated over regularly shaped domains (in general, spheres in three-dimensional problems) and their boundaries. No post-smoothing technique is required for computing the derivatives of the unknown variable, since the original solution, using the moving least squares approximation, is already smooth enough. Several numerical examples are presented in the paper. In the example problems dealing with Laplace & Poisson's equations, high rates of convergence with mesh refinement for the Sobolev norms ∥ · ∥0 and ∥ · ∥1 have been found, and the values of the unknown variable and its derivatives are quite accurate. In essence, the present meshless method based on the LSWF is found to be a simple, efficient, and attractive method with a great potential in engineering applications.