Next Article in Journal
Dark Matter Detection in the Stratosphere
Next Article in Special Issue
Numerical Investigation of Fractional-Order Fornberg–Whitham Equations in the Framework of Aboodh Transformation
Previous Article in Journal
On the Asymptotic Behavior of Class of Third-Order Neutral Differential Equations with Symmetrical and Advanced Argument
Previous Article in Special Issue
Quasi-Monomiality Principle and Certain Properties of Degenerate Hybrid Special Polynomials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Improved Symmetric Numerical Approach for Systems of Second-Order Two-Point BVPs

by
Busyra Latif
1,2,
Md Yushalify Misro
3,
Samsul Ariffin Abdul Karim
4,5,6 and
Ishak Hashim
2,7,*
1
Mathematical Sciences Studies, College of Computing, Informatics and Media, Universiti Teknologi MARA (UiTM) Negeri Sembilan Branch, Seremban 3 Campus, Seremban 70300, Negeri Sembilan, Malaysia
2
Department of Mathematical Sciences, Faculty of Science & Technology, Universiti Kebangsaan Malaysia (UKM), Bangi 43600, Selangor, Malaysia
3
School of Mathematical Sciences, Universiti Sains Malaysia, Gelugor 11800, Pulau Pinang, Malaysia
4
Software Engineering Programme, Faculty of Computing and Informatics, Universiti Malaysia Sabah, Jalan UMS, Kota Kinabalu 88400, Sabah, Malaysia
5
Data Technologies and Applications (DaTA) Research Lab, Faculty of Computing and Informatics, Universiti Malaysia Sabah, Jalan UMS, Kota Kinabalu 88400, Sabah, Malaysia
6
Creative Advanced Machine Intelligence (CAMI) Research Centre, Universiti Malaysia Sabah, Jalan UMS, Kota Kinabalu 88400, Sabah, Malaysia
7
Nonlinear Dynamics Research Center (NDRC), Ajman University, Ajman P.O. Box 346, United Arab Emirates
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(6), 1166; https://doi.org/10.3390/sym15061166
Submission received: 10 April 2023 / Revised: 22 May 2023 / Accepted: 26 May 2023 / Published: 29 May 2023
(This article belongs to the Special Issue Numerical Analysis, Approximation Theory, Differential Equations)

Abstract

:
This study deals with the numerical solution of a class of linear systems of second-order boundary value problems (BVPs) using a new symmetric cubic B-spline method (NCBM). This is a typical cubic B-spline collocation method powered by new approximations for second-order derivatives. The flexibility and high order precision of B-spline functions allow them to approximate the answers. These functions have a symmetrical property. The new second-order approximation plays an important role in producing more accurate results up to a fifth-order accuracy. To verify the proposed method’s accuracy, it is tested on three linear systems of ordinary differential equations with multiple step sizes. The numerical findings by the present method are quite similar to the exact solutions available in the literature. We discovered that when the step size decreased, the computational errors decreased, resulting in better precision. In addition, details of maximum errors are investigated. Moreover, simple implementation and straightforward computations are the main advantages of the offered method. This method yields improved results, even if it does not require using free parameters. Thus, it can be concluded that the offered scheme is reliable and efficient.

1. Introduction

Most problems arising from scientific and engineering applications, especially applications in geodesics, are boundary value problems (BVPs), which are much more difficult to solve than initial value problems (IVPs). Since it is generally difficult to find closed-form solutions for BVPs, many researchers have attempted to develop methods to find approximate or numerical solutions for BVPs. Well-known methods involve the shooting method [1], finite difference methods [2,3,4] and spectral methods [5,6,7]. In some real-life situations, the shooting method produces numerically sensitive systems of algebraic equations, which must be solved using other numerical methods [8].
In the present paper, we consider the following system of linear two-point second-order BVPs:
w ( x ) + p 1 ( x ) w ( x ) + p 2 ( x ) w ( x ) + p 3 ( x ) z ( x ) + p 4 ( x ) z ( x ) + p 5 ( x ) z ( x ) = f 1 ( x ) , z ( x ) + q 1 ( x ) z ( x ) + q 2 ( x ) z ( x ) + q 3 ( x ) w ( x ) + q 4 ( x ) w ( x ) + q 5 ( x ) w ( x ) = f 2 ( x ) ,
with boundary conditions
w ( 0 ) = w ( 1 ) = z ( 0 ) = z ( 1 ) = 0 ,
where 0 x 1 . In particular, f 1 ( x ) and f 2 ( x ) are given functions, and p i and q i with i = 1 , 2 , 3 , 4 , 5 are continuous and sufficiently smooth functions on the interval I = [ 0 , 1 ] . Theorems that systematically list the existence and uniqueness of the problem solutions of (1) and (2) have been studied in [9]. In recent times, applications of linear and non-linear systems of two-point boundary value problems can be found in economics, biology, physics and mathematics. For instance, Nikooeinejad et al. [10] obtained the approximate solution of two-point BVPS for four applications of differential games in economics and management science using a combined numerical algorithm. In biology, the Shortley–Weller scheme has been implemented for a two-point boundary value problem. This numerical scheme was later applied to investigate tumor growth problems in heterogeneous microenvironments [11]. On the other hand, the application of two-point boundary value problems has been addressed in the problem of calculating rocket trajectories in the atmosphere [12].
Several researchers have investigated the linear and non-linear systems of second-order boundary value problems and produced various efficient and accurate numerical methods. These methods include the Laplace homotopy analysis [13,14], continuous genetic algorithm method [15], sinc collocation method [16,17], He’s homotopy perturbation method [18], reproducing kernel space method [19], multistage optimal homotopy asymptotic method [20], variational iteration method (VIM) [21] and Chebyshev finite difference method [22].
Researchers have been interested in the families of B-splines for their potential to approximate the solution of BVPs accurately and efficiently. B-spline methods have several attractive features and flexibility that make them useful in numerical computation to solve BVPs [23]. For example, the B-spline is the smoothest interpolation function compared to other piecewise polynomial interpolation functions [24]. Moreover, B-splines have small local support properties. In recent years, the cubic B-spline collocation method captured the attention of some researchers to solve partial differential equations [25], fractional differential equations [26], fractional partial differential equations [27], etc.
This study focuses on finding the solutions of two-point BVPs using the cubic B-spline method. Bickley was the first to explore cubic splines to approximate the solutions of two-point BVPs [28]. Later, Albasiny and Hoskins enhanced Bickley’s work by solving the two-point BVPs using a tri-diagonal matrix of coefficients [29]. Since then, several researchers have earned more interest in employing spline functions for solving BVPs [30,31,32,33]. Caglar et al. in [34] evaluated the two-point BVPs solutions using the cubic B-spline basis function. Hamid et al. [35,36] considered the ECBM and cubic trigonometric B-spline method for the solution of linear two-point BVPs. Apart from that, Heilat and Ismail [37] used a hybrid cubic B-spline method to evaluate the solutions of non-linear two-point boundary value problems. Recently, a hybrid cubic B-spline method with an optimized parameter was used by Heilat et al. [38] to solve linear two-point BVPs.
The linear system of second-order BVPs has gained attention from Caglar and Caglar [39] and Heilat et al. [40]. They represented the cubic B-spline method (CBM) and ECBM, respectively. The ECBM involved two parameters in boosting the flexibility of the spline curve. Based on the investigation, the ECBM is the best compared to the CBM, He’s homotopy perturbation method [18], Laplace homotopy analysis method [13], reproducing kernel [19] and sinc-collocation method [17]. In recent years, Zhang and Niu [41] found the approximate solution of second-order BVPs using a Lobatto-reproducing kernel and declared the method has high precision accuracy in different spaces.
The new symmetric cubic B-spline method (NCBM) was first studied by Lang and Xu in [42] to solve non-linear second-order BVPs with two dependent variables. The NCBM is a typical CBM, equipped with a new second-order derivative approximation. Then, Iqbal et al. [43] explored the NCBM for solving several third-order Emden–Flower type equations. A year after that, Wasim et al. [44] extended the NCBM and proposed the new extended cubic B-spline method (NECBM) for solving the class of second-order singular BVPs. Moreover, the nonlinear third-order Korteweg–de Vries equations were solved by Abbas et al. [45] using the NCBM to approximate the solutions. Later, Nazir et al. [46] improved the method to a new quintic B-spline approximation technique as a method to approximate the numerical solution of the Boussinesq equation. Recently, Nazir et al. [47] implemented the NCBM for the numerical solutions of coupled viscous Burgers equations.
Thus, motivated by all these works, we aim to figure out whether the proposed method, the NCBM, can perform much better in solving the linear system of two-point second-order BVPs. The rest of this paper is as follows. In Section 2, the typical definition of cubic B-spline basis functions is described. Then, Section 3 presents the descriptions of the numerical method. The convergence analysis of the method has been proven in Section 4. The numerical results and their discussion are summarized in Section 5. Finally, the paper ends in Section 6 with a brief conclusion.

