  We use the linear programming algorithm introduced by Akulin et al. [V. M.
Akulin, G. A. Kabatiansky, and A. Mandilara, Phys. Rev. A 92, 042322 (2015)] to
perform best separable approximation on two-qutrit random density matrices. We
combine the numerical results with theoretical methods in order to generate
random representative families of positive partial transposed bound entangled
(BE) states and analyze their properties. Our results are disclosing that for
the two-qutrit system the BE states have negligible volume and that these form
tiny `islands' sporadically distributed over the surface of the polytope of
separable states. %We devise a method for estimating numerically the average
thickness of these formations and their frequency of occurrence. The detected
families of BE states are found to be located under a layer of pseudo one-copy
undistillable negative partial transposed states with the latter covering the
vast majority of the surface of the separable polytope.
