1. Introduction
Consider Emden–Fowler-type SBVPs with derivative dependence, as expressed by the following equation:
where
and
are real constants. Here,
is the source function dependent on both
and
The conditions
and the allowance of
to be discontinuous at
lead to a reduction of the problem to double SBVPs [
1]. Such problems are prevalent in various areas of astrophysics, including thermal explosion modeling in a rectangular slab [
2,
3], heat source measurements in human heads [
4], oxygen concentration within spherical cells [
5], shallow membrane cap theory [
6], heat conduction problems [
7], unsteady Poiseuille flow in a pipe [
8], electroelastic dynamic problems [
9], and heat explosions [
10].
Solving Emden–Fowler SBVPs with derivative dependence is crucial for predicting system behaviors, such as changes in pressure, density, or temperature within stars or gaseous spheres. The solutions provide valuable insights into the structure and evolution of these systems. Finding numerical solutions for derivative-dependent second-order singular differential equations is particularly challenging due to strong nonlinearity from derivatives in the source function and the singular behavior at the origin. The motivation lies in developing numerical methods that require less computational effort while maintaining high accuracy. The collocation method has gained popularity with the widespread availability and efficiency of computers, being applied to problems in physics, engineering, and other fields.
The existence and uniqueness of the estimated solution of Equation (
1) were readily obtained in [
11,
12,
13,
14,
15,
16] under the conditions of
with
and
for
on (0, 1].
Several numerical methods have been developed to solve Equation (
1) when
, including the cubic spline method [
17], the finite difference method [
18,
19,
20,
21,
22], the Adomian decomposition method (ADM) [
23,
24,
25,
26,
27,
28], the B-spline collocation method [
29,
30], the classical polynomial approximation method [
31], etc. However, there are limited techniques available for solving Emden–Fowler SBVPs with derivative dependence. In 2014, Singh et al. [
26] discussed the Adomian decomposition technique to solve an original utilizing Green’s function. In 2018, Roul [
32] presented an improved normal homotopy analysis method to solve derivative-dependent SBVPs, and in 2019, Roul et al. [
33] discussed quintic spline interpolation. In 2020, Shahni et al. [
34] established an approximate solution for Emden–Fowler-type SBVPs with derivative dependence using Bernstein polynomials. Upon examining existing techniques, limitations were identified, such as a significant amount of computational work, especially for nonlinear singular boundary value problems. Therefore, there is a need for more efficient numerical methods that can overcome these limitations and provide a more accurate solution for nonlinear singular boundary value problems.
This work introduces a constructive approach for solving Emden–Fowler-type SBVPs with derivative dependence. In
Section 2, the differential equation is converted into its equivalent Fredholm integral form. In
Section 3, a collocation technique based on Chebyshev polynomials (CCM) is employed to obtain the system of nonlinear equations upon transformation of the Fredholm integral equation. Subsequently, Newton’s method is implemented to solve the system and obtain the required solution. In
Section 4, the algorithm for the methodology is provided for implementing the method. In
Section 5, error analysis is included to assess the accuracy of the current method. In
Section 6, the maximum absolute error of the current method is computed for various examples using
and
norm. These numerical results are compared with those obtained using the existing BCM method [
34]. The residual errors between CCM and a previously established method, i.e., BCM, are also compared graphically.
6. Numerical Illustration
We utilized MATLAB (R2015a) to determine maximum absolute errors using both
and
norms for various examples to assess the accuracy of the current approach. Subsequently, we compared these results with those obtained using the BCM in the different tables and also graphically compared residual errors. We define
and
norm errors as follows:
and
where
and
represent the approximate and exact solutions, respectively.
Furthermore, the residual error is defined as:
where
and
are obtained by using the CCM and the BCM, respectively.
Example 1. The equivalent integral form iswhere The exact solution of differential Equation (43) is We present a comparison of errors using the
norm of the CCM with the BCM for
and
in
Table 1 and
Table 2, respectively, and a comparison of errors using the
norm of the CCM with the BCM for
and
in
Table 3 and
Table 4, respectively, for Example 1. We note that achieving the desired level of accuracy is more effectively accomplished with the CCM compared to the BCM. Furthermore,
Figure 1 illustrates a comparison of residual errors between the current technique and the BCM, revealing that the residual errors of the CCM are significantly lower than those of the BCM.
Example 2. Its equivalent integral form iswhere The exact solution of differential Equation (44) is We present a comparison of errors of the CCM with those of the BCM for
using the
and
norms in
Table 5 and
Table 6, respectively for Example 2. We note that achieving the desired level of accuracy is more effectively accomplished with the CCM compared to the BCM. Furthermore,
Figure 2 illustrates a comparison of residual errors between the current technique and the BCM, revealing that the residual errors of the CCM are significantly lower than those of the BCM.
Example 3. Its integral form iswhere The exact solution of differential Equation (45) is We present a comparison of errors using the
and
norms of the CCM with those of the BCM for
for Example 3 in
Table 7 and
Table 8, respectively. We clearly observe that achieving the desired level of accuracy is more effectively accomplished with the CCM compared to the BCM. Furthermore,
Figure 3 illustrates a comparison of residual errors between the current technique and the BCM, revealing that the residual errors of the CCM are significantly lower than those of the BCM.
Example 4. The equivalent integral form iswhere The exact solution of differential Equation (46) is We present a comparison of errors using the
and
norms of the CCM with those of the BCM for
and
for Example 4 in
Table 9. We clearly observe that achieving the desired level of accuracy is more effectively accomplished with the CCM compared to the BCM. Furthermore,
Figure 4 illustrates a comparison of residual errors between the current technique and the BCM, revealing that the residual errors of the CCM are significantly lower than those of the BCM.
Example 5. Consider a numerical problem without exact solution Its equivalent form iswhereand , , , . We compare the absolute difference of estimated solutions (
) of the CCM with the BCM in
Table 10. It can be seen from the table that fewer errors in numerical solutions are obtained by the present method than the BCM [
39].