2. Preliminary Concepts

This section describes the classical cubic B-spline approximation and the new second-order approximation invented by Lang and Xu [42]. Let the finite interval [ a , b ] , where a = x 0 < < x N = b is divided into uniform partitions with a mesh point x i = x 0 + i h , i = 0 ( 1 ) N using a step size h = b a N , N Z + . The typical cubic B-spline basis function is defined as [34].
B i ( x ) = ( x x i 2 ) 3 , x [ x i 2 , x i 1 ] , h 3 + 3 h 2 ( x x i 1 ) + 3 h ( x x i 1 ) 2 3 ( x x i 1 ) 3 , x [ x i 1 , x i ] , h 3 + 3 h 2 ( x i + 1 x ) + 3 h ( x i + 1 x ) 2 3 ( x i + 1 x ) 3 , x [ x i , x i + 1 ] , ( x i + 2 x ) 3 , x [ x i + 1 , x i + 2 ] , 0 , otherwise ,
where i = 1 ( 1 ) N + 1 . The cubic B-spline holds the geometric invariability, convex hull property and symmetry [48]. For a sufficiently smooth function w ( x ) and z ( x ) , there exist a unique third-degree spline W ( x ) and Z ( x ) that satisfies the prescribed interpolating conditions given by
W ( x i ) = w ( x i ) , i = 0 , 1 , 2 , . . . , N , W ( a ) = w ( a ) , W ( b ) = w ( b ) , W ( a ) = w ( a ) , W ( b ) = w ( b ) ,
and
Z ( x i ) = z ( x i ) , i = 0 , 1 , 2 , . . . , N , Z ( a ) = z ( a ) , Z ( b ) = z ( b ) , Z ( a ) = z ( a ) , Z ( b ) = z ( b ) ,
in which
W ( x ) = i = 1 N + 1 σ i B i ( x ) ,
Z ( x ) = i = 1 N + 1 η i B i ( x ) ,
where σ i and η i are unknown real coefficients to be computed. The values of B i ( x ) and the first and second derivatives B i ( x ) and B i ( x ) at mesh point x i are tabulated in Table 1. From (4), (5) and Table 1, the cubic B-spline approximations W ( x j ) , W ( x j ) , Z ( x j ) and Z ( x j ) can be simplified as follows:
W j = i = j 1 j + 1 σ i B i ( x ) = 1 6 ( σ j i + 4 σ j + σ j + 1 ) ,
s j ( x ) = i = j 1 j + 1 σ i B i ( x ) = 1 2 h ( σ j i + σ j + 1 ) ,
Z j ( x ) = i = j 1 j + 1 η i B i ( x ) = 1 6 ( η j i + 4 η j + η j + 1 ) ,
r j ( x ) = i = j 1 j + 1 η i B i ( x ) = 1 2 h ( η j i + η j + 1 ) .
The second derivatives, W ( x j ) and Z ( x j ) can be simplified as S j and R j , respectively. Subsequently, the new approximation for second-order derivatives can be represented as follows [42,49]:
S j = 1 12 h 2 14 σ j 1 33 σ j + 28 σ j + 1 14 σ j + 2 + 6 σ j + 3 σ j + 4 , for j = 0 , σ j 2 + 8 σ j 1 18 σ j + 8 σ j + 1 + σ j + 2 , for j = 1 , , N 1 , σ n 4 + 6 σ n 3 14 σ n 2 + 28 σ n 1 33 σ n + 14 σ n + 1 , for j = N ,
and
R j = 1 12 h 2 14 η j 1 33 η j + 28 η j + 1 14 η j + 2 + 6 η j + 3 η j + 4 , for j = 0 , η j 2 + 8 η j 1 18 η j + 8 η j + 1 + η j + 2 , for j = 1 , , N 1 , η n 4 + 6 η n 3 14 η n 2 + 28 η n 1 33 η n + 14 η n + 1 , for j = N .
We note that this NCBM has up to a fifth-order accuracy [49].

3. Implementation of the Method

