In this paper of numerical nature, we test the Lebesgue constant of several pointsets on the disk Ω and propose new ones that enjoy low Lebesgue constant. Furthermore we extend some results in [15], analysing the case of Bos arrays whose radii are nonnegative Gegenbauer-Gauss-Lobatto nodes with exponent α, noticing that the optimal α still allow to achieve pointsets on Ω with low Lebesgue constant Λn for degrees n ≤ 30. Next we introduce an algorithm that through optimization determines pointsets with the best known Lebesgue constants for n ≤ 25. Finally, we determine theoretically a pointset with the best Lebesgue constant for the case n = 1.