  It is shown that difficulties in reconciling the values of the cross section
e+e- -> J/psi+eta_c measured at BELLE and calculated within non-relativistic
QCD (NRQCD) are caused not by some misinterpretation of the data or other
exotic explanations, but by poor applicability of NRQCD to such processes.
  We use general theory of hard exclusive processes in QCD together with more
realistic models of charmonium wave functions, and show that the BELLE result
can be naturally explained.
