Efficient Decomposition Shooting Method for Solving Third-Order Boundary Value Problems

Nawal Al-Zaid, Kholoud Alzahrani, Huda Bakodah^*, Mariam Al-Mazmumy

Department of Mathematics, Faculty of Science, University of Jeddah, Jeddah, Saudi Arabia.

DOI: 10.4236/ijmnta.2023.123006                

Abstract

The present paper proposes a mathematical method to numerically treat a class of third-order linear Boundary Value Problems (BVPs). This method is based on the combination of the Adomian Decomposition Method (ADM) and, the modified shooting method. A complete derivation of the proposed method has been provided, in addition to its numerical implementation and, validation via the utilization of the Runge-Kutta method and, other existing methods. The method has been applied to diverse test problems and turned out to perform remarkably. Lastly, the simulated numerical results have been graphically illustrated and, also supported by some absolute error comparison tables.

Keywords

Linear Third Order BVPs, Shooting Method, Adomian Decomposition Method, Two-Point Boundary Value Problem

1. Introduction

Many physical problems arising in engineering, science, economics, social science and, business are modeled mostly using the Ordinary Differential Equations (ODEs) in the form of either the Boundary Value problems (BVPs) or Initial Value Problems (IVPs), or their combination. Additionally, not all of these ODEs have exact analytical solutions through the known analytical techniques. Thus, this reason has forced various mathematicians and researchers to direct their inquisitiveness in the search for optimal computational procedures that will serve as an alternative approach to at least get hold of approximate solutions whenever the exact solutions are not analytically realistic.

However, the ODEs of particular concern in this study are the third-order linear BVPs. Many attempts have been made in the past literature to study different forms of third-order BVPs, including for instance the coupled systems of third-order BVPs and, laying a solid foundation for the existence and stability of the associated solutions, see [1] - [6] and, the references therewith. Moreover, as the methods for BVPs are very limited, Javeed et al. [7] proposed a numerical approach based on the application of the shooting method to solve different BVPs by coming up with higher-order initial guesses. It is remarkable to state here that the shooting method utilized in [7] is a promising iterative approach that starts off by transforming the governing BVPs to corresponding IVPs. To go further, some coupled systems of third-order BVPs were tackled using different methods like the finite difference method by Noor et al. [8] and, via the application of a meticulous computational procedure by Al-Said [9] . See also the recent work of Nasira et al. [10] for a new numerical procedure to solve third-order BVPs with two-point and multipoint Robin boundary data. Furthermore, the well-known semi-analytical method called the Adomian Decomposition Method (ADM) has equally been greatly used in both the past and recent literature for solving various linear and nonlinear BVPs, see [11] - [20] and the references therewith for some related studies and various modifications of the ADM.

Furthermore, since we aim in this study to combine the ADM and shooting method to numerically examine the third-order linear BVPs, it is notable to mention the work of Attili and Syam [21] where the combination of these methods was proposed and used to study the two-point BVPs of the form

$q\left(t\right)u^{″}\left(t\right)+q^{\prime }\left(t\right)u^{\prime }\left(t\right)=f\left(t,u,u^{\prime }\right),$ (1)

with the following prescribed Dirichlet boundary data

$u\left(a\right)=\alpha ,u\left(b\right)=\beta ,$ (2)

where the functions $q\left(t\right),q^{\prime }\left(t\right)$ and $f\left(t,u,u^{\prime }\right)$ are known continuous functions. Further, the following inverse operator was particularly utilized while studying the problem [21]

$L^{-1}\left(t,D\right)=\underset{a}{\overset{t}{{\displaystyle \int }}}\frac{\text{d}x}{q\left(x\right)}\underset{a}{\overset{x}{{\displaystyle \int }}}\frac{\text{d}s}{\xi \left(s\right)}u\left(s\right),$ (3)

where

$L\left(t,D\right)=\xi \left(t\right)Dq\left(t\right)D,$ (4)

with $D=\frac{d}{d\left(t\right)}$ , and $\frac{1}{q\left(t\right)}$ is locally integrable; while $q\left(t\right)$ and $\xi \left(t\right)$ are smooth functions.

However, the present paper proposes a numerical method to study a class of third-order linear BVPs by utilizing the combination of the ADM and, the modified shooting method. A complete derivation of this method will be provided, executed and, validated using promising methods in the literature.

2. Adomian Decomposition Method (ADM)

