  We extend the piecewise orthogonal collocation method to computing periodic
solutions of coupled renewal and delay differential equations. Through a
rigorous error analysis, we prove convergence of the relevant finite-element
method and provide a theoretical estimate of the error. We conclude with some
numerical experiments to further support the theoretical results.
