We compute angular phase-space integrals with three and four denominators analytically, working within dimensional regularisation via the Mellin-Barnes (MB) representation. The approach converts multifold MB integrals into real parametric integrals and expresses all results in terms of Goncharov polylogarithms (GPLs). For three denominators, all-massless results are obtained to O(ε²) and the single-massive case to O(ε); for four denominators, both the massless and single-massive cases are solved to O(ε⁰). Integrals with multiple massive momenta follow from a partial fraction decomposition reducing them to the single-massive case. Recursion relations relating integrals with higher denominator powers to master integrals are derived. These are essential ingredients to solving full phase-space integrals.