In this section, we extended Caglar’s work [39] by solving the linear system two-point second-order BVP and adopted the new second-order approximation.
Discretizing (1) at the knot x j gives the following expression:
S k + 1 ( x j ) + p 1 ( x j ) s k + 1 ( x j ) + p 2 ( x j ) W k + 1 ( x j ) + p 3 ( x j ) R k + 1 ( x j ) + p 4 ( x j ) r k + 1 ( x j ) + p 5 ( x j ) Z k + 1 ( x j ) = f 1 k ( x j ) ,
R k + 1 ( x j ) + q 1 ( x j ) r k + 1 ( x j ) + q 2 ( x j ) Z k + 1 ( x j ) + q 3 ( x j ) S k + 1 ( x j ) + q 4 ( x j ) s k + 1 ( x j ) + q 5 ( x j ) W k + 1 ( x j ) = f 2 k ( x j ) ,
where k = 0 , 1 , 2 , . By substituting (6)–(11) into (12) for j = 0 , 1 , 2 , , N 1 , N , we obtain the following equation:
For j = 0 ,
14 σ 1 33 σ 0 + 28 σ 1 14 σ 2 + 6 σ 3 σ 4 12 h 2 + p 1 ( x 0 ) σ 1 + σ 1 2 h + p 2 ( x 0 ) σ 1 + 4 σ 0 + σ 1 6 + p 3 ( x 0 ) 14 η 1 33 η 0 + 28 η 1 14 η 2 + 6 η 3 η 4 12 h 2 + p 4 ( x 0 ) η 1 + η 1 2 h + p 5 ( x 0 ) η 1 + 4 η 0 + η 1 6 = f k 1 ( x 0 ) .
For j = 1 , 2 , , N 1 ,
σ j 2 + 8 σ j 1 18 σ j + 8 σ j + 1 + σ j + 2 12 h 2 + p 1 ( x j ) σ j 1 + σ j + 1 2 h + p 2 ( x j ) σ j 1 + 4 σ j + σ j + 1 6 + p 3 ( x j ) η j 2 + 8 η j 1 18 η j + 8 η j + 1 + η j + 2 12 h 2 + p 4 ( x j ) η j 1 + η j + 1 2 h + p 5 ( x j ) η j 1 + 4 η j + η j + 1 6 = f k 1 ( x j ) .
For j = N ,
σ N 4 + 6 σ N 3 14 σ N 2 + 28 σ N 1 33 σ N + 14 σ N + 1 12 h 2 + p 1 ( x N ) σ N 1 + σ N + 1 2 h + p 2 ( x N ) σ N 1 + 4 σ N + σ N + 1 6 + p 4 ( x N ) η N 4 + 6 η N 3 14 η N 2 + 28 η N 1 33 η N + 14 η N + 1 12 h 2 + p 4 ( x N ) η N 1 + η N + 1 2 h + p 5 ( x N ) η N 1 + 4 η N + η N + 1 6 = f k 1 ( x N ) .
By substituting (6)–(11) into (13) for j = 0 , 1 , 2 , , N 1 , N , we obtain the following equation:
For j = 0 ,
14 η 1 33 η 0 + 28 η 1 14 η 2 + 6 η 3 η 4 12 h 2 + q 1 ( x 0 ) η 1 + η 1 2 h + q 2 ( x 0 ) η 1 + 4 η 0 + η 1 6 + q 3 ( x 0 ) 14 σ 1 33 σ 0 + 28 σ 1 14 σ 2 + 6 σ 3 σ 4 12 h 2 + q 4 ( x 0 ) σ 1 + σ 1 2 h + q 5 ( x 0 ) σ 1 + 4 σ 0 + σ 1 6 = f k 1 ( x 0 ) .
For j = 1 , 2 , , N 1 ,
η j 2 + 8 η j 1 18 η j + 8 η j + 1 + η j + 2 12 h 2 + q 1 ( x j ) η j 1 + η j + 1 2 h + q 2 ( x j ) η j 1 + 4 η j + η j + 1 6 + q 3 ( x j ) σ j 2 + 8 σ j 1 18 σ j + 8 σ j + 1 + σ j + 2 12 h 2 + q 4 ( x j ) σ j 1 + σ j + 1 2 h + q 5 ( x j ) σ j 1 + 4 σ j + σ j + 1 6 = f k 1 ( x j ) .
For j = N ,
η N 4 + 6 η N 3 14 η N 2 + 28 η N 1 33 η N + 14 η N + 1 12 h 2 + q 1 ( x N ) η N 1 + η N + 1 2 h + q 2 ( x N ) η N 1 + 4 η N + η N + 1 6 + q 3 ( x N ) σ N 4 + 6 σ N 3 14 σ N 2 + 28 σ N 1 33 σ N + 14 σ N + 1 12 h 2 + q 4 ( x N ) σ N 1 + σ N + 1 2 h + q 5 ( x N ) σ N 1 + 4 σ N + σ N + 1 6 = f k 1 ( x N ) .
Consequently, we have 2 N + 2 linear equations involving 2 N + 6 unknowns. Thus, we need four additional equations, which can be obtained from the boundary conditions in (2) below:
σ 1 + 4 σ 0 + σ 1 = 0 , σ N 1 + 4 σ N + σ N + 1 = 0 , η 1 + 4 η 0 + η 1 = 0 , η N 1 + 4 η N + η N + 1 = 0 .
Hence, the above system will have the ( 2 N + 6 ) × ( 2 N + 6 ) dimensional matrix form that can be expressed as:
A X = B ,
where matrix A is given by:
A = A 1 | A 2 A 4 | A 3 ,
X = [ σ 1 , σ 0 , σ 1 , , σ N 1 , σ N , σ N + 1 , η 1 , η 0 , η 1 , , η N 1 , η N , η N + 1 ] T
and
B = [ 0 , 12 h 2 f 1 ( x 0 ) , 12 h 2 f 1 ( x 1 ) , , 12 h 2 f 1 ( x N ) , 0 , 0 , 12 h 2 f 2 ( x 0 ) , 12 h 2 f 2 ( x 1 ) , , 12 h 2 f 2 ( x N ) , 0 ) ] T .
The four sub-matrices A 1 , A 2 , A 3 and A 4 are represented as follows:
A 1 = 1 4 1 0 0 0 0 0 0 0 0 0 ϵ 1 ϵ 2 ϵ 3 ϵ 4 ϵ 5 ϵ 6 0 0 0 0 0 0 α 1 , 1 α 2 , 1 α 3 , 1 α 4 , 1 α 5 , 1 0 0 0 0 0 0 0 0 α 1 , 2 α 2 , 2 α 3 , 2 α 4 , 2 α 5 , 2 0 0 0 0 0 0 0 0 0 0 0 0 α 1 , N 2 α 2 , N 2 α 3 , N 2 α 4 , N 2 α 5 , N 2 0 0 0 0 0 0 0 0 α 1 , N 1 α 2 , N 1 α 3 , N 1 α 4 , N 1 α 5 , N 1 0 0 0 0 0 0 ϕ 1 ϕ 2 ϕ 3 ϕ 4 ϕ 5 ϕ 6 0 0 0 0 0 0 0 0 0 1 4 1 ,
A 2 = 0 0 0 0 0 0 0 0 0 0 0 0 ϵ 7 ϵ 8 ϵ 9 ϵ 10 ϵ 11 ϵ 12 0 0 0 0 0 0 α 6 , 1 α 7 , 1 α 8 , 1 α 9 , 1 α 10 , 1 0 0 0 0 0 0 0 0 α 6 , 2 α 7 , 2 α 8 , 2 α 9 , 2 α 10 , 2 0 0 0 0 0 0 0 0 0 0 0 0 α 6 , N 2 α 7 , N 2 α 8 , N 2 α 9 , N 2 α 10 , N 2 0 0 0 0 0 0 0 0 α 6 , N 1 α 7 , N 1 α 8 , N 1 α 9 , N 1 α 10 , N 1 0 0 0 0 0 0 ϕ 7 ϕ 8 ϕ 9 ϕ 10 ϕ 11 ϕ 12 0 0 0 0 0 0 0 0 0 0 0 0 ,
A 3 = 1 4 1 0 0 0 0 0 0 0 0 0 β 1 β 2 β 3 β 4 β 5 β 6 0 0 0 0 0 0 γ 1 , 1 γ 2 , 1 γ 3 , 1 γ 4 , 1 γ 5 , 1 0 0 0 0 0 0 0 0 γ 1 , 2 γ 2 , 2 γ 3 , 2 γ 4 , 2 γ 5 , 2 0 0 0 0 0 0 0 0 0 0 0 0 γ 1 , N 2 γ 2 , N 2 γ 3 , N 2 γ 4 , N 2 γ 5 , N 2 0 0 0 0 0 0 0 0 γ 1 , N 1 γ 2 , N 1 γ 3 , N 1 γ 4 , N 1 γ 5 , N 1 0 0 0 0 0 0 δ 1 δ 2 δ 3 δ 4 δ 5 δ 6 0 0 0 0 0 0 0 0 0 1 4 1 ,
and
A 4 = 0 0 0 0 0 0 0 0 0 0 0 0 β 7 β 8 β 9 β 10 β 11 β 12 0 0 0 0 0 0 γ 6 , 1 γ 7 , 1 γ 8 , 1 γ 9 , 1 γ 10 , 1 0 0 0 0 0 0 0 0 γ 6 , 2 γ 7 , 2 γ 8 , 2 γ 9 , 2 γ 10 , 2 0 0 0 0 0 0 0 0 0 0 0 0 γ 6 , N 2 γ 7 , N 2 γ 8 , N 2 γ 9 , N 2 γ 10 , N 2 0 0 0 0 0 0 0 0 γ 6 , N 1 γ 7 , N 1 γ 8 , N 1 γ 9 , N 1 γ 10 , N 1 0 0 0 0 0 0 δ 7 δ 8 δ 9 δ 10 δ 11 δ 12 0 0 0 0 0 0 0 0 0 0 0 0 .
For j = 0 ,
ϵ 1 = 14 6 p 1 ( x 0 ) h + 2 p 2 ( x 0 ) h 2 , ϵ 7 = 14 p 3 ( x 0 ) 6 p 4 ( x 0 ) + 2 p 5 ( x 0 ) h 2 , ϵ 2 = 33 + 8 p 2 ( x 0 ) h 2 , ϵ 8 = 33 p 3 ( x 0 ) + 8 p 5 ( x 0 ) h 2 , ϵ 3 = 28 + 6 p 1 ( x 0 ) h + 2 p 2 ( x 0 ) h 2 , ϵ 9 = 28 p 3 ( x 0 ) + 6 p 4 ( x 0 ) h + 2 p 5 ( x 0 ) h 2 , ϵ 4 = 14 , ϵ 10 = 14 p 3 ( x 0 ) , ϵ 5 = 6 , ϵ 11 = 6 p 3 ( x 0 ) , ϵ 6 = 1 , ϵ 12 = p 3 ( x 0 ) ,
β 1 = 14 6 q 1 ( x 0 ) h + 2 q 2 ( x 0 ) h 2 , β 7 = 14 q 3 ( x 0 ) 6 q 4 ( x 0 ) h + 2 q 5 ( x 0 ) h 2 , β 2 = 33 + 8 q 2 ( x 0 ) h 2 ; β 8 = 33 q 3 ( x 0 ) + 8 q 5 ( x 0 ) h 2 , β 3 = 28 + 6 q 1 ( x 0 ) h + 2 q 2 ( x 0 ) h 2 , β 9 = 28 q 3 ( x 0 ) + 6 q 4 ( x 0 ) h + 2 q 5 ( x 0 ) h 2 , β 4 = 14 , β 10 = 14 q 3 ( x 0 ) , β 5 = 6 , β 11 = 6 q 3 ( x 0 ) , β 6 = 1 , β 12 = q 3 ( x 0 ) .
For j = 1 , 2 , , N 1 ,
α 1 , j = 1 , α 6 , j = p 3 ( x j ) , α 2 , j = 8 6 p 1 ( x j ) h + 2 p 2 ( x j ) h 2 , α 7 , j = 8 p 3 ( x j ) 6 p 4 ( x j ) h + 2 p 5 ( x j ) h 2 , α 3 , j = 18 + 8 p 2 ( x j ) h 2 , α 8 , j = 18 p 3 ( x j ) + 8 p 5 ( x j ) h 2 , α 4 , j = 8 + 6 p 1 ( x j ) h + 2 p 2 ( x j ) h 2 , α 9 , j = 8 p 3 ( x j ) + 6 p 4 ( x j ) h + 2 p 5 ( x j ) h 2 , α 5 , j = 1 , α 10 , j = p 3 ( x j ) ,
γ 1 , j = 1 , γ 6 , j = q 3 ( x j ) , γ 2 , j = 8 6 q 1 ( x j ) h + 2 q 2 h 2 , γ 7 , j = 8 q 3 ( x j ) 6 q 4 ( x j ) h + 2 q 5 ( x j ) h 2 , γ 3 , j = 18 + 8 q 2 ( x j ) h 2 , γ 8 , j = 18 q 3 ( x j ) + 8 q 5 ( x j ) h 2 , γ 4 , j = 8 + 6 q ( x j ) h + 2 q 2 ( x j ) h 2 , γ 9 , j = 8 q 3 ( x j ) + 6 q 4 ( x j ) h + 2 q 5 ( x j ) h 2 , γ 5 , j = 1 , γ 10 , j = q 3 ( x j ) .
For j = N ,
ϕ 1 = 1 , ϕ 7 = p 3 ( x N ) , ϕ 2 = 6 , ϕ 8 = 6 p 3 ( x N ) , ϕ 3 = 14 , ϕ 9 = 14 p 3 ( x N ) , ϕ 4 = 28 6 p 1 ( x N ) h + 2 p 2 ( x N ) h 2 , ϕ 10 = 28 p 3 ( x N ) 6 p 4 ( x N ) h + 2 p 5 ( x N ) h 2 , ϕ 5 = 33 + 8 p 2 ( x N ) h 2 , ϕ 11 = 33 p 3 ( x N ) + 8 p 5 ( x N ) h 2 , ϕ 6 = 14 + 6 p 1 ( x N ) h + 2 p 2 ( x N ) h 2 , ϕ 12 = 14 p 3 ( x N ) h + 6 p 4 ( x N ) h + 2 p 5 ( x N ) h 2 ,
δ 1 = 1 , δ 7 = q 3 ( x N ) , δ 2 = 6 , δ 8 = 6 q 3 ( x N ) , δ 3 = 14 , δ 9 = 14 q 3 ( x N ) , δ 4 = 28 6 q 1 ( x N ) h + 2 q 2 ( x N ) h 2 , δ 10 = 28 q 3 ( x N ) 6 q 4 ( x N ) h + 2 q 5 ( x N ) h 2 , δ 5 = 33 + 8 q 2 ( x N ) h 2 , δ 11 = 33 q 3 ( x N ) + 8 q 5 ( x N ) h 2 , δ 6 = 14 + 6 q 1 ( x N ) h + 2 q 2 ( x N ) h 2 , δ 12 = 14 q 3 ( x N ) + 6 q 4 ( x N ) h + 2 q 5 ( x N ) h 2 .
Since the matrix A is a banded matrix, the system of linear equations is solved using a generalization of the Thomas algorithm. This method has been proposed in [50]. Matlab R2018a running on an Intel(R) CORE(TM) i7-1165G7 CPU 1.30 GHz processor, 8.00 GB RAM, was used to execute the numerical computations.

