A finite differences MATLAB code for the numerical solution of second order singular perturbation problems

We show the main features of the MATLAB code HOFiD_UP for solving second order singular perturbation problems. The code is based on high order finite differences, in particular on the generalized upwind method. Within its simplicity, it uses order variation and continuation for solving any difficult nonlinear scalar problem. Several numerical tests on linear and nonlinear problems are considered. The best performances are reported on problems with perturbation parameters near the machine precision, where most of the codes for two-point BVPs fail.

[1]  U. Ascher,et al.  A New Basis Implementation for a Mixed Order Boundary Value ODE Solver , 1987 .

[2]  Pierluigi Amodio,et al.  A Stepsize Variation Strategy for the Solution of Regular Sturm‐Liouville Problems , 2011 .

[3]  Pierluigi Amodio,et al.  High order generalized upwind schemes and numerical solution of singular perturbation problems , 2007 .

[4]  Pierluigi Amodio,et al.  A Deferred Correction Approach to the Solution of Singularly Perturbed BVPs by High Order Upwind Methods: Implementation Details , 2009 .

[5]  Giuseppina Settanni,et al.  Variable Step/Order Generalized Upwind Methods for the Numerical Solution of Second Order Singular Perturbation Problems 1 2 , 2009 .

[6]  L. Brugnano,et al.  Solving differential problems by multistep initial and boundary value methods , 1998 .

[7]  Donato Trigiante,et al.  Convergence and stability of boundary value methods for ordinary differential equations , 1996 .

[8]  Francesca Mazzia,et al.  Numerical approximation of nonlinear BVPs by means of BVMs , 2002 .

[9]  Pierluigi Amodio,et al.  A Matrix Method for the Solution of Sturm-Liouville Problems 1 , 2011 .

[10]  Donato Trigiante,et al.  A Hybrid Mesh Selection Strategy Based on Conditioning for Boundary Value ODE Problems , 2004, Numerical Algorithms.

[11]  Jeff Cash,et al.  Lobatto deferred correction for stiff two-point boundary value problems , 1998 .

[12]  Margaret H. Wright,et al.  A Deferred Correction Method for Nonlinear Two-Point Boundary Value Problems: Implementation and Numerical Evaluation , 1991, SIAM J. Sci. Comput..

[13]  Francesca Mazzia,et al.  Hybrid Mesh Selection Algorithms Based on Conditioning for Two-Point Boundary Value Problems , 2006 .

[14]  Jeff Cash,et al.  An automatic continuation strategy for the solution of singularly perturbed linear two-point boundary value problems , 1995 .

[15]  Gerald Moore,et al.  An automatic continuation strategy for the solution of singularly perturbed nonlinear boundary value problems , 2001, TOMS.

[16]  Francesca Mazzia,et al.  A new mesh selection algorithm, based on conditioning, for two-point boundary value codes , 2005 .

[17]  M. Stynes,et al.  Numerical methods for singularly perturbed differential equations : convection-diffusion and flow problems , 1996 .

[18]  Robert D. Russell,et al.  Algorithm 569: COLSYS: Collocation Software for Boundary-Value ODEs [D2] , 1981, TOMS.

[19]  Lawrence F. Shampine,et al.  A User-Friendly Fortran BVP Solver , 2006 .

[20]  Pierluigi Amodio,et al.  High-order finite difference schemes for the solution of second-order BVPs , 2005 .