New Cubic B-Spline Approximation for Solving Linear Two-Point Boundary-Value Problems

: In this study, we introduce a new cubic B-spline (CBS) approximation method to solve linear two-point boundary value problems (BVPs). This method is based on cubic B-spline basis functions with a new approximation for the second-order derivative. The theoretical new approximation for a second-order derivative and the error analysis have been successfully derived. We found that the second-order new approximation was ( ) 3 O h accurate. By using this new second-order approximation, the proposed method was ( ) 5 O h accurate. Four numerical problems consisting of linear ordinary differential equations and trigonometric equations with different step sizes were performed to validate the accuracy of the proposed methods. The numerical results were compared with the least squares method, finite difference method, finite element method, finite volume method, B-spline interpolation method, extended cubic B-spline interpolation method and the exact solutions. By finding the maximum errors, the results consistently showed that the proposed method gave the best approximations among the existing methods. We also found that our proposed method involved simple implementation and straightforward computations. Hence, based on the results and the efficiency of our method, we can say that our method is reliable and a promising method for solving linear two-point


Introduction
Boundary value problems (BVPs) have been extensively investigated in the fields of physics, chemistry and engineering.Numerous methods have been implemented throughout the years to approximate the solutions of linear and nonlinear two-point BVPs, such as the variational approach, finite difference (FDM), finite element (FEM), finite volume (FVM) and shooting (LSM) [1][2][3].Bickley first employed cubic splines to solve a simple two-point BVP [4], followed by Albasiny and Hoskins [5].Fyfe continued Bickley's work and established the fourth-order accuracy of the cubic spline interpolation scheme [6].Since then, spline functions have been frequently used for solving BVPs [7][8][9][10].
In 2006, Caglar et al. [11] replaced the cubic spline with a cubic B-spline basis function to solve two-point BVPs and named it the cubic B-spline interpolation method (BSI).Several cubic B-spline based numerical approaches have been widely applied to solve linear and nonlinear BVPs [12][13][14][15].The extended cubic B-spline and cubic trigonometric B-spline were studied by Hamid et al. [16,17] as solutions to linear two-point BVPs.It was found that the cubic trigonometric B-spline provided better results compared with the cubic BSI method if the problems had trigonometric functions.Heilat and Ismail [18] developed a hybrid cubic B-spline method to deal with nonlinear two-point boundary value problems.
This study focuses on numerically solving linear two-point BVPs using a new cubic B-spline (CBS) approximation.These functions are flexible enough to approximate the solution at any point of the domain with high accuracy while maintaining a high degree of smoothness at the knots.In recent years, the new CBS approximation method has been applied for solving BVPs.Iqbal et al. [19] investigated the numerical solutions of secondorder singular BVPs using a new approximation for the second-order derivative.A new approximation for the second-order derivative of the extended cubic B-spline basis was developed by Wasim et al. [20] for solving second-order singular BVPs.Iqbal et al. [21] employed a new CBS approximation for solving several third-order Emden-Flower type equations.Nazir et al. [22] used the new CBS approximation for the numerical solutions of coupled viscous Burgers equations.Abbas et al. [23] proposed a new CBS approximation to approximate the solutions to nonlinear third-order Korteweg-de Vries equations.A new quintic B-spline approximation was studied by Nazir et al. [24] as a method for solving Boussinesq equations.Motivated by these works, this study was conducted to solve linear second-order ordinary differential equations (ODEs) based on a new CBS approximation to figure out whether this method would perform equally well.This paper is organized as follows.In Section 2, we present the typical definition of cubic B-spline basis functions.The formulation of a new approximation to second-order derivatives is developed in Section 3. The numerical method is presented in Section 4. In Section 5, the error analysis is investigated.The results and discussion are demonstrated in Sections 6 and 7, respectively.Section 8 summarizes this paper with a brief conclusion.

Cubic B-Spline Functions
Let there be a finite interval ( ) where For a sufficiently smooth function ( ) y x , there exists a unique third-degree spline ( ) Y x that satisfies the prescribed interpolating conditions: where i σ is the set of unknown real coefficients to be computed.The values of and the first and second derivatives ( ) ( ) ( ) 1) and ( 2) and Table 1, the following relations can be obtained: ( ) ( ) ( ) ( ) Using the relations from Equations (3-5), the following expressions can be established [6]:

S y x hy x h y x h y x
Now consider  j S as the new approximation for ( ) where 1 Z , 2 Z and 3 Z are the three parameters to be chosen so that the error order of  j S is as high as possible.From the linear combination, we have By solving Equations (9-11) simultaneously, we get 1 . Using Equation ( 5) and the parameters' values, Equation ( 8) can be expressed as ( ) Next, at the knot 0 x , the linear combination of four neighboring knots can be written as

S Z S Z S Z S Z S
where

S y x hy x h y x h y x
The linear combination in Equation ( 13) produces the following four expressions: . By substituting the parameter values, Equation ( 13) takes the following form: ( ) Similarly, the approximation to ( ) '' n y x can be written as ( )

