By minimizing the global variance in the 1-reduced local-energy matrix El(xl;xl'), subject to the normalization of the 1-reduced density matrix pI(xI; xl'), one derives an integral matrix equation for El(xl; xl') as a functional of p,(xI;xI') at the location (xl;x,') of an arbitrary member of an N (5 2)-particle system.The implications for the possible local improvement in the accuracy of approximate wave functions through the imposition of global constraints are briefly discussed.