In this paper, we show that the set of quaternionic polynomials is dense in the Bergman spaces of slice regular functions in the unit ball, both of the first and of the second kind. Several proofs are presented, including constructive methods based on the Taylor expansion and on the convolution polynomials. In the last case, quantitative estimates in terms of higher‐order moduli of smoothness and of best approximation quantity are obtained.