4. Convergence Analysis

In this section, we will prove the order of convergence of our method.
Theorem 1.
Let p i ( x ) C 2 [ 0 , 1 ] , where i = 1 , 2 , 3 , 4 , 5 are continuous and sufficiently smooth functions. Then, let w ˜ be the known exact solution of the boundary value problems (1), (2) and also m ˜ be the cubic B-spline approximation to w ˜ . Thus, the uniform error is stated by
w ˜ m ˜ ϕ h 2 .
Proof. 
Let w ˜ be the exact solution of the boundary value problems (1), (2) and m ˜ be the cubic B-spline approximation to w ˜ given by:
w ˜ = m ˜ = i = 1 N + 1 σ i ˜ B i ( x ) i
where
σ ˜ = σ i ˜ = [ σ ˜ 1 , σ ˜ 0 , , σ ˜ N + 1 ] T .
Furthermore, suppose m ^ ( x ) is the computed cubic B-spline approximation to m ˜ ( x ) given by
w ^ ( x i ) = m ^ ( x i ) = i = 1 N + 1 σ i ^ B i ( x ) i
where
σ ^ = σ i ^ = [ σ ^ 1 , σ ^ 0 , , σ ^ N + 1 ] T .
To approximate the error
w ˜ ( x i ) m ˜ ( x i ) ,
we need to estimate the error
w ˜ ( x i ) m ^ ( x i )
and
w ^ ( x i ) m ˜ ( x i )
differently. We know that the system of the ( n + 3 ) × ( n + 3 ) matrix leads to
A σ = F .
It follows that
A σ ˜ = F ˜
and
A σ ^ = F ^ .
Then, by substracting (25) and (26), we have
A ( σ ^ σ ˜ ) = F ^ F ˜ ,
where A is an ( n + 3 ) × ( n + 3 ) dimensional matrix, and
F = [ F 1 , F 0 , , F N + 1 ] T ,
where T stands for transpose. Hence, from (27), we have
( σ ^ σ ˜ ) = A 1 ( F ^ F ˜ ) .
Now, consider taking the infinity norm from (29), and we have
σ ^ σ ˜ = A 1 F ^ F ˜ .
Note that the B-spline basis B 1 , B 0 , B 1 , , B N + 1 satisfies the following inequality
i = 1 N + 1 σ i ^ B i ( x ) i 1 .
Adopted from [51,52,53], we have
A 1 F ^ F ˜ ϕ h 2 ,
σ ^ σ ˜ ϕ h 2 .
Additionally,
m ^ ( x i ) m ˜ ( x i ) = ( σ ^ σ ˜ ) i = 1 N + 1 B i ( x ) i ,
m ^ ( x i ) m ˜ ( x i ) = ( σ ^ σ ˜ ) i = 1 N + 1 B i ( x ) i .
Now, consider
m ^ ( x i ) m ˜ ( x i ) ( σ ^ σ ˜ ) i = 1 N + 1 B i ( x ) i ϕ h 2 ,
w ˜ ( x i ) m ^ ( x i ) ρ h 4
and
w ˜ ( x i ) m ˜ ( x i ) w ˜ ( x i ) m ^ ( x i ) + m ^ ( x i ) m ˜ ( x i ) .
Substituting (35) and (36) into (37), we have
w ˜ ( x i ) m ˜ ( x i ) ϕ h 2 + ρ h 4 = α h 2 ,
where α = ϕ + ρ h 2 . □
Thus, this method is second-order convergent, given by
w ˜ ( x i ) m ˜ ( x i ) α h 2 .

5. Numerical Examples

In this section, three numerical problems of the linear system of ordinary differential equations are compared with the exact solutions and existing methods to demonstrate the efficiency and accuracy of the proposed method. The numerical errors are measured using the error norm L , defined as follows:
L = max j W ( x j ) w ( x j ) or L = max j Z ( x j ) z ( x j ) .

5.1. Problem 1

