Matlab code to fill a brick-shaped (or other shape) 3D domain with random spheres in contact with each other, with spherical radii following a prescribed probability distribution. For details, see the preprint arXiv:2107.01233 or the corresponding paper in SoftwareX