Theoretical estimations for the astrophysical S-factor and the d(alpha,gamma)6Li reaction rates are obtained on the base of the two-body model with the alpha-d potential of a simple Gaussian form, which describes correctly the phase-shifts in the S-,
P-, and D-waves, the binding energy and the asymptotic normalization constant in the final S-state. Wave functions of the bound and continuum states are calculated by using the Numerov algorithm of a high accuracy. A good convergence of the results for the E1- and E2- components of the transition is shown when increasing the upper limit of effective integrals up to 40 fm. The obtained results for the S-factor and reaction rates in the temperature interval 10E+6 K < T < 10E+10 K are in a good agreement with the results of Ref. A.M. Mukhamedzhanov, et.al., Phys. Rev., C 83, 055805 (2011), where the authors used the known asymptotical form of wave function at low energies and a complicated potential at higher energies.
