In this report an efficient numerical algorithm for solving linear two-point boundary value problem for loaded differential equations with delay is proposed. The Dzhumabaev parameterization method is applied to construct this efficient numerical algorithm. Based on this method, necessary and sufficient conditions for solvability of considering problem have been established. The application of the Dzhumabaev parameterization method reduces the linear two-point boundary value problem for loaded differential equations with delay to the solving a system of algebraic equations and Cauchy problems for ordinary differential equations.