We discuss the classical problem of how to pick N weighted points on a d−dimensional manifold so as to obtain a reasonable quadrature ruleThis problem, naturally, has a long history; the purpose of our paper is to propose selecting points and weights so as to minimize the energy functional d(x, y) is the geodesic distance and d is the dimension of the manifold. This yields point sets that are theoretically guaranteed, via spectral theoretic properties of the Laplacian −∆, to have good properties. One nice aspect is that the energy functional is universal and independent of the underlying manifold; we show several numerical examples.