Numerical solution of ordinary differential equations of the second order with variable coefficients by the finite differences method using difference schemes of the fourth order of accuracy

Often when we solve the problems of mathematical physics, we obtain differential equations, the analytical solution of which is impossible. Then it is necessary to solve the problem using numerical methods. We can develop a universal method for solving such equations using difference schemes. We will use difference schemes of the fourth order of accuracy and solve equations with variable coefficients.