Chebfun Solutions to a Class of 1D Singular and Nonlinear Boundary Value Problems


The Chebyshev collocation method implemented in Chebfun is used in order to solve a class of second order one-dimensional singular and genuinely nonlinear boundary value problems. Efforts to solve these problems with conventional ChC have generally failed, and the outcomes obtained by finite differences or finite elements are seldom satisfactory. We try to fix this situation using the new Chebfun programming environment. However, for tough problems, we have to loosen the default Chebfun tolerance in Newton’s solver as the ChC runs into trouble with ill-conditioning of the spectral differentiation matrices. Although in such cases the convergence is not quadratic, the Newton updates decrease monotonically. This fact, along with the decreasing behaviour of Chebyshev coefficients of solutions, suggests that the outcomes are trustworthy, i.e., the collocation method has exponential (geometric) rate of convergence or at least an algebraic rate. We consider first a set of problems that have exact solutions or prime integrals and then another set of benchmark problems that do not possess these properties. Actually, for each test problem carried out we have determined how the Chebfun solution converges, its length, the accuracy of the Newton method and especially how well the numerical results overlap with the analytical ones (existence and uniqueness).


Calin-Ioan Gheorghiu
Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy


Chebfun; differential equation; non-linearity; singularity; convergence; Bernstein growth;
improper integrals; boundary layer

Paper coordinates

C.-I. Gheorghiu, Chebfun solutions to a class of 1D singular and nonlinear boundary value problems, Computation, 10 (2022) no. 7, art. id. 116,


About this paper



Publisher Name


Print ISSN
Online ISSN


google scholar link

1. Gheorghiu, C.I. Pseudospectral solutions to some singular nonlinear BVPs. Numer. Algor. 2015, 68, 1–14. [CrossRef]
2. Gheorghiu, C.I. A Third-Order Nonlinear BVP on the Half-Line. Available online: (accessed on 30 January 2020).
3. Agarwal, R.P.; O’Regan, D. A singular Homann differential equation. ZAMM Z. Angew. Math. Mech. 2003, 83, 344–350. [CrossRef]
4. Agarwal, R.P.; O’Regan, D. Singular initial and boundary value problems with sign changing nonlinearity. IMA J. Appl. Math 2000, 65, 173–198. [CrossRef]
5. Trefethen, L.N.; Birkisson, A.; Driscoll, T.A. Exploring ODEs; SIAM: Philadelphia, PA, USA, 2018; pp. 51–62.
6. Driscoll, T.A.; Bornemann, F.; Trefethen, L.N. The chebop system for automatic solution of differential equations. Bit Numer. Math. 2008, 48, 701–723. [CrossRef]
7. Russell, R.D.; Shampine, L.F. Numerical methods for singular boundary value problems. SIAM J. Numer. Anal. 1975, 12, 13–36. [CrossRef]
8. Singh, R.; Guleria, V.; Singh, M. Haar wavelet quasi-linearization method for numerical solution of Emden–Fowler type equations. Math. Comput. Simul. 2020, 174, 123–133. [CrossRef]
9. Boyd, J.P. Chebyshev and Fourier Spectral Methods, 2nd ed.; Dover: Mineola, NY, USA, 2001; pp. 19–57.
10. Kierzenka, J.; Shampine, L.F. A BVP solver based on residual control and the Maltab PSE. ACM Trans. Math. Softw. 2001, 27, 299–316. [CrossRef]
11. Granas, A.; Guenther, R.B.; Lee, J.W. Continuation and shooting methods for boundary value problems of Bernstein type. J. Fixed Point Theory Appl. 2009, 6, 27–61. [CrossRef]
12. Van Gorder, R.A. Exact fist integral for a Lane-Emden equation of the second kind modeling a thermal explosion in a rectangular slab. New Astron. 2011, 16, 492–497. [CrossRef]
13. Agarwal, R.P.; O’Regan, D. Singular Problems Modelling Phenomena in the Theory of Pseudoplastic Fluids. ANZIAM J. 2003, 45, 167–179. [CrossRef]
14. Flagg, R.C.; Luning, C.D.; Perry, W.L. Implementation of New Iterative Techniques for Solutions of Thomas-Fermi and EmdenFowler Equations. J. Comput. Phys. 1980, 38, 396–405. [CrossRef]
15. Gheorghiu, C.-I. Accurate Spectral Collocation Computation of High Order Eigenvalues for Singular Schrödinger Equations. Computation 2021, 9, 2. [CrossRef]
16. Duggan, R.C.; Goodman, A.M. Pointwise Bounds for a Nonlinear Heat Conduction Model of the Human Head. Bull. Math. Biol. 1986, 48, 229–236. [CrossRef]
17. Callegari, A.J.; Reiss, E.L.; Keller, H.B. Membrane Buckling: A Study of Solution Multiplicity. Commun. Pure Appl. Math. 1971, 24, 499–527. [CrossRef]
18. Agarwal, R.P.; O’Regan, D. Nonlinear Super-linear Singular and Non-singular Second Order Boundary Value Problems. J. Differ. Equ. 1998, 143, 60–95. [CrossRef]
19. Mazzia, F.; Sgura, I. Numerical approximation of nonlinear BVPs by means of BVMs. Appl. Numer. Math. 2002, 42, 337–352. [CrossRef]

Related Posts