摘要:We calculate the double charmonium production cross-section within the framework of Bethe-Salpeter Equation in the electron-positron annihilation, at center of mass energy GeV, that proceeds through the exchange of a single virtual photon. In this calculation, we make use of the full Dirac structure of 4D BS wave functions of these charmonia, with the incorporation of all the Dirac covariants (both leading and subleading). The calculated cross-sections for the double charmonium productions for final states, (), (), (), and (), are close to experimental data and in broad agreement with results of other theoretical models.