In this paper, we study the numerical method for solving singularly perturbed partial delay differential equations with integral boundary conditions. Due to a small perturbation parameter acted on the higher derivative, solution of the problem exhibit a boundary layers at left and right end plane of the domain. An interior layer is also formed because of the presence of large delay on the space variable. A hybrid numerical method is proposed in the spatial direction, and an implicit Euler method is used in temporal direction. The proposed hybrid scheme constitute of cubic spline method in the boundary layer region and a classical finite difference method in the outer layer region. The integral boundary condition is treated using Simpson’s 1/3 rule. The uniform stability and convergence analysis for the proposed scheme is studied. The developed method is uniformly convergent with second-order in space and first-order in time. Two numerical test examples are considered to validate a theoretical result. The numerical results are in agreement with the theoretical estimations.
MSC Classification: 65M06 , 65M12 , 65M22 , 65M25