This paper deals with the solution of the spherically symmetric time-dependent Hartree-Fock approximation applied to nuclear giant monopole resonances in the small amplitude regime. The problem is spatially unbounded as the resonance state is in the continuum. The practical requirement to perform the calculation in a finite-sized spatial region yields an artificial boundary, which is not present physically. The question of how to ensure the boundary does not interfere with the internal solution, while keeping the overall calculation time low is studied. Here we propose an absorbing boundary condition scheme to handle the conflict. The derivation, via a Laplace transform method, and implementation is described. An inverse Laplace transform required by the absorbing boundaries is calculated using a method of nonlinear least squares. The accuracy and efficiency of the scheme is tested and results presented to support the case that they are an effective way of handling the artificial boundary.