Charmonium Production in e^+e^- Annihilation at sqrt{s}=10.6GeV

We review the charmonium production in $e^+e^-$ annihilation at $\sqrt{s}=10.6GeV$. NRQCD predictions for the $J/\psi$ production are compared with the recent measurements by BaBar and Belle. A NRQCD calculation for the $D$-wave charmonium production is reported to further test the color-octet mechanism and detect the $J^{PC}=2^{--}$ state. The issue of double $c \bar c$ production in $e^+e^-$ annihilation observed by Belle is also discussed.