Numerical method for singularly perturbed fourth order ordinary differential equations of convection-diffusion type

Abstract

In this paper, we have proposed a numerical method for singularly perturbed fourth order ordinary differential equations of convection-diffusion type. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and finite difference method. In order to get a numerical solution for the derivative of the solution, the given interval is divided into two subintervals called inner region (boundary layer region) and outer region. The shooting method is applied to inner region whereas for the outer region, standard finite difference method is applied. Necessary error estimates are derived. Computational efficiency and accuracy are verified through numerical examples.