Consider the system of the linear two-point BVPs equation [40]
w ( x ) + x w ( x ) + x z ( x ) = 2 , z ( x ) + 2 x z ( x ) + 2 x w ( x ) = 2 ,
with boundary conditions
w ( 0 ) = w ( 1 ) = z ( 0 ) = 0 = z ( 1 ) = 0 .
The exact solutions are w ( x ) = x 2 x and z ( x ) = x x 2 .
In Table 2, the approximate solutions, the exact solutions and the absolute errors for Problem 1 when N = 5 are reported. It clearly shows that the approximate solutions are in good agreement with the exact solution. Table 3 lists the comparison of error norm between ECBM and the proposed method for Problem 1 with different step sizes, N = 5 and N = 21 . Two free parameters that are involved in ECBM for this case in Table 3 are obtained by trial and error. The truncation error of the proposed method was O ( h 5 ) accurate. On the other hand, CBM and ECBM in [40] were O ( h 2 ) accurate. Our presented method produced more accurate results compared to the earlier methods.

5.2. Problem 2

Consider the system of the linear two-point BVPs equation [40]
w ( x ) + w ( x ) + x w ( x ) + z ( x ) + 2 x z ( x ) = f 1 ( x ) , z ( x ) + z ( x ) + 2 w ( x ) + x 2 w ( x ) = f 2 ( x ) ,
with boundary conditions
w ( 0 ) = w ( 1 ) = z ( 0 ) = 0 = z ( 1 ) = 0 ,
where x [ 0 , 1 ] , f 1 ( x ) = 2 ( x + 1 ) cos ( x ) + π cos ( π x ) + 2 x sin ( π x ) + ( 4 x 2 x 2 4 ) sin ( x ) and f 2 ( x ) = 4 ( x 1 ) cos ( x ) 2 ( 2 x 2 + x 3 ) sin ( x ) ( π 2 1 ) sin ( π x ) . The true solutions are w ( x ) = 2 ( 1 x ) sin ( x ) and z ( x ) = sin ( π x ) .
Table 4 lists the approximate solution, the exact solution and the absolute errors for Problem 2 when N = 5 . It clearly shows that the approximate solutions promise a good agreement with the exact solution.
The comparison of the absolute errors between the present methods and the proposed method are shown for Problem 2 in Table 5 and Table 6. The reproducing kernel method in [19] was solved using eleven points in [ 0 , 1 ] , while the sinc-collocation method in [17] was solved for N = 5 with the same number of points in [ 0 , 1 ] . The error bounds for the reproducing kernel method are at least O ( h 2 ) and O ( h ) . Two trial and error-free parameters involved in ECBM for this case are λ 1 = 1.0 × 10 3 and λ 2 = 1.0 × 10 3 , respectively. Evidently, our proposed method produced better approximations compared with the earlier methods.
In Table 7 and Table 8, the comparison of error norm between ECBM and the presented method are tabulated for Problem 2 when N = 5 and N = 25 , respectively. This clearly shows our presented method is more powerful. Two free parameters that are involved in ECBM for this case in Table 7 and Table 8 are obtained from the optimization technique. Table 9 reports the L error norm with different N for Problem 2.

5.3. Problem 3

Consider the system of the linear two-point BVPs equation [21]
w ( x ) + ( 2 x 1 ) w ( x ) + cos ( π x ) z ( x ) = f 1 ( x ) , z ( x ) + x w ( x ) = f 2 ( x ) ,
with boundary conditions
w ( 0 ) = w ( 1 ) = z ( 0 ) = z ( 1 ) = 0 ,
where x ( 0 , 1 ) , f 1 ( x ) = π 2 sin ( π x ) + ( 2 x 1 ) π cos ( π x ) + ( 2 x 1 ) cos ( π x ) and f 2 ( x ) = 2 + x sin ( π x ) . The analytical solutions are w ( x ) = sin ( π x ) and z ( x ) = x 2 x .
Table 10 lists the approximate, the exact solution and the absolute error when N = 5 . It clearly shows that the approximate solutions exhibit a good agreement with the exact solution.
The comparison of the absolute errors between the existing methods and the proposed method is shown for Problem 3 when N = 41 and N = 20 in Table 11 and Table 12. From both tables, we noted that for N = 20 , NCBM can already match the accuracy of the VIM, CBM and ECBM. Two free parameters involved in ECBM in Table 11 and Table 12 are obtained from trial and error. The VIM in [21] was solved using one iteration step. It is observed that our proposed method is more precise compared to all earlier methods.
Table 13 compares the error norm between He’s Homotopy Perturbation, Laplace Homotopy, ECBM and the presented method for Problem 3 when N = 5 , while Table 14 compares the error norm between ECBM and the presented method for N = 41 . Two free parameters that are involved in ECBM in Table 13 are obtained from the optimization technique. Conversely, trial and error are applied to find two free parameters in Table 14. The L error norm with different N for Problem 3 are tabulated in Table 15. It is noted that our presented method is found to be reasonably good.
Table 2, Table 4, Table 5 and Table 6 and Table 10, Table 11 and Table 12 list the numerical results, the exact solution and the absolute errors for each problem at uniform mesh. Table 9 and Table 15 present the L norm with different N for Problem 2 and Problem 3, respectively. Additionally, for each problem, the details of error norm L of the existing and the proposed methods at different values of N are reported in Table 3, Table 7, Table 8, Table 13 and Table 14. Consequently, the approximation obtained by the proposed method was more precise compared to others. A larger N offers greater precision but at the cost of a longer computation. This method also does not require a free parameter, but it is still the most superior and reliable method compared to the stated existing methods.

6. Conclusions

The NCBM has been applied and analyzed to numerically solve a linear system of two-point boundary value problems in this study. The method presented was based on a typical cubic B-spline, a CBS basis function that engages with the new approximation for the second-order derivative. Theoretically, it has been discovered that our method is second-order convergence. Three numerical examples were presented, and error norms, L , were calculated. We found that as the step size decreased, the error decreased, resulting in higher accuracy. Thus, it is concluded that our method gives comparable results to the stated existing methods. This method has the following three advantages: (a) it can avoid the unnecessary calculation in finding the unknown parameter; (b) it can produce up to a fifth-order accuracy; and (c) it can solve the linear system of a two-point boundary value problem accurately and efficiently. In the future, the proposed method can be used to solve more difficult problems in engineering and sciences through a graphics processing unit, GPU.

Author Contributions

Conceptualization, B.L., S.A.A.K. and I.H.; methodology, B.L., M.Y.M., S.A.A.K. and I.H.; software, B.L., M.Y.M. and S.A.A.K.; validation, B.L., M.Y.M., S.A.A.K. and I.H.; formal analysis, B.L., M.Y.M., S.A.A.K. and I.H.; investigation, B.L.; writing the original draft preparation, B.L.; writing the review and editing, B.L., S.A.A.K. and I.H.; supervision, S.A.A.K. and I.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by I.H.’s UKM grant number DIP-2021-018.

Data Availability Statement

The authors confirm that the data supporting the findings of this study are available within the article.

Acknowledgments

The authors gratefully acknowledge the anonymous referees for their insightful comments and suggestions that helped to enhance this paper. B. Latif is thankful for the full scholarship (SLAB) awarded by the Malaysian Ministry of Higher Education (MOHE) and Universiti Teknologi MARA (UiTM) Malaysia.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BVPsBoundary value-problems
NCBMNew symmetric cubic B-spline method
IVPsInitial value-problems
CBSCubic B-spline
CBMCubic B-spline method
ECBMExtended cubic B-spline method
VIMVariational iteration method

