We calculate the double charmonium production cross-section within the framework of 4×4 Bethe-Salpeter Equation in the electron-positron annihilation, at center of mass energy s=10.6GeV, 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, (J/Ψ,ηc), (Ψ′,ηc), (J/Ψ,ηc′), and (Ψ′,ηc′), are close to experimental data and in broad agreement with results of other theoretical models.