Corner-force obstruction for a cubic polytrope (source code)

= Corner-force obstruction for a cubic polytrope
{title2=$\nabla\Phi_{\rm int}(0)=0\ne\nabla\Phi_N(0)$}

The sine-product <cubic polytropic interior> has zero <mass density> gradient at every corner. Its hydrostatic potential $C-2K\rho$ therefore has zero gradient there. In contrast, at a corner of a cube containing a positive <mass density>, every interior mass element contributes a gravitational-potential gradient pointing into the same coordinate octant. Each component of the gradient is strictly nonzero. This directly proves failure of global isolated gravitational matching, even though the interior hydrostatic and Poisson equations hold.