classical gravity, but now with Lambda

what about a Lagrange density \(\mathcal{L}(\Phi,\nabla\Phi)\) for Newtonian gravity including the cosmological constant, which leads to a solution \(\Phi = -GM/r + \lambda r^2\) for a spherically symmetric matter distribution?