In this paper, a system of two circular concentric non-touching tori is considered. An electrostatic charge can be applied to each of the tori. The problem is to find the density of charge distribution over the surfaces of the tori, taking into account the Coulomb interaction between the surfaces. The required density is found numerically by the method of successive approximations based on the fact that in the static case, the total tangential strength of the Coulomb forces is zero at each of the points on the surface of each of the tori. The corresponding functional is constructed, the problem of numerical minimization of which is solved by the gradient descent method.