In this paper, we develop a numerical approximation to the solution of a system of integral equations arising from a nonlocal hyperbolic reliability model. Convergence of this numerical method is proved. And this numerical scheme is used to study the behavior of the solution. [Copyright &y& Elsevier]