References

  1. Ha, S.N. A nonlinear shooting method for two-point boundary value problems. Comput. Math. Appl. 2001, 42, 1411–1420. [Google Scholar] [CrossRef]
  2. Chawla, M.M. An eighth order tridiagonal finite difference method for nonlinear two-point boundary value problems. Bit Numer. Math. 1977, 17, 281–285. [Google Scholar] [CrossRef]
  3. Usmani, R.A. A method of high-order accuracy for the numerical integration of boundary value problems. BIT 1973, 13, 458–469. [Google Scholar] [CrossRef]
  4. Tirmizi, I.A.; Twizelle, E.H. Higher-order finite-difference methods for nonlinear second-order two-point boundary-value problems. Appl. Math. Lett. 1977, 15, 897–902. [Google Scholar] [CrossRef]
  5. Boyd, J.P. Chebyshev and Fourier Spectral Methods, 2nd ed.; Courier Corporation: Mineola, NY, USA, 2000. [Google Scholar]
  6. Mehrpouya, M.A.; Shamsi, M. Gauss pseudospectral and continuation methods for solving two-point boundary value problems in optimal control theory. Appl. Math. Model. 2015, 39, 5047–5057. [Google Scholar] [CrossRef]
  7. Gheorghiu, C.I. Chebfun solutions to a class of 1D singular and nonlinear boundary value problems. Computation 2022, 10, 116. [Google Scholar] [CrossRef]
  8. Betts, J.T. Survey of numerical methods for trajectory optimization. J. Guid. Control. Dyn. 1998, 21, 193–207. [Google Scholar] [CrossRef]
  9. Wang, H.Z.; Li, Y. Existence and uniqueness of solutions to two-point boundary value problems for ordinary differential equations. Z. Angew. Math. Phys. ZAMP 1996, 47, 373–384. [Google Scholar] [CrossRef]
  10. 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. 2022, 200, 404–427. [Google Scholar] [CrossRef]
  11. Sweidan, M.; Chen, X.; Zheng, X. The Shortley–Weller scheme for variable coefficient two-point boundary value problems and its application to tumor growth problem with heterogeneous microenvironment. J. Comput. Appl. Math. 2020, 376, 112874. [Google Scholar] [CrossRef]
  12. Campos, L.M.B.C.; Gil, P.J.S. The two-point boundary-value problem for rocket trajectories. Aerospace 2020, 7, 131. [Google Scholar] [CrossRef]
  13. Ogunlaran, O.M.; Ademola, A.T. On the Laplace homotopy analysis method for a non-linear system of second-order boundary value problems. Gen. Math. Notes 2015, 26, 11–22. [Google Scholar]
  14. Sami Bataineh, A.; Noorani, M.S.M.; Hashim, I. Modified homotopy analysis method for solving systems of second-order BVPs. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 430–442. [Google Scholar] [CrossRef]
  15. Arqub, O.A.; Abo-Hammour, Z. Numerical solution of systems of second-order boundary value problems using continuous genetic algorithm. Inf. Sci. 2014, 279, 396–415. [Google Scholar] [CrossRef]
  16. El-Gamel, M. Sinc-collocation method for solving linear and non-linear system of second-order boundary value problems. Appl. Math. 2012, 3, 1627–1633. [Google Scholar] [CrossRef]
  17. Dehghan, M.; Saadatmandi, A. The numerical solution of a non-linear system of second-order boundary value problems using the sinc-collocation method. Math. Comput. Model. 2007, 46, 1434–1441. [Google Scholar] [CrossRef]
  18. Saadatmandi, A.; Dehghan, M.; Eftekhari, A. Application of H e’s homotopy perturbation method for non-linear system of second-order boundary value problems. Nonlinear Anal. Real World Appl. 2009, 10, 1912–1922. [Google Scholar] [CrossRef]
  19. Geng, F.; Cui, M. Solving a non-linear system of second order boundary value problems. J. Math. Anal. Appl. 2007, 327, 1167–1181. [Google Scholar] [CrossRef]
  20. Ratib Anakira, N.; Alomari, A.K.; Jameel, A.F.; Hashim, I. Multistage optimal homotopy asymptotic method for solving initial–value problems. J. Nonlinear Sci. Appl. 2016, 9, 1826–1843. [Google Scholar] [CrossRef]
  21. Lu, J. Variational iteration method for solving a non-linear system of second-order boundary value problems. Comput. Math. Appl. 2007, 54, 1133–1138. [Google Scholar] [CrossRef]
  22. Saadatmandi, A.; Askari Farsangi, J. Chebyshev finite difference method for a non-linear system of second-order boundary value problems. Appl. Math. Comput. 2007, 192, 586–591. [Google Scholar]
  23. Akram, T.; Abbas, M.; Iqbal, A.; Baleanu, D.; Asad, H.H. Novel numerical approach based on modified extended cubic B–spline functions for solving non-linear time-fractional telegraph equation. Symmetry 2020, 12, 1154. [Google Scholar] [CrossRef]
  24. Tayebi, S.; Momani, S.; Abu Arqub, O. The cubic B-spline interpolation method for numerical point solutions of conformable boundary value problems. Alex. Eng. J. 2022, 61, 1519–1528. [Google Scholar] [CrossRef]
  25. Parumasur, N.; Adetona, R.A.; Singh, P. Efficient solution of burgers’, modified burgers’ and KdV–burgers’ equations using B-spline approximation functions. Mathematics 2023, 11, 1847. [Google Scholar] [CrossRef]
  26. Abbas, M.; Bibi, A.; Alzaidi, A.S.M.; Nazir, T.; Majeed, A.; Akram, G. Numerical solutions of third-order time-fractional differential equations using cubic B-Spline functions. Fractal Fract. 2022, 6, 528. [Google Scholar] [CrossRef]
  27. Arqub, O.A.; Tayebi, S.; Baleanu, D.; Osman, M.S.; Mahmoud, W.; Alsulam, H. A numerical combined algorithm in cubic B-spline method and finite difference technique for the time-fractional nonlinear diffusion wave equation with reaction and damping terms. Results Phys. 2022, 6, 105912. [Google Scholar] [CrossRef]
  28. Sakai, M. Piecewise cubic interpolation and two-point boundary value problems. Publ. Res. Inst. Math. 1971, 7, 345–362. [Google Scholar] [CrossRef]
  29. Albasiny, E.L.; Hoskins, W.D. Cubic spline solutions to two-point boundary value problems. Comput. J. 1969, 12, 151–153. [Google Scholar] [CrossRef]
  30. Al-Said, E.A. Cubic spline method for solving two-point boundary-value problems. Korean J. Comput. Appl. Math. 1998, 5, 669–680. [Google Scholar] [CrossRef]
  31. Müllenheim, G. Solving two-point boundary value problems with spline functions. Ima J. Numer. Anal. 1992, 12, 503–518. [Google Scholar] [CrossRef]
  32. Khalifa, A.K.A.; Eilbeck, J.C. Collocation with quadratic and cubic splines. Ima J. Numer. Anal. 1982, 2, 111–121. [Google Scholar] [CrossRef]
  33. Al-Said, E.A.; Noor, M.A.; Al-Shejari, A.A. Numerical solutions for system of second-order boundary value problems. Korean J. Comput. Appl. Math. 1998, 5, 659–667. [Google Scholar] [CrossRef]
  34. Caglar, H.; Caglar, N.; Elfaituri, K. B-spline interpolation compared with finite difference, finite element and finite volume methods which applied to two-point boundary value problems. Appl. Math. Comput. 2006, 175, 72–79. [Google Scholar] [CrossRef]
  35. Hamid, N.N.A.; Majid, A.A.; Ismail, A.I.M. Extended cubic B-spline method for linear two-point boundary value problems. Sains Malays. 2011, 40, 1285–1290. [Google Scholar]
  36. Hamid, N.N.A.; Majid, A.A.; Izani, A. Extended cubic B-spline interpolation method applied to linear two-point boundary value problems. World Acad. Sci. Eng. Technol. 2010, 62, 566–568. [Google Scholar]
  37. Salem Heilat, A.; Ismail, A.I.M. Hybrid cubic B-spline method for solving non-linear two-point boundary value problems. Int. J. Pure Appl. Math. 2016, 110, 369–381. [Google Scholar]
  38. Heilat, A.; Zureigat, H.; Batiha, B. New spline method for solving linear two-point boundary value problems. Eur. J. Pure Appl. Math. 2021, 14, 1283–1294. [Google Scholar] [CrossRef]
  39. Caglar, N.; Caglar, H. B-spline method for solving linear system of second-order boundary value problems. Comput. Math. Appl. 2009, 57, 757–762. [Google Scholar] [CrossRef]
  40. Heilat, A.S.; Hamid, N.A.A.; Ismail, A.I.M. Extended cubic B-spline method for solving a linear system of second-order boundary value problems. SpringerPlus 2016, 5, 1–18. [Google Scholar] [CrossRef]
  41. Zhang, J.; Niu, J. Lobatto-reproducing kernel method for solving a linear system of second order boundary value problems. J. Appl. Math. Comput. 2016, 68, 3631–3653. [Google Scholar] [CrossRef]
  42. Lang, F.-G.; Xu, X.-P. A new cubic B-spline method for approximating the solution of a class of non-linear second-order boundary value problem with two dependent variables. Sci. Asia 2014, 40, 444–450. [Google Scholar] [CrossRef]
  43. Kashif Iqbal, M.; Abbas, M.; Wasim, I. New cubic B-spline approximation for solving third order Emden–Flower type equations. Appl. Math. Comput. 2018, 331, 319–333. [Google Scholar]
  44. Wasim, I.; Abbas, M.; Kashif Iqbal, M. A new extended B-spline approximation technique for second order singular boundary value problems arising in physiology. J. Math. Comput. Sci. 2019, 19, 258–267. [Google Scholar] [CrossRef]
  45. Abbas, M.; Iqbal, M.K.; Zafar, B.; Zin, S.B.M. New cubic B-spline approximations for solving non-linear third-order Korteweg-De Vries equation. Indian J. Sci. Technol. 2019, 12, 1–9. [Google Scholar] [CrossRef]
  46. Nazir, T.; Abbas, M.; Iqbal, M.K. A new quintic B-spline approximation for numerical treatment of Boussinesq equation. J. Math. Comput. Sci. 2020, 20, 30–42. [Google Scholar] [CrossRef]
  47. Nazir, T.; Abbas, M.; Iqbal, M.K. New cubic B-spline approximation technique for numerical solutions of coupled viscous Burgers equations. Eng. Comput. 2020, 38, 83–106. [Google Scholar] [CrossRef]
  48. Goh, J.; Majid, A.A.; Ismail, A.I.M. Extended cubic uniform B-spline for a class of singular boundary value problems. Scienceasia 2011, 37, 79–82. [Google Scholar] [CrossRef]
  49. Latif, B.; Karim, S.A.A.; Hashim, I. New cubic B-spline approximation for solving linear two-point boundary-value problems. Mathematics 2021, 9, 1250. [Google Scholar] [CrossRef]
  50. Tolmachev, D.; Chertovskih, R.; Zheligovsky, V. Algorithmic aspects of simulation of magnetic field generation by thermal convection in a plane layer of fluid. Mathematics 2023, 11, 808. [Google Scholar] [CrossRef]
  51. De Boor, C. On the convergence of odd-degree spline interpolation. J. Approx. Theory 1968, 1, 452–463. [Google Scholar] [CrossRef]
  52. Hall, C.A. On error bounds for spline interpolation. J. Approx. Theory 1968, 1, 209–218. [Google Scholar] [CrossRef]
  53. Prenter, P.M. Splines and Variational Methods; Courier Corporation: Dover, NY, USA, 2008. [Google Scholar]
