ACM Transactions on Mathematical Software (TOMS)
Algorithm 743: WAPR--a Fortran routine for calculating real values of the W-function
ACM Transactions on Mathematical Software (TOMS)
Spectral methods in MatLab
Solution of the transcendental equation wew = x
Communications of the ACM
Journal of Computational and Applied Mathematics
New perturbation-iteration solutions for Bratu-type equations
Computers & Mathematics with Applications
The Legendre wavelet method for solving initial value problems of Bratu-type
Computers & Mathematics with Applications
Sinc-Galerkin method for numerical solution of the Bratu's problems
Numerical Algorithms
Letter to the editor: On the integral solution of the one-dimensional Bratu problem
Journal of Computational and Applied Mathematics
A simple solution of the Bratu problem
Computers & Mathematics with Applications
Hi-index | 0.48 |
We introduce two methods for simultaneously approximating both branches of a two-branched function using Chebyshev polynomials. Both schemes remove the pernicious, convergence-wrecking effects of the square root singularity at the limit point where the two branches meet. The "Chebyshev-Shafer" method gives the approximants as the solution to a quadratic equation; the "mapping" algorithm makes a quadratic change of variable. For both, the only input information is the set of values of f(x) at a set of discrete points. There is little to choose between the two schemes in accuracy, but the single expansion/mapping method is more flexible in that it can accommodate different ranges on the two branches. The eigenrelation for the one-dimensional Bratu equation is an interesting example because the upper branch is also singular at another point besides the limit point; this, too, can be removed by subtracting the asymptotic solution, which is the Lambert W-function, from the upper branch only. When the domain of the variable is infinite, the quadratic change of variable can still be applied by substituting rational Chebyshev functions, which are basis functions for an unbounded interval, for Chebyshev polynomials. We illustrate this by approximating the real-valued root of the Brill quintic, u5 - u - λ, λ ∈ [-∞, ∞], which was first solved by Hermite using elliptic modular functions more than a century ago.