Numerical method of solution to loaded nonlocal boundary value problems for ordinary differential equations

A numerical method is suggested for solving systems of nonautonomous loaded linear ordinary differential equations with nonseparated multipoint and integral conditions. The method is based on the convolution of integral conditions into local ones. As a result, the original problem is reduced to an initial value (Cauchy) problem for systems of ordinary differential equations and linear algebraic equations. The approach proposed is used in combination with the linearization method to solve systems of loaded nonlinear ordinary differential equations with nonlocal conditions. An example of a loaded parabolic equation with nonlocal initial and boundary conditions is used to show that the approach can be applied to partial differential equations. Numerous numerical experiments on test problems were performed with the use of the numerical formulas and schemes proposed.