The ADM is an efficient semi-analytical method that is widely used to solve a variety of functional equations. To demonstrate how this method is utilized, let us consider the following generalized third-order linear ODE

$y^{‴}=p\left(x\right)y^{″}+q\left(x\right)y^{\prime }+r\left(x\right)y+s\left(x\right),a\le x\le b,$ (5)

subject to the following initial conditions

$y\left(a\right)=\lambda ,y^{\prime }\left(a\right)=\alpha ,y^{″}\left(a\right)=t.$ (6)

Now, expressing the above equation in an operator form, we write

$Ly=p\left(x\right)y^{″}+q\left(x\right)y^{\prime }+r\left(x\right)y+s\left(x\right),a\le x\le b,$ (7)

where $L=\frac{{\text{d}}^3}{\text{d}{x}^3}$ , and the functions $p\left(x\right),q\left(x\right),r\left(x\right),s\left(x\right)$ are known prescribed functions; while $\lambda ,\alpha$ and t are also known constants. The inverse operator is expressed as

$L^{-1}(.)={\displaystyle \underset{a}{\overset{x}{\int }}{\displaystyle \underset{a}{\overset{x}{\int }}{\displaystyle \underset{a}{\overset{x}{\int }}(.)\text{d}x\text{d}x\text{d}x}}}.$ (8)

Therefore, applying the inverse operators given in Equation (8) to Equation (7), we have

$y\left(x\right)=\varphi \left(x\right)+L^{-1}\left(p\left(x\right)y^{″}+q\left(x\right)y^{\prime }+r\left(x\right)y+s\left(x\right)\right),$ (9)

such that

$L\varphi \left(x\right)=0.$ (10)

Additionally, the ADM decomposed the solution $y\left(x\right)$ as follows

$y\left(x\right)={\displaystyle \underset{n=0}{\overset{\infty }{\sum }}{y}_n}\left(x\right),$

such that upon substituting it into Equation (9) yields

${\displaystyle \underset{n=0}{\overset{\infty }{\sum }}{y}_n}=\varphi \left(x\right)+L^{-1}\left[p\left(x\right)\left(\frac{{\text{d}}^2}{\text{d}{x}^2}{\displaystyle \underset{n=0}{\overset{\infty }{\sum }}{y}_n}\right)+q\left(x\right)\left(\frac{\text{d}}{\text{d}x}{\displaystyle \underset{n=0}{\overset{\infty }{\sum }}{y}_n}\right)+r\left(x\right)\left({\displaystyle \underset{n=0}{\overset{\infty }{\sum }}{y}_n}\right)+s\left(x\right)\right],$

through which the ADM recurrently gives the following iterates

$y_0=\varphi \left(x\right)+L^{-1}\left(s\left(x\right)\right),$

$y_{n+1}=L^{-1}\left(p\left(x\right){y^{″}}_n+q\left(x\right){y^{\prime }}_n+r\left(x\right)y_n\right),n\ge 0.$ (11)

such that the overall closed-form solution (approximate solution) of the governing equation given in Equation (5) together with the conditions given in Equation (6) is expressed as

$y_m\left(x\right)={\displaystyle \underset{n=0}{\overset{m}{\sum }}{y}_n}\left(x\right).$

3. Efficient Decomposition Shooting Method (EDSM)

The shooting method is a promising iterative approach that starts off by transforming the governing BVP to a corresponding system of IVPs with specified initial value conditions. These unknown initial conditions are guessed to solve the IVPs; whereas the accurateness of the guessed missing initial conditions is thereafter ascertained by putting both the given value at the terminal point and, the computed value of the dependent variable side by side. Furthermore, when these two values differ, a new value should be guessed continuously until an agreed value is attained between the two under a specified degree of accurateness.

At the moment, let us consider the following third-order linear two-point BVP

$y^{‴}=p\left(x\right)y^{″}+q\left(x\right)y^{\prime }+r\left(x\right)y+s\left(x\right),a\le x\le b,$ (12)

$y\left(a\right)=\lambda ,y^{\prime }\left(a\right)=\alpha ,y\left(b\right)=\beta ,$ (13)

Now, based on the shooting method, the third-order BVP given above will turn into two IVPs where the boundary conditions given in Equation (13) will be replaced with specific initial conditions for each IVP as follows

$u^{‴}=p\left(x\right)u^{″}+q\left(x\right)u^{\prime }+r\left(x\right)u+s\left(x\right),a\le x\le b,$ (14)

