Solving two-point boundary value problems using combined homotopy perturbation method and Green’s function method

In this paper, an algorithm is proposed for the solution of second-order boundary value problems with two-point boundary conditions. The Green’s function method is applied first to transform the ordinary differential equation into an equivalent integral one, which has already satisfied the boundary conditions. And then, the homotopy perturbation method is used to the resulting equation to construct the numerical solution for such problems. Numerical examples demonstrate the efficiency and reliability of the algorithm developed, it is quite accurate and readily implemented for both linear and nonlinear differential equations with homogeneous and nonhomogeneous boundary conditions. Furthermore, the lower order approximation is of higher accuracy for most cases. Some other extended applications of this algorithm are also exhibited.