In this study, a bi-Gaussian approximation to the number density function (NDF) is used to approximate the both the maximum-entriopy NDF and closing moment fluxes of fourth-order 14-moment maximum-entropy closure associated with fully multi-dimensional kinetic theory.