Аннотация:We obtain a quadrature formula for the harmonic double layer potential with continuous density defined on a smooth surface (closed or open). Numerical tests confirm that this formula has higher accuracy than the standard quadrature formula. The advantage of the new formula is especially noticeable near the surface, where the values calculated by the standard formula tend to infinity, while the new formula provides acceptable accuracy.