The article proposes a numerically approximate method for solving a boundary value problem for an integro-differential equation with a parameter and considers its convergence, stability, and accuracy. The integro-differential equation with a parameter is approximated by a loaded differential equation with a parameter. A new general solution to the loaded differential equation with a parameter is introduced and its properties are described. The solvability of the boundary value problem for the loaded differential equation with a parameter is reduced to the solvability of a system of linear algebraic equations with respect to arbitrary vectors of the introduced general solution. The coefficients and the right-hand sides of the system are compiled through solutions of the Cauchy problems for ordinary differential equations. Algorithms are proposed for solving the boundary value problem for the loaded differential equation with a parameter. The relationship between the qualitative properties of the initial and approximate problems is established, and estimates of the differences between their solutions are given.