Making use of two different analytical-numerical methods for capacity computation, we obtain matching to a very high precision numerical values for capacities of a wide family of planar condensers. These two methods are based respectively on the use of the Lauricella function [BeV1, Be1, Be2] and Riemann theta functions [B, Gr, BG1, BG2]. We apply these results to benchmark the performance of numerical algorithms, which are based on adaptive hp-finite element method [HRV1, HRV2, HRV3] and boundary integral method [Ts, JS, AC]. (Andrei Bogatyrëv), sergei.goreinov@ya.ru (Sergei Goreinov), guelpho@mail.ru (Oleg Grigoriev), harri.hakula@aalto.fi (Harri Hakula), vuorinen@utu.fi (Matti Vuorinen)