A fourth/second order accurate collocation method for singularly perturbed two-point boundary value problems using tension splines (Q5943702)

The concept of tension spline is derived from the classical Euler-Bernoulli model of an elastic beam which is subjected to a tensile force proportional to the square of the tension parameters. Tension splines have been applied to the problem of convex interpolation [cf. \textit{D. G. Schweikert}, J. Math. Phys. 45, 312-317 (1966; Zbl 0146.14102)]. In order to avoid ill-conditioned collocation matrices for large tension parameters, the present paper combines the collocation method for two-point boundary value problems [cf. \textit{J. E. Flaherty} and \textit{W. Mathon}, SIAM J. Sci. Stat. Comput. 1, 260-289 (1980; Zbl 0465.65045)] with an algorithm based on the B-spline representation of the tension splines. This approach provides a regular collocation matrix, even in the limiting case when the perturbation parameter \(\varepsilon > 0\) tends to zero. In addition, the B-spline representation enables to establish a global convergence result for the collocation method [cf. \textit{M. Marušić} and \textit{M. Rogina}, Adv. Comput. Math. 6, No. 1, 65-76 (1996; Zbl 0876.65058)]. Numerical experiments suggest that the error analysis could be helpful in developing an algorithm for mesh size selection.