$u\left(a\right)=\lambda ,u^{\prime }\left(a\right)=\alpha ,u^{″}\left(a\right)=0,$ (15)

and

$v^{‴}=p\left(x\right)v^{″}+q\left(x\right)v^{\prime }+r\left(x\right)v,a\le x\le b,$ (16)

$v\left(a\right)=0,v^{\prime }\left(a\right)=0,v^{″}\left(a\right)=1.$ (17)

Additionally, Equation (14) can be written as

$Lu=p\left(x\right)u^{″}+q\left(x\right)u^{\prime }+r\left(x\right)u+s\left(x\right),$ (18)

where the proposed differential operator may take the form of L earlier defined; together with inverse $L^{-1}$ given in Equation (8).

By applying $L^{-1}$ on Equation (18), one gets

$u\left(x\right)={\varphi }_1\left(x\right)+L^{-1}\left(p\left(x\right)u^{″}\right)+L^{-1}\left(q\left(x\right)u^{\prime }\right)+L^{-1}\left(r\left(x\right)u\right)+L^{-1}\left(s\left(x\right)\right),$ (19)

such that $L{\varphi }_1\left(x\right)=0$ .

Next, the ADM decomposes the solution $u\left(x\right)$ by the following infinite series

$u\left(x\right)={\displaystyle \underset{n=0}{\overset{\infty }{\sum }}{u}_n}\left(x\right),$ (20)

where the component $u_n\left(x\right)$ are recurrently determined. More, substituting the above series solution into Equation (19) yields

$\begin{array}{c}{\displaystyle \underset{n=0}{\overset{\infty }{\sum }}{u}_n}={\varphi }_1\left(x\right)+L^{-1}\left(p\left(x\right){\displaystyle \underset{n=0}{\overset{\infty }{\sum }}{u^{″}}_n}\left(x\right)\right)+L^{-1}\left(q\left(x\right){\displaystyle \underset{n=0}{\overset{\infty }{\sum }}{u^{\prime }}_n}\left(x\right)\right)\\ +L^{-1}\left(r\left(x\right){\displaystyle \underset{n=0}{\overset{\infty }{\sum }}{u}_n}\left(x\right)\right)+L^{-1}\left(s\left(x\right)\right),\end{array}$ (21)

of which the components $u_n\left(x\right)$ are recurrently found as follows

$\begin{array}{l}{u}_0={\varphi }_1\left(x\right)+L^{-1}\left(s\left(x\right)\right),\\ u_1=L^{-1}\left(p\left(x\right){u^{″}}_0\left(x\right)+q\left(x\right){u^{\prime }}_0\left(x\right)+r\left(x\right)u_0\left(x\right)\right),\\ u_2=L^{-1}\left(p\left(x\right){u^{″}}_1\left(x\right)+q\left(x\right){u^{\prime }}_1\left(x\right)+r\left(x\right)u_1\left(x\right)\right),\\ ⋮\end{array}$ (22)

and so on. Also, for the sake of numerical computation, the m-term approximant are considered as

$u={\psi }_{1,m}\left(x\right)={\displaystyle \underset{k=0}{\overset{m}{\sum }}{u}_k}\left(x\right),$ (23)

Similarly, corresponding to the second BVP given in Equation (16), we can compute

$v={\psi }_{2,m}\left(x\right)={\displaystyle \underset{k=0}{\overset{m}{\sum }}{v}_k}\left(x\right),$ (24)

where the components $v_k\left(x\right)$ can be determined as

$\begin{array}{l}{v}_0={\varphi }_2\left(x\right),\\ v_1=L^{-1}\left(p\left(x\right){v^{″}}_0\left(x\right)+q\left(x\right){v^{\prime }}_0\left(x\right)+r\left(x\right)v_0\left(x\right)\right),\\ v_2=L^{-1}\left(p\left(x\right){v^{″}}_1\left(x\right)+q\left(x\right){v^{\prime }}_1\left(x\right)+r\left(x\right)v_1\left(x\right)\right),\\ ⋮\end{array}$ (25)

and so on. Subsequently, we construct z where z can be written as

$\begin{array}{l}{z}_1\left(x\right)=u_1\left(x\right)+\theta v_1\left(x\right),\\ z_2\left(x\right)=u_2\left(x\right)+\theta v_2\left(x\right),\\ ⋮\end{array}$ (26)

and so on, where $\theta$ is a constant number.

