We study a variational multiparticle correlation theory for the ground state of inhomogeneous Bose fluids at zerotemperature. The ground state wave function of a many-boson system is written as a generalized Jastrow–Feenberg function including multiparticle correlation functions. Euler–Lagrange equations for the one, two, three and four particle correlation functions are derived using a systematic functional theoretical approach within the framework of the generalized hypernetted chain theory (HNC2 equations). The Euler–Lagrange equations for the optimal one and two particle correlation functions are analyzed including the infinite series of the elementary diagrams summed by integral equations. For the determination of the optimal triplet correlation factor we propose a novel iterative procedure in which the HNC2 equation is used to calculate the triplet distribution function.