The cross section for the scattering of very slow neutrons by deuterons is calculated by numerical methods. Polarization is completely neglected and the wave equation for the process is set up in such a form as to take correctly into account exchange effects between the incident neutron and the neutron initially in the deuteron. This wave equation is then replaced by an integral equation the solution of which is correctly symmetrized and has the right asymptotic value to describe the scattering process. The numerical integration is performed by replacing the integral equation by a finite set of simultaneous linear algebraic equations. The work is greatly simplified by the use of a sum of two Gauss functions to approximate the ground state deuteron wave function. It is assumed throughout this paper that the interactions between like and unlike particles are equal and are of the general form ${V}_{\mathrm{ij}}=\ensuremath{-}[(1\ensuremath{-}g\ensuremath{-}{g}_{1}\ensuremath{-}{g}_{2}){P}_{\mathrm{ij}}+g{P}_{\mathrm{ij}}{Q}_{\mathrm{ij}}+{g}_{1}+{g}_{2}{Q}_{\mathrm{ij}}]J({r}_{\mathrm{ij}}),$ where the symbols have their usual meanings and where $J({r}_{\mathrm{ij}})$ is a Gauss function. The calculation is carried out for two sets of $g$'s. For the first set, ${g}_{1}={g}_{2}=0$, $g=0.2$, the cross section is found to be equal to 4.57\ifmmode\times\else\texttimes\fi{}${10}^{\ensuremath{-}24}$ ${\mathrm{cm}}^{2}$, and for the second set of $g$'s, ${g}_{2}=2$, $g=0.22\ensuremath{-}{g}_{2}$, ${g}_{1}=0.25\ensuremath{-}0.8{g}_{2}$, the value of the cross section is found to be equal to 6.91\ifmmode\times\else\texttimes\fi{}${10}^{\ensuremath{-}24}$ ${\mathrm{cm}}^{2}$. The experimental value is at least 20 percent smaller than the first of these values.