Hence, let $u\left(x\right)$ and $v\left(x\right)$ denote the solutions to the third-order linear IVPs given in Equations (14)-(17), respectively, then, we define

$z\left(x\right)=u\left(x\right)+\frac{\beta -u\left(b\right)}{v\left(b\right)}v\left(x\right),v\left(b\right)\ne 0,$ (27)

where $z\left(x\right)$ is the solution to the third-order linear BVP given Equations ((12), (13)).

4. Numerical Examples

This section examines the proposed methodology for the third-order linear BVPs by demonstrating its applicability to some selected test problems. The method is also compared with the mixture of the fourth-order Runge-Kutta approach and the shooting method to further assess the performance of the proposed method. Similarly, we will compare our results with other numerical methods in the literature [22] - [28] according to the value of m used in them. Therefore, each method used in the comparison was indicated by an abbreviated symbol as follows:

· LSRKM4: Fourth-Order Linear Shooting Runge-Kutta Method in [22] .

· FDM: Finite Difference Method in [22] .

· PAM (3, 3): Pade Approximation (3, 3) Method in [22] .

· RCAM: Rational Chebyshev Approximation Method in [22] .

· QBSM: Quartic B-Spline Method in [22] .

· HPM2 Second-Order Homotopy Perturbation Method [23] .

· LSHPM2: Second-Order Least Square Homotopy Perturbation Method [23] .

· RKT3s4: Three-Stage Fourth-Order Explicit Runge-Kutta Type Method [24] .

· GJGOMM: Generalized Jacobi-Galerkin Operational Matrix Method in [25] .

· OCM: Operational Collocation Method in [26] .

· ESM: Exponential Spline Method [27] .

· QNSM: Quintic Non-polynomial Spline Method [28] .

Further, we present certain supportive Tables 1-8 and Figures 1-6 reporting the absolute error difference between the exact analytical solution and, on the other hand, the obtained numerical results using the proposed Efficient Decomposition Shooting Method (E_EDSM) and, further validated with the fourth-order Runge-Kutta method (E_SRKM4). Furthermore, it is worth mentioning here that, in all the graphical comparative illustrations given in Figures 1-6, the exact analytical solution is portrayed using a blue line, the proposed approximate solution E_EDSM is shown using a dashed dotted red line and, finally, the approximate benchmark solution E_SRKM4 is depicted using asterisked curve black.

Example 1. Consider the third-order linear BVP given by [22]

$y^{‴}=2x^2{y}^{″}-3x{y}^{\prime }-5x^{2}y+{\text{e}}^{2x}\left(3x^3-x^2-5x-4\right),0\le x\le 1,$

$y\left(0\right)=1,y^{\prime }\left(0\right)=1,y\left(1\right)=0.$

The exact analytical solution is given by $y\left(x\right)={\text{e}}^{2x}\left(1-x\right)$ .

Now, based on the proposed modified shooting decomposition method, we make consideration to the following two IVPs

$u^{‴}=2x^2{u}^{″}-3x{u}^{\prime }-5x^{2}u+{\text{e}}^{2x}\left(3x^3-x^2-5x-4\right),$ (28)

$u\left(0\right)=1,u^{\prime }\left(0\right)=1,u^{″}\left(0\right)=0,$ (29)

and

$v^{‴}=2x^2{v}^{″}-3x{v}^{\prime }-5x^{2}v,$ (30)

$v\left(0\right)=0,v^{\prime }\left(0\right)=0,v^{″}\left(0\right)=1.$ (31)

The operator versions of Equations (28) and (30) can be written as

$Lu=2x^2{u}^{″}-3x{u}^{\prime }-5x^{2}u+{\text{e}}^{2x}\left(3x^3-x^2-5x-4\right),$ (32)

$Lv=2x^2{v}^{″}-3x{v}^{\prime }-5x^{2}v.$ (33)

Applying $L^{-1}$ to both sides in Equation (32) using the conditionsin Equation (29), we get

$\begin{array}{c}u\left(x\right)=\frac{15}{4}+\frac{27}{8}x+\frac{17}{16}{x}^2+\frac{3}{8}{x}^3{\text{e}}^{2x}-\frac{29}{16}{x}^2{\text{e}}^{2x}+\frac{25}{8}x{\text{e}}^{2x}\\ -\frac{11}{4}{\text{e}}^{2x}+L^{-1}\left(2x^2{u}^{″}-3x{u}^{\prime }-5x^{2}u\right).\end{array}$ (34)

