Troesch’s problem arises in the investigation of the confinement of plasma column by radiation pressure and is defined by
subject to the boundary conditions
The main difficulty associated with Troesch’s problem is the boundary layer near . Accordingly, many researchers try to solve this equation numerically, some of these methods are: Sinc-Collocation method  , the modified homotopy perturbation technique  , a smart nonstandard finite difference for second order nonlinear boundary value problem  , the finite element method and discontinuous Galerkin methods   , and the shifted Jacobi-Gauss collocation method  . A cubic spline collocation method is given in  . All of these methods provide a numerical solution of Troesch’s problem for moderate values of λ. Recently, an accurate asymptotic approximation of Troesch’s problem with large values up to is reported in  . For more details see   .
This paper is devoted to derive a logarithmic finite difference method as a new method to solve the Troesch’s problem for arbitrarily large parameter λ.
The paper remaining of this paper is organized as follows. In Section 2, we derive the numerical scheme as well as a fourth order finite difference scheme. Section 3 is devoted for the numerical results and comparisons with some existing methods. Concluding remarks are given in Section 4.
2. Numerical Method
In order to derive the numerical method for solving Troesch’s, we present the following definition.
Definition: The inverse sine hyperbolic function is defined by
and its derivative is defined by
By using the previous definition, the Troesch’s problem (1) can be written as
which is the main point in this paper.
Using uniform mesh including points , such that , where . We denote the exact and the numerical solutions respectively by and at the grid point . Now by using the second order central finite difference approximation for the second derivative
into the Equation (5), this will lead us to the logarithmic finite difference method
subject to the boundary conditions
The resulting system in (7) represents a nonlinear tridiagonal system in the unknown solution . The details of the logarithmic can be given displayed by the following algorithm.
Algorithm 1. Logarithmic method
1) For using (7) and the boundary condition , we obtain
2) For we use
3) For using (7) and the boundary condition to obtain
Newton’s method is used to solve the nonlinear system which can be described in the following way.
Algorithm 2. Newton’s method
1) For .
2) Solve the tridiagonal linear system
for the unknown vector Z using Crout’s method. The elements of the function and the Jacobian matrix J are defined as follows
and the Jacobian matrix has the tridiagonal structure
with elements defined by
3) Update your solution by using
4) Repeat steps (1)-(3) till the following condition
is satisfied. The constant is assumed to be small.
The initial guess vector is taken as the linear interpolation between the given boundary conditions, and has the following form
The previous method is of second order accuracy and it works for wide range values of λ.
A fourth order finite difference method can be used for solving the Troche’s method for limited values of λ. This method can be given as follows.
By using the forth order approximation of the second derivative
Using this approximation and after some manipulation we will end with the following scheme
The resulting method (14) is of fourth order accuracy. The numerical solution in this case can be also obtained by solving the nonlinear tridiagonal system. Newton’s method is used to solve this system as we did before.
3. Numerical Results
To test the efficiency of the proposed method, we choose and , for the moderate values of . Comparison with some existing methods is given in Tables 1-5. In Figure 1 we display the numerical solution for . The numerical results indicate that the proposed method is highly accurate comparing to the published works.
For large values of λ, we choose . In Table 6 and Table 7, we display the numerical solution for . To the knowledge of the authors no published work discussed these values of λ before. Again we test our method severely for and 106, the method won in this case as well and the results are given in Table 8, these results are given for the first time for such values of λ up to my knowledge. We display the numerical solution for in Figures 2-5.
Finally, we tested the fourth order method and make some comparison with the logarithmic method, we have noticed that, the results produced by this method are quite good for the small values of λ only, see Table 9, and the method is handicapped and fail for larger values of λ.
Table 1. Troesch’s problem with λ = 0.5 and h = 0.0005.
Table 2. Troesch’s problem with λ = 1.0 and h = 0.0005.
Table 3. Troesch’s problem with λ = 5.0 and h = 0.0005.
Table 4. Troesch’s problem with λ = 10 and h = 0.0005.
Table 5. Troesch’s problem with λ = 60 and h = 0.0001.
Figure 1. Numerical solution for λ = 3, 5, 8, 10.
Table 6. Numerical solution using the logarithmic method and h = 0.0001.
Table 7. Numerical solution using the logarithmic method and h = 0.0001.
Table 8. Numerical solution using Logarithmic method for large values of λ and h = 0.0001.
Table 9. Comparison of the proposed methods at λ = 0.9 with h = 0.0005.
Figure 2. Numerical solution with λ = 100.
Figure 3. Numerical solution with λ = 200.
Figure 4. Numerical solution with λ = 1000.
Figure 5. Numerical solution with λ = 10,000.
In this work, we have derived a logarithmic finite difference method to overcome the difficulty in solving the Troesch’s problem for large values of λ. This progress is very important, since all existing methods were trying to obtain the numerical solution for Troesch’s problem for large values of λ. The logarithmic method which we have derived is able and succeeds to get the numerical solution of Troesch’s problem for . To recap things, a new logarithmic finite difference method is derived and can provide the numerical solutions for large values of λ. I think that, to the best of our knowledge, the method and some of the given results are published for the first time.
 Erdogan, U. and Ozis, T. (2011) A Smart Nonstandard Finite Difference Scheme for Second Order Nonlinear Boundary Value Problems. Journal of Computational Physics, 230, 6464-6474.
 Temimi, H. and Kurkcu, H. (2014) An Accurate Asymptotic Approximation and Precise Numerical Solution of Highly Sensitive Troesch’s Problem. Applied Mathematics and Computation, 235, 253-260.
 Scott, M (1975) On the Conversion of Boundary-Value Problems into Stable Initial Value Problems via Several Invariant Imbedding Algorithms. In: Aziz, A.K., Ed., Numerical Solution of Boundary-Value Problems for Ordinary Differential Equations.