  We study the processes, $e^- e^+ \rightarrow \gamma^* \gamma^* \rightarrow
J/\psi +J/\psi$, and $e^- e^+ \rightarrow \gamma^* \gamma^* \rightarrow \eta_c+
\eta_c$ at $\sqrt{s}=10.6$ GeV in the framework of $4\times 4$ Bethe-Salpeter
equation. For $J/\psi+J/\psi$ production, the dominant contribution is through
fragmentation process, while for $\eta_c+\eta_c$ production, the quark
rearrangement diagrams contribute. Our results of cross section for
$J/\psi+J/\psi$ and $\psi(2S)+\psi(2S)$ are compatible with the experimental
upper limits set by Belle Collaboration, while in the absence of experimental
data for $\eta_c(1S)+\eta_c(1S)$, and $\eta_c(2S)+\eta_c(2S)$ production, we
have given theoretical prediction of their cross sections, and compared with
NRQCD prediction.
