Abstract

In this paper, an efficient modification is introduced into the well-known homotopy perturbation method which proves very effective to control the convergence region of approximate solution. By using this scheme, explicit exact solution is calculated in the form of a convergent power series with easily computable components. The proposed algorithm is tested on some neutral stiff functional-differential equations with proportional delays (DDEs). Numerical results explicitly reveal the complete reliability, efficiency and accuracy of the suggested technique. It is observed that the approach may be implemented on other models of physical nature.

Key words: Homotopy perturbation method, systems of stiff delay differential equations, Pade approximation.