Bs linear multistep methods on non uniform mesh

BS methods are a special class of Linear Multistep Methods dened using B-spline functions. These methods are always convergent and have good stability properties when used as Boundary Value Methods. In addition, if k is the number of steps, a C^k spline of degree k + 1 can be computed with low computational cost and this serves as a continuous extension to the solution. It is shown that the continuous solution and the discrete solution both share the same order of convergence. In this paper we introduce this class of methods in the general case of a non-uniform mesh and we present numerical results showing their performance when dealing with some singularly perturbed Boundary Value Ordinary Differential Equations.