  We present a practical method to solve Faddeev three-body equations at
energies above three-body breakup threshold as integral equations in coordinate
space. This is an extension of previously used method for bound states and
scattering states below three-body breakup threshold energy. We show that
breakup components in three-body reactions produce long-range effects on
Faddeev integral kernels in coordinate space, and propose numerical procedures
to treat these effects. Using these techniques, we solve Faddeev equations for
neutron-deuteron scattering to compare with benchmark solutions.