Similarly, applying $L^{-1}$ to both sides of Equation (33) using the conditions in Equation (31), we get

$v\left(x\right)=\frac{x^2}{2}+L^{-1}\left(2x^2{v}^{″}-3x{v}^{\prime }-5x^{2}v\right).$ (35)

What’s more, decomposing the respective solutions in Equations (34) and (35) based on the ADM, then respective recursive relationships are thus given by

$\{\begin{array}{l}{u}_0\left(x\right)=\frac{15}{4}+\frac{27}{8}x+\frac{17}{16}{x}^2+\frac{3}{8}{x}^3{\text{e}}^{2x}-\frac{29}{16}{x}^2{\text{e}}^{2x}+\frac{25}{8}x{\text{e}}^{2x}-\frac{11}{4}{\text{e}}^{2x},\\ u_{n+1}\left(x\right)=2L^{-1}\left(x^2{u^{″}}_n\right)-3L^{-1}\left(x{u^{\prime }}_n\right)-5L^{-1}\left(x^2{u}_n\right),n\ge 0,\end{array}$

and

$\{\begin{array}{l}{v}_0\left(x\right)=\frac{x^2}{2},\\ v_{n+1}\left(x\right)=2L^{-1}\left(x^2{v^{″}}_n\right)-3L^{-1}\left(x{v^{\prime }}_n\right)-5L^{-1}\left(x^2{v}_n\right),n\ge 0.\end{array}$

In addition, expressing some of the terms from the above recursive relationships, we get

$\{\begin{array}{l}{u}_0\left(x\right)=\frac{15}{4}+\frac{27}{8}x+\frac{17}{16}{x}^2+\frac{3}{8}{x}^3{\text{e}}^{2x}-\frac{29}{16}{x}^2{\text{e}}^{2x}+\frac{25}{8}x{\text{e}}^{2x}-\frac{11}{4}{\text{e}}^{2x},\\ u_1\left(x\right)=-\frac{87}{128}x-\frac{17}{672}{x}^7-\frac{9}{64}{x}^6-\frac{167}{480}{x}^5-\frac{27}{64}{x}^4+\frac{9}{64}{x}^5{\text{e}}^{x^2}-\frac{57}{64}{x}^4{\text{e}}^{x^2}\\ +\frac{73}{32}{x}^3{\text{e}}^{x^2}-\frac{397}{128}{x}^2{\text{e}}^{x^2}+\frac{333}{128}x{\text{e}}^{x^2}-\frac{123}{128}{\text{e}}^{x^2}-\frac{23}{128}{x}^2+\frac{123}{128},\\ ⋮\end{array}$

and

$\{\begin{array}{l}{v}_0\left(x\right)=\frac{x^2}{2},\\ v_1\left(x\right)=-\frac{1}{60}{x}^5-\frac{1}{84}{x}^7,\\ ⋮\end{array}$

Then, the solutions of Equations (28) and (30) with $m=10$ are obtained in a series form. Finally, the approximate solution $z\left(x\right)$ with $m=10$ , $h=\frac{1}{10}$ , and $x_k=kh$ for $k=0,1,\cdots ,m$ , is given by

$z\left(x_k\right)=u\left(x_k\right)+\frac{-u\left(1\right)}{v\left(1\right)}v\left(x_k\right).$

In Table 1, we report the absolute error difference between the exact analytical solution and, the proposed solution E_EDSM and, further validated with E_SRKM4 and the other methods described in [22] . In E_SRKM4 and E_LSRKM4, the same method was used, which is the shooting method with Runge-Kutta method of order four, but using different hypotheses, so that in E_SRKM4, hypotheses (29) and (31) were used, but in E_LSRKM4 other hypotheses were used, which are as follows

$u\left(0\right)=1,u^{\prime }\left(0\right)=0,u^{″}\left(0\right)=0,$

and

$v\left(0\right)=0,v^{\prime }\left(0\right)=1,v^{″}\left(0\right)=0,$

thus, we note that E_SRKM4 was more accurate than E_LSRKM4.

Table 1. The absolute error for EDSM, SRKM4 and the methods described in [22] when $m=10$ .

Table 2. Comparison between different methods when $m=10$ .

From Table 2, we can see that EDSM is the most efficient method for solving example 1, compare with the results of the five methods in ref. [22] and SRKM4.

Again, we portray the exact analytical and approximate solutions in Figure 1; one would notice an ideal agreement between these solutions.

Example 2. Consider the third-order linear BVP given by [23] [24] [25] [26]

$y^{‴}-xy={\text{e}}^x\left(x^3-2x^2-5x-3\right),0\le x\le 1$

$y\left(0\right)=0,y^{\prime }\left(0\right)=1,y\left(1\right)=0.$

The exact analytical solution is given by $y\left(x\right)=x\left(1-x\right){\text{e}}^x$ .

The approximate solution $z\left(x\right)$ with $m=10$ , $h=\frac{1}{10}$ , and $x_k=kh$ for $k=0,1,\cdots ,m$ , is given by

$z\left(x_k\right)=u\left(x_k\right)+\frac{-u\left(1\right)}{v\left(1\right)}v\left(x_k\right).$

In Table 3, we report the absolute error difference between the exact analytical solution and, the proposed solution E_EDSM and, further validated with E_SRKM4 and the other methods described in [23] .

From Table 4, we can see that EDSM is the most efficient method for solving example 2, compare with the results of the all methods in [23] [24] [25] [26] and SRKM4.

Again, we portray the exact analytical and approximate solutions in Figure 2; one would notice an ideal agreement between these solutions.

Example 3. Consider the third-order linear singularly perturbed BVP given by [24] [25] [27] :

$-ϵy^{‴}+y=81{ϵ}^2\cos \left(3x\right)+3ϵ\sin \left(3x\right),0\le x\le 1$

Figure 1. Graphical comparative depictions of the exact and approximate solutions with $m=10$ .

Table 3. The absolute error for EDSM, SRKM4 and the methods described in [23] when $m=10$ .

Table 4. Comparison between different methods.

Figure 2. Graphical comparative depictions of the exact and approximate solutions with $m=10$ .

$y\left(0\right)=0,y^{\prime }\left(0\right)=9ϵ,y\left(1\right)=3ϵ\sin \left(3\right).$

The exact analytical solution is given by $y\left(x\right)=3ϵ\sin \left(3x\right)$ .

The approximate solution $z\left(x\right)$ with $m=10$ , $h=\frac{1}{10}$ , and $x_k=kh$ for $k=0,1,\cdots ,m$ , is given by

$z\left(x_k\right)=u\left(x_k\right)+\frac{3ϵ\sin \left(3\right)-u\left(1\right)}{v\left(1\right)}v\left(x_k\right).$

In Table 5, we report the absolute error difference between the exact analytical solution and, the proposed solution E_EDSM and, further validated with E_SRKM4.

From Table 6, we can see that EDSM is the most efficient method for solving example 3, compare with the results of the all methods in [24] [25] [27] and SRKM4.

Again, we portray the exact analytical and approximate solutions in Figures 3-5; one would notice an ideal agreement between these solutions.

Example 4. Consider the third-order linear BVP given by [24] [27] [28] :

$\begin{array}{l}-ϵy^{‴}+y=6ϵ{\left(1-x\right)}^5{x}^3-6{ϵ}^2(6{\left(1-x\right)}^5-90{\left(1-x\right)}^{4}x\\ +180{\left(1-x\right)}^3{x}^2-60{\left(1-x\right)}^2{x}^3),0\le x\le 1,\end{array}$

$y\left(0\right)=0,y^{\prime }\left(0\right)=0,y\left(1\right)=0.$

The exact analytical solution is given by $y\left(x\right)=6x^3ϵ{\left(1-x\right)}^5$ .

The approximate solution $z\left(x\right)$ with $m=10$ , $h=10$ , and $x_k=kh$ for

$k=0,1,\cdots ,m$ , is given by

Table 5. The absolute error for EDSM and SRKM4 when $m=10$ and different values $ϵ$ .

Table 6. Comparison between different methods when different values m and $ϵ$ .

Figure 3. Graphical comparative depictions of the exact and approximate solutions with $m=10$ and $ϵ=\frac{1}{16}$ .

Figure 4. Graphical comparative depictions of the exact and approximate solutions with $m=10$ and $ϵ=\frac{1}{32}$ .

Figure 5. Graphical comparative depictions of the exact and approximate solutions with $m=10$ and $ϵ=\frac{1}{64}$ .

$z\left(x_k\right)=u\left(x_k\right)+\frac{-u\left(1\right)}{v\left(1\right)}v\left(x_k\right).$

In Table 7, we report the absolute error difference between the exact analytical solution and, the proposed solution E_EDSM and, further validated with E_SRKM4.

From Table 8, we can see that EDSM is the most efficient method for solving example 4, comparison with the results of the all methods in [24] [27] [28] and SRKM4 when $m=10$ . Also, the EDSM was more accurate when compared with the same methods at $m\ge 20$ , where we obtained the absolute error value of zero

Table 7. The absolute error for our method when $m=10$ and different values $ϵ$ .

Table 8. Comparison between different methods when $m=10$ and different values $ϵ$ .

for all values of $ϵ$ .

Again, we portray the exact analytical and approximate solutions in Figures 6-8; one would notice an ideal agreement between these solutions.

5. Conclusion

In conclusion, the present paper proposed a numerical method to treat a particular class of third-order BVPs based on the combination of the shooting method and, the Adomian decomposition method (EDSM). A complete derivation of the method has been provided, in addition to its numerical implementation and, validation with the help of the shooting method with the fourth-order Runge-Kutta method (SRKM4). The proposed method was further applied to certain test problems and turned out to outperform the SRKM4 and, other available

Figure 6. Graphical comparative depictions of the exact and approximate solutions with $m=10$ and $ϵ=\frac{1}{16}$ .

Figure 7. Graphical comparative depictions of the exact and approximate solutions with $m=10$ and $ϵ=\frac{1}{32}$ .

Figure 8. Graphical comparative depictions of the exact and approximate solutions with $m=10$ and $ϵ=\frac{1}{64}$ .

methods in the literature. Lastly, we reported the simulated numerical results via graphical illustrations and, comparison tables.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Tunç, C. (2009) On the Stability and Boundedness of Solutions of Nonlinear Vector Differential Equations of Third Order. Nonlinear Analysis: Theory, Methods & Applications, 70, 2232-2236. https://doi.org/10.1016/j.na.2008.03.002
[2]	Ezeilo, J.O.C. (1967) A Generalization of a Boundedness Theorem for a Certain Third-Order Differential Equation. Mathematical Proceedings of the Cambridge Philosophical Society, 63, 735-742. https://doi.org/10.1017/S0305004100041736
[3]	Ezeilo, J.O.C. (1962) A Property of the Phase-Space Trajectories of a Third Order Non-Linear Differential Equation. Journal of the London Mathematical Society, 37, 33-41. https://doi.org/10.1112/jlms/s1-37.1.33
[4]	Reissig, R., Sansone, G. and Conti, R. (1974) Nonlinear Differential Equations of Higher Order. Noordhoff, Groningen.
[5]	Tunç, C. and Ales, M. (2006) Stability and Boundedness Results for Solutions of Certain Third Order Nonlinear Vector Differential Equations. Nonlinear Dynamics, 45, 273-281. https://doi.org/10.1007/s11071-006-1437-3
[6]	Rauch, L.L. (1950) Oscillation of a Third Order Nonlinear Autonomous System. In: Lefschetz, S., Ed., Contributions to Theory of Nonlinear Oscillations, Vol. 1, Princeton University Press, Princeton, 39-88. https://doi.org/10.1515/981400882632-003
[7]	Javeed, S., Shabnam, A. and Baleanu, D. (2019) An Improved Shooting Technique for Solving Boundary Value Problems Using Higher Order Initial Approximation Algorithms. Punjab University Journal of Mathematics, 51, 101-113.
[8]	Noor, M.A., Al-Said, E. and Noor, K. (2012) Finite Difference Method for Solving a System of Third-Order Boundary Value Problem. Journal of Applied Mathematics, 2012, Article ID: 351764. https://doi.org/10.1155/2012/351764
[9]	Al-Said, E.A. (2000) Numerical Solutions for System of Third-Order Boundary Value Problems. International Journal of Computer Mathematics, 78, 111-121. https://doi.org/10.1080/00207160108805100
[10]	Nasir, N.M., Majid, Z.A., Ismail, F. and Bachok, N. (2021) Direct Integration of the Third-Order Two Point and Multipoint Robin Type Boundary Value Problems. Mathematics and Computers in Simulation, 182, 411-427. https://doi.org/10.1016/j.matcom.2020.10.028
[11]	Adomian, G. (1994) Solving Frontier Problems of Physics: The Decomposition Method. Kluwer, Boston. https://doi.org/10.1007/978-94-015-8289-6
[12]	Adomian, G. (1988) A Review of the Decomposition Method in Applied Mathematics. Journal of Mathematical Analysis and Applications, 135, 501-544. https://doi.org/10.1016/0022-247X(88)90170-9
[13]	Singh, N. and Kumar, M. (2011) Adomian Decomposition Method for Solving Higher Order Boundary Value Problems. Mathematical Theory and Modeling, 2, 11-22.
[14]	Adomain, G. and Rach, R. (1992) Noise Terms in Decomposition Solution Series. Computers & Mathematics with Applications, 24, 61-64. https://doi.org/10.1016/0898-1221(92)90031-C
[15]	Adomain, G. and Rach, R. (1994) Modified Decomposition Solution of Linear and Nonlinear Boundary-Value Problems. Nonlinear Analysis: Theory, Methods & Applications, 23, 615-619. https://doi.org/10.1016/0362-546X(94)90240-2
[16]	Bakodah, H.O. (2012) Some Modifications of Adomian Decomposition Method Applied to Nonlinear System of Fredholm Integral Equations of the Second Kind. International Journal of Contemporary Mathematical Sciences, 7, 929-942.
[17]	Bakodah, H.O. (2013) Modified Adomian Decomposition Method for the Generalized Fifth Order KdV Equations. American Journal of Computational Mathematics, 3, 53-58. https://doi.org/10.4236/ajcm.2013.31008
[18]	Al-Zaid, N.A., Bakodah, H.O. and Hendi, F.A. (2013) Numerical Solutions of the Regularized Long-Wave (RLW) Equation Using New Modification of Laplace-Decomposition Method. Advances in Pure Mathematics, 3, 159-163. https://doi.org/10.4236/apm.2013.31A022
[19]	Al-Zaid, N.A., Bakodah, H.O. and Ebaid, A. (2018) Solving a Class of Partial Differential Equations with Different Types of Boundary Conditions by Using a Generalized Inverse Operator: Decomposition Method. Nonlinear Analysis and Differential Equations, 6, 25-41. https://doi.org/10.12988/nade.2018.843
[20]	Bakodah, H.O., Hendi, F.A. and Al-Zaid, N. (2012) Application of the New Modified Decomposition Method to the Regularized Long-Wave Equation. Life Science Journal, 9, 5862-5866.
[21]	Attili, B.S. and Syam, M.I. (2008) Efficient Shooting Method for Solving Two Point Boundary Value Problems. Chaos, Solitons & Fractals, 35, 895-903. https://doi.org/10.1016/j.chaos.2006.05.094
[22]	Shanab, S. (2017) Numerical Methods for Solving Third Order Two-Point Boundary Value Problems. An-Najah National University, Nablus.
[23]	Qayyum, M. and Oscar, O. (2021) Least Square Homotopy Perturbation Method for Ordinary Differential Equations. Journal of Mathematics, 2021, Article ID: 7059194. https://doi.org/10.1155/2021/7059194
[24]	Abdulsalam, A., Senu, N. and Majid, Z.A. (2019) Direct One-Step Method for Solving Third-Order Boundary Value Problems. International Journal of Applied Mathematics, 32, 155-176. https://doi.org/10.12732/ijam.v32i2.1
[25]	Abd-Elhameed, W.M. (2015) Some Algorithms for Solving Third-Order Boundary Value Problems Using Novel Operational Matrices of Generalized Jacobi Polynomials. Abstract and Applied Analysis, 2015, Article ID: 672703. https://doi.org/10.1155/2015/672703
[26]	Abd-Elhameed, W.M. and Napoli, A. (2020) A Unified Approach for Solving Linear and Nonlinear Odd-Order Two-Point Boundary Value Problems. Bulletin of the Malaysian Mathematical Sciences Society, 43, 2835-2849. https://doi.org/10.1007/s40840-019-00840-7
[27]	Wakjira, Y.A. and Duressa, G.F. (2020) Exponential Spline Method for Singularly Perturbed Third-Order Boundary Value Problems. Demonstratio Mathematica, 53, 360-372. https://doi.org/10.1515/dema-2020-0024
[28]	Wakjira, Y.A., Duressa, G.F. and Bullo, T.A. (2018) Quintic Non-Polynomial Spline Methods for Third Order Singularly Perturbed Boundary Value Problems. Journal of King Saud University-Science, 30, 131-137. https://doi.org/10.1016/j.jksus.2017.01.008