  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.