Description of the Numerical Method
In this section, the numerical scheme for solving linear second-order ordinary differential equations using the new CBS approximation is discussed.Let us consider the general form of a linear second order of a two-point boundary value problem [26]: with boundary conditions where ≤ ≤ a x b , 1 η and 2 η are constants.In addition, ( ) continuous and sufficiently smooth functions on the interval  [26].Assume the CBS solution to Equation ( 15) is given as where σ i is the constants to be calculated.

x h p x p x f x h
In addition, we need two additional equations, which can be obtained from the boundary conditions of Equation ( 16) as shown below: Hence, this will give the following ( ) ( ) and for 1,2, , 1, ( ) The unknown column vector c is determined by a well-known Thomas algorithm, and its components are plugged into Equation ( 17) to obtain the CBS approximation.The solutions to the tridiagonal linear system of equations given in Equation ( 17) are guaranteed to be unique.This is because the matrix is strictly dominant and nonsingular.
Let μ E be the shift operator; in other words, ( ) can be expressed as [6] ( ) ( ) . Thus, we have

Y x y x h y x h y x h y x
Thus, the following theorem have been proven.
Theorem 1.Let ( ) y x be a sufficiently smooth function for ≤ ≤ a x b and further assume that ( ) Y x is the cubic B-spline approximation defined by Equation ( 17) for ( ) y x .Then, at the mesh points i x for 0,1,2, ,

Y x y x h y x h y x h y x
, Let us consider ( ) ( ) ( ) = − e x Y x y x as the error term.By substituting Equations (24)(25)(26) into the Taylor series expansion of ( ) where 0,1 Hence, it is clear that the proposed new CBS approximation is ( )

Results
In this section, the accuracy of the new CBS approximation method was measured and compared with exact solutions and existing methods.The accuracy was measured using the maximum error, defined as ( ) ( ) The numerical computations were performed using Matlab R2018a running on an Intel(R) CORE(TM) i7 CPU 1.30GHz processor, 8.00GB RAM.

Problem 1
Consider the following linear second-order differential equation [11]: which is subject to the following boundary conditions: The exact solution is ( ) ( ) Table 2 lists the absolute error when . It is clear that the proposed method agreed with the exact solutions.In Table 3, the comparison of maximum errors between the proposed method, FDM, FEM, FVM, BSI, extended cubic B-Spline interpolation using Newton's method (ECBI(N)) and extended cubic B-Spline interpolation using a built-in function (ECBI(B)) are given with different values of N to show the impact of N .Evidently, our results were better than the others.

Problem 2
Consider the following problem [3]: having the boundary conditions ( ) ( ) The true solution is ( ) .We can say that the proposed method fit well with the exact solutions.Problem 2 was studied using LSM, FDM and BSI.As can be seen in Table 5, a comparison of the maximum errors between the proposed method and the mentioned methods is given with different values of N .Obviously, our method generated more accurate results compared with the others.6.We can say that the proposed method was in good agreement with the exact solutions.The maximum errors between the proposed method and BSI method for 20, 50 N = and 100 are presented in Table 7.Clearly, our method produced better approximations compared with BSI.

Problem 4
Consider the following equation [27]: with boundary conditions ( ) The analytical solution is ( ) ( ) ( ) ( ) In Table 8, the absolute errors for 1 / 20 = h are provided.Clearly, the proposed method worked well with the exact solutions.The maximum errors between the proposed method and the BSI method for 20, 50 N = and 100 are demonstrated in Table 9.Again, our method produced better approximations compared with BSI.

Discussion
The numerical results, exact solution and absolute error for each problem at specific values of subintervals are presented in Tables 2, 4, 6 and 8. Additionally, for each problem, the details of the maximum errors of the existing methods at different values of N are given in Tables 3, 5, 7 and 9. From our derivations, the truncation error of the proposed method was ( ) 5 O h accurate.On the other hand, LSM produced truncation errors that were ( )

4
O h accurate, while FDM, FEM, FVM, BSI and ECBI were ( ) 2 O h accurate.Due to this fact, the approximation obtained by the proposed method was more accurate compared with the others.One potential extension to the present study is the generation of optimal trajectories, obtained from variational calculus to optimal control problems when the boundary conditions are given.

Conclusions
This study set out to numerically solve linear two-point BVPs.The method presented was based on a CBS basis function that used the new approximation for the second-order derivative.Four examples were considered and compared with the exact solutions to validate the accuracy of the proposed method.Investigation of the proposed method as a solution to the four problems also showed that it gave more precise and accurate results compared with FDM, FEM, FVM, ECBI(N), ECBI(B), BSI and LSM by calculating the maximum error.We found that the errors decreased, which led to higher accuracy as the step size increased.By performing the error analysis, we found that the proposed method provided fifth-order accuracy.Hence, it was concluded that the proposed method was effective to solve linear two-point BVPs.
typical CBS basis function is defined as[23]

Table 1 .
Coefficients of

Table 2 .
Absolute errors for Problem 1 at

Table 3 .
Comparison of maximum errors for Problem 1.

Table 4 .
Absolute errors for Problem 2 at

Table 5 .
Comparison of maximum errors for Problem 2.

Table 6 .
Absolute errors for Problem 3 at

Table 7 .
Comparison of maximum errors for Problem 3.

Table 8 .
Absolute errors for Problem 4 at

Table 9 .
Comparison of maximum errors for Problem 4.