Table 1. Coefficient of B i ( x ) , B i ( x ) and B i ( x ) at the nodes.
Table 1. Coefficient of B i ( x ) , B i ( x ) and B i ( x ) at the nodes.
x i 2 x i 1 x i x i + 1 x i + 2
B i ( x ) 01410
B i ( x ) 0 1 / 2 h 0 1 / 2 h 0
B i ( x ) 0 1 / h 2 2 / h 2 1 / h 2 0
Table 2. Absolute errors for Problem 1 when N = 5 .
Table 2. Absolute errors for Problem 1 when N = 5 .
xNCBMExactAbsolute ErrorNCBMExactAbsolute Error
W ( x ) Solution w ( x ) W ( x ) w ( x ) Z ( x ) Solution z ( x ) Z ( x ) z ( x )
0.2 0.16 0.16 2.78 × 10 17 0.160.16 2.22 × 10 16
0.4 0.24 0.24 1.11 × 10 16 0.240.24 5.55 × 10 17
0.6 0.24 0.24 1.67 × 10 16 0.240.24 1.11 × 10 16
0.8 0.16 0.16 1.39 × 10 16 0.160.16 2.78 × 10 17
Table 3. The L error norm for Problem 1 when N = 5 and N = 21 .
Table 3. The L error norm for Problem 1 when N = 5 and N = 21 .
ECBM [40]NCBMECBM [40]ECBM [40]NCBM
N = 5 N = 5 N = 21 N = 21 N = 21
λ 1 = λ 2 = 0 λ 1 = λ 2 = 0 λ 1 = λ 2 = 1.25 × 10 14
W ( x ) 3.47 × 10 15 1.67 × 10 16 3.72 × 10 13 1.73 × 10 13 1.07 × 10 15
Z ( x ) 3.69 × 10 15 2.22 × 10 16 2.53 × 10 13 1.67 × 10 13 6.94 × 10 16
Table 4. Absolute errors for Problem 2 when N = 5 .
Table 4. Absolute errors for Problem 2 when N = 5 .
xNCBMExactAbsolute ErrorNCBMExactAbsolute Error
W ( x ) Solution w ( x ) W ( x ) w ( x ) Z ( x ) Solution z ( x ) Z ( x ) z ( x )
0.20.3177940.317871 7.71 × 10 5 0.5877780.587785 6.92 × 10 6
0.40.4672740.467302 2.81 × 10 5 0.9509010.951057 1.55 × 10 4
0.60.4517760.451714 6.19 × 10 5 0.9509010.951057 1.56 × 10 4
0.80.2870360.286942 9.37 × 10 5 0.5877780.587785 6.98 × 10 6
Table 5. Absolute errors for Problem 2 for w ( x ) .
Table 5. Absolute errors for Problem 2 for w ( x ) .
x ECBM [40]ECBM [40]
ReproducingSinc− λ 1 = 0 λ 1 = 1.0 × 10 3 NCBM
Kernel [19]Collocation [17] λ 2 = 0 λ 2 = 1.0 × 10 3 N = 25
N = 5 N = 25 N = 25
0.08 3.3 × 10 3 3.2 × 10 3 1.3 × 10 4 1.4 × 10 5 9.6 × 10 8
0.24 7.7 × 10 3 9.2 × 10 4 2.7 × 10 4 1.1 × 10 5 1.6 × 10 7
0.4 9.7 × 10 3 2.0 × 10 3 2.7 × 10 4 2.1 × 10 5 9.2 × 10 8
0.56 9.5 × 10 3 2.2 × 10 4 2.0 × 10 4 5.9 × 10 5 1.9 × 10 8
0.72 7.3 × 10 3 4.1 × 10 3 9.4 × 10 5 7.8 × 10 5 1.0 × 10 7
0.88 3.4 × 10 3 1.0 × 10 2 1.6 × 10 5 5.6 × 10 5 9.3 × 10 8
0.96 1.1 × 10 3 2.1 × 10 3 3.6 × 10 8 2.3 × 10 5 4.0 × 10 8
Table 6. Absolute errors for Problem 2 for z ( x ) .
Table 6. Absolute errors for Problem 2 for z ( x ) .
x ECBM [40]ECBM [40]
ReproducingSinc− λ 1 = 0 λ 1 = 1.0 × 10 3 NCBM
Kernel [19]Collocation [17] λ 2 = 0 λ 2 = 1.0 × 10 3 N = 25
N = 5 N = 25 N = 25
0.08 7.7 × 10 3 1.5 × 10 3 3.8 × 10 4 2.2 × 10 4 1.8 × 10 7
0.24 2.0 × 10 2 7.0 × 10 3 9.9 × 10 4 6.0 × 10 4 4.7 × 10 7
0.4 2.7 × 10 2 7.4 × 10 3 1.3 × 10 3 8.3 × 10 4 6.3 × 10 7
0.56 2.7 × 10 2 1.0 × 10 2 1.4 × 10 3 8.6 × 10 4 6.4 × 10 7
0.72 2.0 × 10 2 4.4 × 10 3 1.1 × 10 3 6.8 × 10 4 5.1 × 10 7
0.88 9.4 × 10 3 2.1 × 10 2 5.0 × 10 4 3.3 × 10 4 2.5 × 10 7
0.96 3.1 × 10 3 6.9 × 10 3 1.7 × 10 4 1.1 × 10 4 8.5 × 10 8
Table 7. The L error norm for Problem 2 when N = 5 .
Table 7. The L error norm for Problem 2 when N = 5 .
ECBM [40]ECBM [40]NCBM
λ 1 = λ 2 = 0 λ 1 = 1.269208 × 10 2
λ 2 = 6.634523 × 10 2
W ( x ) 2.09 × 10 3 1.80 × 10 5 9.37 × 10 5
Z ( x ) 1.75 × 10 4 1.75 × 10 4 1.56 × 10 4
Table 8. The L error norm for Problem 2 when N = 25 .
Table 8. The L error norm for Problem 2 when N = 25 .
ECBM [40]ECBM [40]NCBM
λ 1 = λ 2 = 0 λ 1 = 1.0 × 10 3
λ 2 = 1.0 × 10 3
W ( x ) 2.72 × 10 4 7.80 × 10 5 1.56 × 10 7
Z ( x ) 1.36 × 10 3 8.60 × 10 4 6.53 × 10 7
Table 9. The L error norm with different N for Problem 2.
Table 9. The L error norm with different N for Problem 2.
NNCBM
W ( x )
NCBM
Z ( x )
40 2.40 × 10 8 1.00 × 10 7
80 1.50 × 10 9 6.27 × 10 9
100 6.14 × 10 10 2.57 × 10 9
Table 10. Absolute errors for Problem 3 when N = 5 .
Table 10. Absolute errors for Problem 3 when N = 5 .
xNCBMExactAbsolute ErrorNCBMExactAbsolute Error
W ( x ) Solution w ( x ) W ( x ) w ( x ) Z ( x ) Solution z ( x ) Z ( x ) z ( x )
0.20.5878450.587785 5.98 × 10 5 0.160001 0.160000 1.10 × 10 6
0.40.9509590.951057 9.72 × 10 5 0.240002 0.240000 2.47 × 10 6
0.60.9509590.951057 9.72 × 10 5 0.240002 0.240000 2.37 × 10 6
0.80.5878450.587785 5.97 × 10 5 0.160000 0.160000 4.86 × 10 7
Table 11. Absolute errors for Problem 3 for w ( x ) .
Table 11. Absolute errors for Problem 3 for w ( x ) .
x ECBM [40]
VIMCBMECBM [40] λ 1 = 1.0 × 10 3 NCBMNCBM
[21][40] λ 1 = 0 λ 2 = 0 N = 41 N = 20
N = 41 λ 2 = 0 N = 41
0.1 3.30 × 10 4 1.40 × 10 4 1.30 × 10 4 2.83 × 10 6 1.63 × 10 8 2.83 × 10 7
0.2 2.51 × 10 3 2.80 × 10 4 2.56 × 10 4 5.55 × 10 6 3.57 × 10 8 6.34 × 10 7
0.3 7.84 × 10 3 3.90 × 10 4 3.60 × 10 4 7.81 × 10 6 5.39 × 10 8 9.59 × 10 7
0.4 1.66 × 10 2 4.60 × 10 4 4.28 × 10 4 9.30 × 10 6 6.71 × 10 8 1.19 × 10 6
0.5 2.77 × 10 2 4.80 × 10 2 4.52 × 10 4 9.82 × 10 6 7.26 × 10 8 1.27 × 10 6
0.6 3.87 × 10 2 4.60 × 10 2 4.28 × 10 4 9.30 × 10 6 6.93 × 10 8 1.19 × 10 6
0.7 4.59 × 10 2 3.90 × 10 2 3.60 × 10 4 7.81 × 10 6 5.78 × 10 8 9.59 × 10 7
0.8 4.49 × 10 2 2.80 × 10 2 2.56 × 10 4 5.56 × 10 6 4.05 × 10 8 6.34 × 10 7
0.9 3.09 × 10 2 1.50 × 10 2 1.30 × 10 4 2.83 × 10 6 2.10 × 10 8 6.83 × 10 7
Table 12. Absolute errors for Problem 3 for z ( x ) .
Table 12. Absolute errors for Problem 3 for z ( x ) .
x ECBM [40]ECBM [40]
CBM λ 1 = 0 λ 1 = 1.0 × 10 3 NCBMNCBM
[40] λ 2 = 0 λ 2 = 0 N = 41 N = 20
N = 41 N = 41 N = 41
0.1 5.74 × 10 6 5.74 × 10 6 1.25 × 10 7 8.64 × 10 10 1.54 × 10 8
0.2 1.13 × 10 5 1.13 × 10 5 2.46 × 10 7 1.71 × 10 9 3.05 × 10 8
0.3 1.64 × 10 5 1.64 × 10 5 3.56 × 10 7 2.49 × 10 9 4.42 × 10 8
0.4 2.03 × 10 5 2.03 × 10 5 4.42 × 10 7 3.11 × 10 9 5.51 × 10 8
0.5 2.26 × 10 5 2.26 × 10 5 4.91 × 10 7 3.49 × 10 9 6.12 × 10 8
0.6 2.26 × 10 5 2.26 × 10 5 4.92 × 10 7 3.53 × 10 9 6.10 × 10 8
0.7 2.01 × 10 5 2.01 × 10 5 4.37 × 10 7 3.20 × 10 9 5.39 × 10 8
0.8 1.51 × 10 5 1.51 × 10 5 3.29 × 10 7 2.49 × 10 9 4.01 × 10 8
0.9 8.14 × 10 6 8.14 × 10 6 1.76 × 10 7 1.48 × 10 9 2.13 × 10 8
Table 13. The L error norm for Problem 3 when N = 5 .
Table 13. The L error norm for Problem 3 when N = 5 .
ECBM [40]ECBM [40]
He’s HomotopyLaplace λ 1 = 0 λ 1 = 6.639145 × 10 2 NCBM
Pertubation [18]Homotopy [13] λ 2 = 0 λ 2 = 1.161882 × 10 6
W ( x ) 2.1 × 10 4 2.2 × 10 5 2.8 × 10 2 1.4 × 10 4 9.7 × 10 5
Z ( x ) 3.2 × 10 4 1.1 × 10 5 1.4 × 10 3 7.2 × 10 6 2.5 × 10 6
Table 14. The L error norm for Problem 3 when N = 41 .
Table 14. The L error norm for Problem 3 when N = 41 .
ECBM [40]ECBM [40]
λ 1 = 0 λ 1 = 1.0 × 10 3 NCBM
λ 2 = 0 λ 2 = 0
W ( x ) 4.52 × 10 4 9.82 × 10 6 7.26 × 10 8
Z ( x ) 2.26 × 10 5 4.92 × 10 7 3.56 × 10 9
Table 15. The L error norm with different N for Problem 3.
Table 15. The L error norm with different N for Problem 3.
NNCBMNCBM
W ( x ) Z ( x )
60 1.59 × 10 8 7.77 × 10 10
80 5.02 × 10 9 2.46 × 10 10
100 2.06 × 10 9 1.01 × 10 10
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Latif, B.; Misro, M.Y.; Abdul Karim, S.A.; Hashim, I. An Improved Symmetric Numerical Approach for Systems of Second-Order Two-Point BVPs. Symmetry 2023, 15, 1166. https://doi.org/10.3390/sym15061166

AMA Style

Latif B, Misro MY, Abdul Karim SA, Hashim I. An Improved Symmetric Numerical Approach for Systems of Second-Order Two-Point BVPs. Symmetry. 2023; 15(6):1166. https://doi.org/10.3390/sym15061166

Chicago/Turabian Style

Latif, Busyra, Md Yushalify Misro, Samsul Ariffin Abdul Karim, and Ishak Hashim. 2023. "An Improved Symmetric Numerical Approach for Systems of Second-Order Two-Point BVPs" Symmetry 15, no. 6: 1166. https://doi.org/10.3390/sym15061166

APA Style

Latif, B., Misro, M. Y., Abdul Karim, S. A., & Hashim, I. (2023). An Improved Symmetric Numerical Approach for Systems of Second-Order Two-Point BVPs. Symmetry, 15(6), 1166. https://doi.org/10.3390/sym15061166

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop