Abstract
This paper presents a novel numerical approach to addressing three types of high-order singular boundary value problems. We introduce and consider three modified Chebyshev polynomials (CPs) of the third kind as proposed basis functions for these problems. We develop new derivative operational matrices for the three modified CPs of the third kind by deriving formulas for their first derivatives. Our approach follows a unified method for numerically handling singular differential equations (DEs). To transform these equations into algebraic systems suitable for numerical treatment, we employ the collocation method in combination with the introduced operational matrices of derivatives of the modified CPs of the third kind. We address the convergence examination for the three expansions in a unified manner. We present numerous numerical examples to demonstrate the accuracy and efficiency of our unified numerical approach.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
In a single dimension, a boundary value problem involves an ordinary differential equation (ODE) and conditions that specify the solution and/or its derivatives at a minimum of two points. The quantity of governing conditions matches the order of the ODE. Boundary value problems (BVPs) may be found in various scientific and technological contexts, including but not limited to boundary layer theory in fluid mechanics, heat power transfer theory, and optimization theory. The important book of Agarwal [1] was devoted to studying theoretically the existence and uniqueness theorems for BVPs. In addition, some applications in which these types of equations appear were also presented. Due to the importance of BVPs, they were numerically treated by different techniques. Some algorithms were proposed to treat the linear and non-linear BVPs. For example, numerical algorithms based on Galerkin methods were proposed in [2, 3] to treat some linear BVPs. In contrast, some other algorithms were handled to treat the non-linear BVPs (see, for example, [4, 5]). Among the other schemes that were proposed to treat various BVPs are the differential transform method in [6], operational methods [7], Sinc-collocation method in [8], Green’s function method in [9], iterative methods [10], and spline methods [11]. Orthogonal polynomials (OPs), a topic linked to numerous crucial areas of analysis, have witnessed significant development in recent years. The fields of DEs and integral equations rely heavily on OPs. Additionally, they provide rather general and illustrative examples of specific circumstances in the theory of orthogonal systems. Among the most utilized polynomials in different applications are the Jacobi polynomials (see, for example, [12, 13]). Four kinds of CPs are special polynomials of Jacobi polynomials. It is well-known that these polynomials have their crucial parts in approximation theory. For some contributions that approximate solutions are expressed in terms of different CPs, see, for example, [14, 15].
The use of spectral approaches dramatically influences the extent of numerical analysis; see [16,17,18] for an illustration of their importance in engineering and fluid dynamics. The spectral method’s fundamental principle is to treat approximations as the sum of a few chosen special functions or polynomials. There are three methods to get numerical solutions to differential and integral equations, and their variants are widely used. The three variations differ in the choice of trial and test functions. Using the underlying conditions of the problem as a guide, the Galerkin technique is conditionally applied by choosing the trial functions that meet those conditions; see, for example, [19,20,21]. A critical difference between the tau and Galerkin approaches is that the tau approach is not restricted in choosing basis functions. The collocation method is widely used because it may be used for every differential or integral equation; see for example [22,23,24].
Due to the great importance of the different OPs in various applications, particularly in numerical analysis, their investigations are vital. Recently, Abd-Elhameed in [25] studied a type of OPs, generalized Jacobi polynomials. This paper derived connections between these kinds of OPs and other polynomials. It is worth noting that some studies were performed on some modified CPs to be capable of treating some types of DEs. In [26], the authors introduced certain modified CPs of the second kind and used them to treat the third-order Emden-Fowler type DEs, while the authors in [27] introduced modified CPs of the third-kind (MCPs3K) and employed them to provide solutions to some singular high-order DEs.
Operational matrices of integrals and derivatives are useful in numerical analysis. These matrices are essential to finding numerical solutions to different DEs. (see, for example, [28, 29]). Obtaining an operational matrix of derivatives of certain basis functions relies on expressing the first derivative of this basis in terms of their originals.
Currently, the article introduces three kinds of MCPs3K that will be utilized to solve three types of singular BVPs using the typical collocation method. A unified approach will be proposed for solving the three singular BVPs. Consequently, this paper contributes to:
-
Introducing three new types of modified CPs.
-
Establishing the expressions of the first-order derivative of the modified polynomials.
-
Constructing three new operational matrices of derivatives of the modified polynomials.
-
Designing a unified approach based on applying the typical collocation algorithm to treat three types of singular high-order BVPs.
-
Investigating the convergence of the three expansions of the modified polynomials.
-
Examining the accuracy of the presented algorithm via presenting some illustrative examples.
To the best of our knowledge, the established operational matrices are novel. This gives us motivation to utilize them to solve some singular differential equations.
The paper will proceed as described below. An account of third-kind CPs is provided in Sect. 2, along with introducing three modified versions of third-kind CPs. Three Galekin operational metrics of the three modified CPs are presented in Sect. 3. To deal with the three distinct kinds of singular high-order BVPs, the standard collocation technique is applied in conjunction with the three Galerkin operational matrices proposed in Sect. 4. The convergence and error analysis of the three modified Chebyshev expansions are examined in Sect. 5 via a unified approach. In Sect. 6, some numerical results are displayed. In the final section, we offer some final thoughts.
2 An Account on Third-Kind of CPs and Three Modified Ones
In this section, we will introduce some characteristics of the MCPs3K. In addition, three kinds of modified CPs of the third-kind will be considered. They will be used as the basis functions in our proposed numerical algorithm.
2.1 An Account on Third-Kind CPs
It is well-known CPs of the third-kind are particular polynomials of the Jacobi polynomials that can be expressed as
Below is a recurrence relation that can be utilized to construct these polynomials:
The shifted CPs on [a, b] can be defined as
\({\mathcal {V}}^{*}_{m}(t),\ m\ge 0\) are orthogonal on [a, b] whose orthogonality relation is
where \({w_{0}}(t)=\sqrt{\displaystyle \frac{t-a}{b-t}}\).
2.2 Three Kinds of MCPs3K
This section introduces three kinds of OPs that are considered modifications of CPs. The main reason for introducing these polynomials is that these polynomials will meet the underlying conditions of the singular BVPs under investigation.
For our present purposes, and for every positive integer n with \(n\ge 2\), we introduce the three sets of OPs \(\left\{ \phi _{n,j}(t)\right\} _{j\ge 0}\), \(\left\{ \psi _{n,j}(t)\right\} _{j\ge 0}\) and \(\left\{ \chi _{n,j}(t)\right\} _{j\ge 0}\) defined as
Remark 1
The three selected basis functions that are defined in (5), (6), and (7) will meet the boundary conditions of the BVPs that will be considered and treated in the subsequent section.
Remark 2
To give a unified approach for our proposed numerical algorithm, we will denote \(\Upsilon _{n,j}(t)\) to be one of the three polynomials \(\phi _{n,j}(t)\), \(\psi _{n,j}(t)\) and \(\chi _{n,j}(t)\).
Remark 3
According to Remark 2 and the orthogonality relation (4), it can be seen \(\Upsilon _{n,j}(t)\) are orthogonal on [a, b] regarding \({w_{n}}(t)=\displaystyle \frac{{w_{0}}(t)}{\Omega ^{2}_{n}(t)}\), with \({w_{0}}(t)=\sqrt{\displaystyle \frac{t-a}{b-t}}\), and
Furthermore, we have the following unified orthogonality relation:
3 Introducing Three New Operational Matrices of Derivatives
Here, we create three new operational matrices of derivatives for the three MCPs3K introduced in (5), (6), and (7). The formulas for these polynomials’ first derivatives provide a path to these operational matrices. Three theorems will be In this concern, we will state and prove three theorems.
Theorem 1
Let \(\phi _{n,i}(t)\) be the modified CPs that are defined in (5). For any integer i with \(i\ge 1\), the first-derivative \(D \phi _{n,i}(t)\) can be expressed as
where
and
with
Proof
Let us first suppose \([a,b]=[-1,1]\). In such case, we have
A process of induction is followed on i. In the case of \(i=1\), it is easy to prove that relation (10) holds. Suppose (10) applies to i and \((i-1)\). The next step is to prove that
With reference to (2), it is obvious that
Differentiating (15) with respect to t yields
If the induction hypothesis is applied to \(D\phi _{n,i}(t)\) and \(D\phi _{n,i-1}(t)\) in (17), then we get
The substitution of (15) written in the form
into relation (18), and using \({\bar{\vartheta }}_{i,i}(n)=0\) for all i, will yield-after some algebraic computations- the following formula
Formulae (26) and (27) imply the following relations:
It is not difficult to demonstrate the above relations leads to the following formula:
and accordingly, it can be seen that Eq. (18) turns to the form (14).Now, we can finish the proof by substituting \(\dfrac{2t-a-b}{b-a}\) for t in (14). \(\square\)
The following two theorems give the expressions for the first-order derivative of the two modified Cps that are defined in (6) and (7).
Theorem 2
Let \(\psi _{n,i}(t)\) be the modified CPs that are defined in (5). For any integer i with \(i\ge 1\), \(D\psi _{n,i}(t)\) can be expressed as
where
and
with
Proof
Similarly to the proof of Theorem 1. \(\square\)
Theorem 3
Let \(D \chi _{n,i}(t)\) be the modified CPs that are defined in (5). For any integer i with \(i\ge 1\), \(D \chi _{n,i}(t)\) can be expressed as
where
and
with
Proof
Similar to the proof of Theorem 1. \(\square\)
Remark 4
Aiming to utilize a unified numerical algorithm for three different BVPs using the three modified CPs in (5), (6), and (7), we will unify the results of Theorems 1,2 and 3 in the following theorem.
Theorem 4
Let \(\Upsilon _{n,i}(t)\) be one of the three modified polynomials that are defined in (5), (6), and (7). \(D \Upsilon _{n,i}(t),\,i\ge 1,\) can be expressed as
where
and
Proof
Theorem 4 is a direct consequence of Theorems 1, 2 and 3. \(\square\)
3.1 Function Approximation
Now, if consider a function \(u(t)\in \, L^{2}_{w_{n}}[a,b]\) that can be expanded as:
where
Let us assume that we can approximate the series in Eq. (31) as
where
3.2 Three Operational Matrices of Modified CPs
Based on the three derivative formulas of the basis functions considered in (5),(6), and (7) that are unified in Theorem 4, the following corollary presents the operational matrix of the derivative of \({\varvec{\Upsilon _{n,N}}}(t)\).
Corollary 1
In terms of \({\varvec{\Upsilon _{n,N}}}(t),\ \displaystyle \frac{d{\varvec{\Upsilon _{n,N}}}(t)}{dt}\), the following formula holds:
where \({\varvec{\delta _{n}}}(t)=\left( \delta _{n,0}(t),\delta _{n,1}(t),\dots , \delta _{n,N}(t)\right) ^T and \ G_{n}=\big (g^{(n)}_{ij}\big )_{0\le i,j\le N}\). The entries of \(G_{n}\) have the following form:
Now, we give the operational matrices of derivatives for the three considered modified CPs as examples. For \(N=6\), the corresponding operational matrices for every \(n\ge 2,\) have the following forms:
-
(i)
If \(\Upsilon _{n,i}(t)=\phi _{n,i}(t)\)
$$\begin{aligned} \hspace{-40pt}G_{n}=\displaystyle \frac{4}{b-a} \begin{pmatrix} 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ n+1 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 5-2 n &{} n+2 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 3 n+2 &{} 5-2 n &{} n+3 &{} 0 &{} 0 &{} 0 &{} 0 \\ 10-4 n &{} 3 (n+1) &{} 5-2 n &{} n+4 &{} 0 &{} 0 &{} 0 \\ 5 n+3 &{} 10-4 n &{} 3 n+4 &{} 5-2 n &{} n+5 &{} 0 &{} 0 \\ 15-6 n &{} 5 n+4 &{} 10-4 n &{} 3 n+5 &{} 5-2 n &{} n+6 &{} 0 \\ \end{pmatrix} _{7\times 7}, \end{aligned}$$(36) -
(ii)
If \(\Upsilon _{n,i}(t)=\psi _{n,i}(t)\)
$$\begin{aligned} G_{n}=\displaystyle \frac{4}{b-a} \begin{pmatrix} 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 2 n+1 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 1 &{} 2 (n+1) &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 6 n+2 &{} 1 &{} 2 n+3 &{} 0 &{} 0 &{} 0 &{} 0 \\ 2 &{} 6 n+3 &{} 1 &{} 2 (n+2) &{} 0 &{} 0 &{} 0 \\ 10 n+3 &{} 2 &{} 6 n+4 &{} 1 &{} 2 n+5 &{} 0 &{} 0 \\ 3 &{} 10 n+4 &{} 2 &{} 6 n+5 &{} 1 &{} 2 (n+3) &{} 0 \\ \end{pmatrix} _{7\times 7}, \end{aligned}$$(37) -
(iii)
If \(\Upsilon _{n,i}(t)=\chi _{n,i}(t)\)
$$\begin{aligned} G_{n}=\displaystyle \frac{4}{b-a} \begin{pmatrix} 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 2 (n+1) &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ -1 &{} 2 n+3 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 6 n+5 &{} -1 &{} 2 (n+2) &{} 0 &{} 0 &{} 0 &{} 0 \\ -2 &{} 6 (n+1) &{} -1 &{} 2 n+5 &{} 0 &{} 0 &{} 0 \\ 10 n+8 &{} -2 &{} 6 n+7 &{} -1 &{} 2 (n+3) &{} 0 &{} 0 \\ -3 &{} 10 n+9 &{} -2 &{} 6 n+8 &{} -1 &{} 2 n+7 &{} 0 \\ \end{pmatrix}_{7\times 7}. \end{aligned}$$(38)
Due to Corollary 1, \(\displaystyle \frac{d^{r}\varvec{\Upsilon _{n,N}}(t)}{dt^{r}}\) can be easily evaluated. This result can be seen in the corollary below.
Corollary 2
For every positive integer r, \(\displaystyle \frac{d^{r}\varvec{\Upsilon _{n,N}}(t)}{dt^{r}}\) can be represented as
with \(\varvec{\eta _{n}}^{(r)}(t)=\sum \limits _{k=0}^{r-1}G^{k}_{n}\,\varvec{\delta _{n}}^{(r-k-1)}(t)\).
Proof
Repeated application of (35) yields (39). \(\square\)
Remark 5
In the upcoming section, we will treat three types of singular BVPs with domain [0, 1]. So, we will employ the three operational matrices established in Sect. 3 for \(a=0,b=1\).
4 Galerkin Operational Matrix Method for Treating Three Types of Singular High-order BVPs
In this section, we show how to use the Galerkin operational matrices of derivatives given in (39) to numerically solve the high-order singular BVPs in the form [30,31,32]
subject to one of the following three types of the BCs:
where L is a nonlinear differential operator of an order less than n.
The main idea in the proposed algorithm is to convert (40) subject to the non-homogeneous BCs conditions (41), (42) or (43) into a modified one subject to the homogeneous BCs. For this purpose, we set
The determination of the polynomials \(P_{n-1}(t)\) the type of the boundary condition governed by the BVP (40). For the BCs (41), \(P_{n-1}(t)\) takes the form
For the BCs (42), and (43), to determine \(P_{n-1}(t)\), it can be assumed that
and the coefficients \(A_{i}\, (i=\frac{n}{2},\dots ,n-1)\) and \(B_{i}\, (i=\frac{n+1}{2},\dots ,n-1)\) can be computed by solving the following two systems of equations:
and
Now, making use of the transformation (44), it can be seen Eq. (40) governed by the on-homogeneous BCs (41), (42) or (43) turns into a modified one governed by the homogeneous BCs. More precisely, the modified equation is
subject to one of the following the homogeneous BCs:
with
Now, we will propose a unified collocation approach namely shifted Chebyshev third kind Galerkin collocation operational matrix methods, (SC3GCOMMB1), (SC3GCOMMB2) and (SC3GCOMMB3) using the three considered basis functions to solve the singular BVPs (49)
for the purpose of finding approximate solutions to Eq. (49) directed with the three homogeneous BCs (50), (51) or (52), respectively. Through the use of Eqs. (35) and (39), we are able to determine the high-order derivatives of the approximate solution \(\tilde{{u}}_{N}(t)\). Hence, the residual of Eq. (49) can be written as
where \(q\in \left\{ n,\,\frac{n}{2},\,\frac{n-1}{2}\right\}\), according to the choice of basis (53). The application of shifted Chebyshev third kind collocation operational matrix method to solve numerically (49) subject to one of BCs (50), (51) or (52). The collocation points are selected as the \((N+1)\) zeros of \({\mathcal {V}}^{*}_{N+1}(t)\). Consider them \(t_{i},\, i=0,1,\dots ,N\). Thus
Any appropriate algorithm may be used to find the solution \(\tilde{{u}}_{N}(t)\) to the \((N+1)\) nonlinear algebraic equations generated by the system (55).
5 Convergence and Error Analysis of the Three Modified Chebyshev Expansions
Here, we offer the results of our research into the error analysis of the proposed algorithm. Suppose that
where \(q\in \left\{ n,\tfrac{n}{2},\,\tfrac{n-1}{2}\right\}\) according to the considered BCs (50)–(52), respectively, with \(g(t)\in C^{n}[0,1]\le M,\) for all \(x\in [0,1]\), and M is a positive integer.
Theorem 5
Assume that \({\tilde{u}}(t)\in L^n_{w_{q}}[0,1]\) be the exact solution of the differential Eq. (49) subject to one of the BCs (50), (51) or (52). let \({\tilde{u}}(t)\) be expanded as
The series converges uniformly to \({\tilde{u}}(t)\), and the coefficients \(c_{i}\) satisfy the the following inequality
Proof
The application of the orthogonality relation (4) yields
The substitution by: \(2\, t-1=\cos \theta\) in (57) yields
The integration by parts n times leads to
where \(R_{i}(\theta )\) is a trigonometric polynomial in \(\cos \theta\), \(\sin \theta\). So, using \(|\cos \theta |\le 1, |\sin \theta |\le 1\) and after some manipulations we get the desired result. \(\square\)
Theorem 6
If we consider an approximate solution of (49) in the form \({\tilde{u}}_{N}(t)=\displaystyle \sum _{i=0}^{N}c_{i}\, \Upsilon _{n,i}(t)\). Under the same conditions stated in Theorem 5, then the global error satisfies the following inequality:
Proof
As a consequence of the orthogonality property of \(\Upsilon _{n,i}(t),\, i\ge 0\), we have
and accordingly, Theorem 5 leads to
and therefore, the following inequality also holds (see, [33])
In conclusion, the last inequality suggests that
Theorem 6 is now proved. \(\square\)
6 Numerical Results
This section aims to demonstrate the usefulness of our suggested collocation method by providing some instances. Additionally, we provide some comparisons with other approaches. Furthermore, in every instance within this section, the mistakes are assessed using the maximum norm, namely \(E_{N}=\max \limits _{0 \le t \le 1}|u(t)-u_{N}(t)|\).
Remark 6
We point out that the numerical examples presented in this section to solve Eq. (40) subject to the given boundary conditions will be obtained using Mathematica 13.3 on a computer system equipped with an Intel(R) Core(TM) i9-10850 CPU operating at 3.60 GHz, featuring 10 cores and 20 logical processors.
Example 1
Consider the following two singular BVPs:
subject to the BCs:
or
The exact solution of (60) is: \(u(t)=1+t^{6}.\) Our unified collocation algorithm is applied, and the approximate solution is given by
where
and \(c_{i}=0, i=2,3,\dots ,N\). In view of (44), we obtain \(u_{N}(t)=1+t^6\), which is the exact solution.
Example 2
Consider the singular BVP: [30, 31, 34, 35]
subject to the BCs:
The exact solution is \(u(t)=t^{3}e^{t}\).
We apply our suggested method SC3GCOMMB1 for various values of N. Table 1 shows the maximum errors that occur accompanied by the CPU time (in seconds), and Table 2 compares the results of various built techniques. Furthermore, Fig. 1 presents logerror for various N to show the stability of solutions of Example 2. The application of our proposed method SC3GCOMMB1.
Example 3
Consider the following singular BVPs:
subject to the BCs:
where g(t) is chosen such that the exact solution of (66) is
and the exact solution of (67) is \(u(t)=\dfrac{1}{\sqrt{4+t^2}}\).
The application of proposed method SC3GCOMMB2 to solve (66) gives
where \(\Upsilon _{i}(t)= t^{2}(1-t)^2\, V^{*}_{j}(t)\). In such case, we have
The greatest errors obtained by solving (67) using our suggested method SC3GCOMMB2 with various values of N are shown in Table 3, accompanied by the CPU (in seconds), and a comparison of the various techniques is shown in Table 4. Furthermore, Fig. 2a shows the computational error norm with \(N = 13\), and Fig. 2b shows the approximate and precise solutions.
Example 4
Consider the fourth-order Emden-Flower type equation [32, 36, 37]
subject to the initial-boundary conditions
with the exact solution: \(u(t)=t^{2}e^{t}.\)
Table 6 compares our technique to the methods proposed in [32, 36, 37], while Table 5 shows the highest errors that arise from using our method (SC3GCOMMB2). Further, Fig. 3a shows the computational error norm with \(N = 11\), and Fig. 3b shows the approximate and precise solutions.
Example 5
Consider the non-linear Emden-Flower type equation
subject to the initial-boundary conditions
where f(t) is selected to meet the exact solution \(u(t)=t^4(1-t)^3\,e^{t}.\)
The greatest errors that resulted from using our SC3GCOMMB3 technique are shown in Table 7. Figure 4a shows the computational error norm with \(N = 14\), and Fig. 4b shows the approximate and precise solutions.
Example 6
Consider the non-linear Lane-Emden equation
subject to the initial-boundary conditions
where the explicit exact solution is not available, so the following error norm is used to check the accuracy in this case:
Utilizing SC3GCOMMB2 and using \(N=1,3,6,9,12,14,\) give the numerical results shown in Table 8.
7 Results and Discussions
New spectral solutions of the singular high-order type equations were implemented and analyzed. Three types of high-order BVPs were solved by applying a unified spectral collocation algorithm. Three modified CPs were utilized as basis functions. The core of the unified algorithm is the establishment of three operational matrices of the modified CPs. The convergence analysis of the three modified Chebyshev expansions was investigated. Several singular BVPs of different orders and various BCs were solved via the proposed technique. The numerical results show the agreement of the resulting approximate solution with the exact ones. As an expected future work, we aim to introduce other modified OPs to treat other types of DEs. The presented method can be developed to solve partial differential equations, such as diffusion-wave and heat equations. The computed operational matrices can also be generalized for solving fractional differential equations.
Data Availability
Not Applicable.
References
Agarwal, R.P.: Boundary Value problems From Higher Order Differential Equations. World Scientific, (1986)
Abd-Elhameed, W.M., Alkenedri, A.M.: Spectral solutions of linear and nonlinear BVPs using certain Jacobi polynomials generalizing third-and fourth-kinds of Chebyshev polynomials. CMES Comput. Model. Eng. Sci. 126(3), 955–989 (2021)
Doha, E.H., Abd-Elhameed, W.M., Bhrawy, A.H.: New spectral-Galerkin algorithms for direct solution of high even-order differential equations using symmetric generalized Jacobi polynomials. Collect. Math. 64(3), 373–394 (2013)
Abdelhakem, M., Fawzy, M., El-Kady, M., Moussa, H.: An efficient technique for approximated BVPs via the second derivative Legendre polynomials pseudo-Galerkin method: Certain types of applications. Results Phys. 43, 106067 (2022)
Verma, A.K., Pandit, B., Verma, L., Agarwal, R.P.: A review on a class of second order nonlinear singular BVPs. Mathematics 8(7), 1045 (2020)
Rezapour, S., Tellab, B., Deressa, C.T., Etemad, S., Nonlaopon, K.: Hu-type stability and numerical solutions for a nonlinear model of the coupled systems of Navier BVPs via the generalized differential transform method. Fractal Fract. 5(4), 166 (2021)
Abdelhakem, M., Alaa-Eldeen, T., Baleanu, D., Alshehri, M.G., El-Kady, M.: Approximating real-life BVPs via Chebyshev polynomials’ first derivative pseudo–Galerkin method. Fractal Fract. 5(4), 165 (2021)
Qiu, Wenlin, Da, Xu., Zhou, Jun, Guo, Jing: An efficient sinc-collocation method via the de transformation for eighth-order boundary value problems. J. Comput. Appl. Math. 408, 114136 (2022)
Ali, F., Ali, J., Uddin, I.: A novel approach for the solution of BVPs via Green’ function and fixed point iterative method. J. Appl. Math. Comput. 66, 167–181 (2021)
Nikooeinejad, Z., Heydari, M., Loghmani, G.B.: A numerical iterative method for solving two-point BVPs in infinite-horizon nonzero-sum differential games: Economic applications. Math. Comput. Simul. 200, 404–427 (2022)
Khalid, A., Ghaffar, A., Naeem, M.N., Nisar, K.S., Baleanu, D.: Solutions of BVPs arising in hydrodynamic and magnetohydro-dynamic stability theory using polynomial and non-polynomial splines. Alex. Eng. J. 60(1), 941–953 (2021)
Doha, E.H., Abd-Elhameed, W.M., Ahmed, H.M.: The coefficients of differentiated expansions of double and triple Jacobi polynomials. Bull. Iranian Math. Soc. 38(3), 739–765 (2012)
Alsuyuti, M.M., Doha, E.H., Ezz-Eldien, S.S.: Galerkin operational approach for multi-dimensions fractional differential equations. Commun. Nonlinear Sci. Numer. Simul. 114, 106608 (2022)
Kumar, S., Nieto, J.J., Ahmad, B.: Chebyshev spectral method for solving fuzzy fractional Fredholm–Volterra integro-differential equation. Math. Comput. Simul. 192, 501–513 (2022)
Ji, T., Hou, J., Yang, C.: The operational matrix of Chebyshev polynomials for solving pantograph-type Volterra integro-differential equations. Adv. Continuous Discrete Models 2022(1), 1–16 (2022)
Boyd, J.P.: Chebyshev and Fourier Spectral Methods. Courier Corporation, (2001)
Trefethen, L.N.: Spectral Methods in MATLAB, volume 10. SIAM, (2000)
Hesthaven, J., Gottlieb, S., Gottlieb, D.: Spectral Methods for Time-Dependent Problems, volume 21. Cambridge University Press, (2007)
Abd-Elhameed, W.M., Alsuyuti, M.M.: Numerical treatment of multi-term fractional differential equations via new kind of generalized Chebyshev polynomials. Fractal Fract. 7(1), 74 (2023)
Hafez, R.M., Zaky, M.A.: High-order continuous Galerkin methods for multi-dimensional advection-reaction-diffusion problems. Eng. Comput. 36(4), 1813–1829 (2020)
Abdelhamied, D., Abdelhakem, M., El-Kady, M., Youssri, Y.H.: Modified shifted Chebyshev residual spectral scheme for even-order BVPs. Math. Sci. Lett. 12(1), 14–18 (2023)
Moghadam, A.A., Soheili, A.R., Bagherzadeh, A.S.: Numerical solution of fourth-order BVps by using Lidstone-collocation method. Appl. Math. Comput. 425, 127055 (2022)
Mohebbi, A.: Crank- Nicolson and Legendre spectral collocation methods for a partial integro-differential equation with a singular kernel. J. Comput. Appl. Math. 349, 197–206 (2019)
Khader, M.M., Eid, A., Adel, M.: Implementing the Vieta–Lucas collocation optimization method for MHD Casson and Williamson model under the effects of heat generation and viscous dissipation. J. Math., 2022, (2022)
Abd-Elhameed, W.M.: Novel formulae of certain generalized Jacobi polynomials. Mathematics 10(22), 4237 (2022)
Abd-Elhameed, W.M., Ahmed, H.M.: Tau and Galerkin operational matrices of derivatives for treating singular and Emden-Fowler third-order-type equations. Int. J. Mod. Phys. C 33(5), 2250061–17 (2022)
Abd-Elhameed, W.M., Al-Harbi, M.S., Amin, A.K., Ahmed, H.M.: Spectral treatment of high-order Emden-Fowler equations based on modified Chebyshev polynomials. Axioms 12(2), 99 (2023)
Saadatmandi, A., Dehghan, M.: A new operational matrix for solving fractional-order differential equations. Comput. Math. Appl. 59(3), 1326–1336 (2010)
Doha, E.H., Abd-Elhameed, W.M., Bassuony, M.A.: On using third and fourth kinds Chebyshev operational matrices for solving Lane-Emden type equations. Rom. J. Phys. 60(3–4), 281–292 (2015)
Aruna, K., Ravi Kanth, A.S.V.: A novel approach for a class of higher order nonlinear singular boundary value problems. Inter. J. Pure Apll. Math. 84(4), 321–329 (2013)
Hasan, Y.Q., Zhu, L.M.: Solving singular boundary value problems of higher-order ordinary differential equations by modified Adomian decomposition method. Commun. Nonlinear Sci. Numer. Simul. 14(6), 2592–2596 (2009)
Iqbal, M.K., Abbas, M., Zafar, B.: New quartic B-spline approximations for numerical solution of fourth order singular boundary value problems. J. Math. 52(3), 47–63 (2020)
Stewart, J.: Single variable essential calculus: early transcendentals. Cengage Learning, (2012)
Mishra, H.K., Saini, S.: Quartic B- spline method for solving a singular singularly perturbed third-order boundary value problems. Am. J. Numer. Anal. 3(1), 18–24 (2015)
Iqbal, M.K., Abbas, M., Wasim, I.: New cubic B-spline approximation for solving third order Emden- Flower type equations. Appl. Math. Comput. 331, 319–333 (2018)
Kim, W., Chun, C.: A modified Adomian decomposition method for solving higher-order singular boundary value problems. Z. Naturforsch. 65(12), 1093–1100 (2010)
Taiwo, O.A., Hassan, M.O.: Approximation of higher-order singular initial and boundary value problems by iterative decomposition and Bernstein polynomial methods. J. Adv. Math. Comput. Sci., pages 498–515, (2015)
Acknowledgements
The authors express their gratitude to the editor for providing valuable assistance throughout the publication process. Furthermore, the comments and recommendations provided by the referees have significantly improved the article.
Funding
No funding was received to assist with the preparation of this manuscript.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no Conflict of interest.
Ethical Approval
Hereby we confirm that article is not under consideration in other journals.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ahmed, H.M., Abd-Elhameed, W.M. Spectral Solutions of Specific Singular Differential Equations Using A Unified Spectral Galerkin-Collocation Algorithm. J Nonlinear Math Phys 31, 42 (2024). https://doi.org/10.1007/s44198-024-00194-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s44198-024-00194-0