JAMPJournal of Applied Mathematics and Physics2327-4352Scientific Research Publishing10.4236/jamp.2016.48149JAMP-69599ArticlesPhysics&Mathematics Septic B-Spline Solution of Fifth-Order Boundary Value Problems BinLin1School of Mathematics and Computation Science, Lingnan Normal University, Zhanjiang, China* E-mail:040820160408144614548 July 2016accepted 6 August 9 August 2016© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

A numerical method based on septic B-spline function is presented for the solution of linear and nonlinear fifth-order boundary value problems. The method is fourth order convergent. We use the quesilinearization technique to reduce the nonlinear problems to linear problems and use B-spline collocation method, which leads to a seven nonzero bands linear system. Illustrative example is included to demonstrate the validity and applicability of the proposed techniques.

Septic B-Spline Function Fifth-Order Boundary Value Problems B-Spline Collocation Method Nonlinear Problems
1. Introduction

Consider the following fifth-order boundary value problem.

With boundary conditions

where are known real constants, and are continuous on. This problem arising in the mathematical modeling of viscoelastic flows   has been studied by several authors  -  . A. Lamnii, H. Mraoui, D. Sbibih and A. Tijini studied the fifth-order boundary value problem based on splines quasi-interpolants and proved to be second order convergent.

B-spline functions based on piece polynomials are useful wavelet basis functions, the resulting matrices are sparse, but always, banded. And that possess attractive properties: piecewise smooth, compact support, symmetry, rapidly decaying, differentiability, linear combination, B-splines were introduced by Schoenberg in 1946  . Up to now, B-spline approximation method for numerical solutions has been researched by various researchers  -  .

In this paper, the septic B-spline function is used as a basis function and the B-spline collocation method is studied to solve the linear and nonlinear fifth-order boundary value problems. The method is fourth order convergent. We use the quesilinearization technique to reduce the nonlinear problems to linear problems. The present method is tested for its efficiency by considering two examples.

2. Septic B-Spline Interpolation

An arbitrary Nth order spline function with compact support of N. It is a concatenation of N sections of (N-1)th order polynomials, continuous at the junctions or “knots”, and gives continuous (N-1)th derivatives at the junctions.

Let be a uniform partition of such that, , where. Let the septic B-spline function with knots at the points be given by

The set of splines forms a basis for the functions defined over. The values of and its derivatives are as shown in Table 1.

We seek the approximation to the exact solution, which uses these septic B-splines:

which satisfies the following interpolation conditions:

where are unknown real coefficients.

Using the septic B-spline function Equation (3) and the approximate solution Equation (4), the nodal values and at the node are given in terms of element parameters by

The values of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1720638x33.png" xlink:type="simple"/></inline-formula> and its derivatives with knots
0112011912416119112010
000
00
000
0000
000

From Equations (4)-(7), we have

Using operator notations, we obtain

Expanding them in powers of, we obtain

Hence we get

3. Spline Collocation Method3.1. Linear Problems

From Equation (1) and Equation (12), we can get

Using the boundary conditions and by neglecting the error of Equation (13), we can obtain following linear equations

Or

where

where

T denoting transpose.

In which B is a square matrix of order N + 7 with seven nonzero bands. Since B is nonsingular, after solving the linear system Equation (15) for, we can obtain the septic spline approximate

solution with the accuracy being.

3.2. Nonlinear Problems

Consider the nonlinear fifth order boundary value problem

with boundary conditions

We use the quesilinearization technique to reduce the above nonlinear problem to a sequence of linear problems. Expanding the right hand side of Equation (16), we have

Equation (18) can be rewritten as

where

Equation (19) once the initial values (k = 0, , ,) has been computed from the initial conditions, Equation (19) becomes into a linear equations with constant coefficients. Equation (19) can be solved by using iterative method.

Subject to the boundary conditions

Instead of solving nonlinear problem (16) with boundary conditions (17), we solve a sequence of linear problems (19) with boundary conditions (20), we consider as the numerical solution to nonlinear problem (16) with boundary conditions (17).

4. Computation of Error

The relative error of numerical solution is given by

The pointwise errors are given by

The maximum pointwise errors are given by

5. Numerical Tests

In the section, we illustrate the numerical techniques discussed in the previous section by the following problems.

Example 1. Consider the following equation  -  :

With boundary conditions

The exact solution is given by

The numerical results are shown in Table 2, the comparison of maximum absolute errors are given by Table 3. The relative errors for different values of h are seen in Figure 1. The pointwise errors of example are given in Figure 2. The maximum pointwise errors for different values of h are given in Figure 3.

Example 2. Consider the following nonlinear equation    .

Maximum absolute errors, relative error for example 1
hCPU time (seconds)
1/81.718286042191597e−0043.294552154146086e−0043.218
1/106.447839923018339e−0051.351847113984665e−0048.172
1/161.028885985859818e−0052.068672045465237e−0052.953
1/204.192815893866442e−0068.479871033239017e−0065.391
1/326.368030709968942e−0071.296221338426021e−0066.922
1/402.475231232201836e−0075.045908672756208e−0078.688
1/501.671824484961171e−0073.420027458407504e−0078.172
1/643.108106372273767e−0086.343766826003454e−0086.921
Comparison of maximum absolute errors for example 1
h
Our methodCaglar et al. Shahid.et al. Khan et al. 
1/106.447839923018339E−50.15702.259E−44.025E−3
1/204.192815893866442E−60.07471.33E−53.911E−3
1/402.475231232201836E−70.02085.2812E−71.145E−2
The relative errors of example 1 for different values of h The pointwise errors of example 1 The maximum pointwise errors of example 1 for different values of h

With boundary conditions

The exact solution is given by.

Comparison of numerical results and pointwise errors are given in Table 4. The numerical result is found in good agreement with exact solution.

Example 2. Comparison of results and pointwise errors
XNumericalExactOur errorsErrors of Errors of Errors of Errors of 
0.11.105170891343271.105170918075652.673237986527965e−0087.0e−41.3e−72.3e−70
0.21.221402661375461.221402758160179.678471002416700e−0087.2e−44.2e−71.6e−61.0e−5
0.31.349858642596011.349858807576001.649799901137783e−0074.1e−47.2e−74.6e−61.0e−5
0.41.491824482622821.491824697641272.150184499338792e−0074.6e−49.4e−78.9e−61.0e−4
0.51.648721270700131.648721034678922.360212099095094e−0074.7e−41.0e−61.3e−53.2e−4
0.61.822118582738881.822118800390512.176516300522735e−0074.8e−49.3e−71.6e−53.6e−4
0.72.013752540823812.013752707470481.666466702410219e−0073.9e−47.1e−71.6e−61.4e−4
0.82.225540831634892.225540928492479.685758017852209e−0083.1e−44.1e−71.2e−53.1e−4
0.92.459603080489082.459603111156953.066787002126148e−0081.6e−41.3e−75.1e−65.8e−4
6. Conclusion

In the paper, the fifth-order boundary value problems are solved by means of septic B-splines collocation method. We use the quesilinearization technique to reduce the nonlinear problems to linear problems and reduce a boundary value problem to the solution of algebraic equations with seven nonzero bands. The numerical results show that the present method is relatively simple to collocate the solution at the mesh points and easily carried out by a computer and approximates the exact solution very well.

Acknowledgements

The authors would like to thank the editor and the reviewers for their valuable comments and suggestions to improve the results of this paper. This work was supported by the Natural Science Foundation of Guangdong (2015A030313827).

Cite this paper

Bin Lin, (2016) Septic B-Spline Solution of Fifth-Order Boundary Value Problems. Journal of Applied Mathematics and Physics,04,1446-1454. doi: 10.4236/jamp.2016.48149

ReferencesKarageorghis, A., Phillips, T.N. and Davies, A.R. (1988) Spectral Collocation Methods for the Primary Two-Point Boundary Value Problem in Modeling Viscoelastic Flows. International Journal for Numerical Methods in Engineering, 26, 805-813. http://dx.doi.org/10.1002/nme.1620260404Davies, A.R., Karageorghis, A. and Phillips, T.N. (1988) Spectral Galerkin Methods for the Primary Two-Point Boundary Value Problem in Modeling Viscoelastic Flows. International Journal for Numerical Methods in Engineering, 26, 647-662. http://dx.doi.org/10.1002/nme.1620260309Caglar, H.N., Caglar, S.H. and Twizell, E.H. (1999) The Numerical Solution of Fifth-Order Boundary Value Problems with Sixth-Degree B-Spline Functions. Applied Mathematics Letters, 12, 25-30. http://dx.doi.org/10.1016/S0893-9659(99)00052-XSiddiqi, S.S. and Akram, G. (2006) Solutions of Fifth Order Boundary Value Problems Using Nonpolynomial Spline Technique. Applied Mathematics and Computation, 175, 1574-1581. http://dx.doi.org/10.1016/j.amc.2005.09.004Lamnii, A., Mraoui, H., Sbibih, D. and Tijini, A. (2008) Sextic Spline Solution of Fifth Order Boundary Value Problems. Mathematics and Computers in Simulation, 77, 237-246. http://dx.doi.org/10.1016/j.matcom.2007.09.008De Boor, C. (1978) A Practical Guide to Splines. Springer-Verlag, pp. 54, 105. http://dx.doi.org/10.1007/978-1-4612-6333-3Saka, B. and Dag, I. (2007) Quartic B-Spline Collocation Method to the Numerical Solutions of the Burgers’ Equation. Chaos, Solitons & Fractals, 32, 1125-1137. http://dx.doi.org/10.1016/j.chaos.2005.11.037Ramadan, M.A., EI-Danaf, T.S. and Alaal, F. (2005) A Numerical Solution of the Burgers’ Equation Using Septic B-Splines. Chaos, Solitons & Fractals, 26, 795-804. http://dx.doi.org/10.1016/j.chaos.2005.01.054Ramadan, M.A., EI-Danaf, T.S. and Abd Alaal, F.E.I. (2005) A Numerical Solution of the Burgers’ Equation Using Septic B-Splines. Chaos, Solitons & Fractals, 26, 1249-1258. http://dx.doi.org/10.1016/j.chaos.2005.02.019&Ccedilaglar, H., &Ccedilaglar, N. and &Oumlzer, M. (2009) B-Spline Solution of Non-Linear Singular Boundary Value Problems Arising in Physiology. Chaos, Solitons & Fractals, 39, 1232-1237. http://dx.doi.org/10.1016/j.chaos.2007.06.007&Ccedilaglar, H., &Oumlzer, M. and &Ccedilaglar, N. (2008) The Numerical Solution of the One-Dimensional Heat Equation by Using Third Degree B-Spline Functions. Chaos, Solitons & Fractals, 38, 1197-1201. http://dx.doi.org/10.1016/j.chaos.2007.01.056Saka, B. and Dag, I. (2007) Quartic B-Spline Collocation Method to the Numerical Solutions of the Burgers’ Equation. Chaos, Solitons & Fractals, 32, 1125-1137. http://dx.doi.org/10.1016/j.chaos.2005.11.037Jain, P.C., Shankar, R. and Bhardwaj, D. (1997) Numerical Solution of the Korteweg-de Vries (KdV) Equation. Chaos, Solitons & Fractals, 8, 943-951. http://dx.doi.org/10.1016/S0960-0779(96)00135-XLin, B., Li, K.T. and Cheng, Z.X. (2009) B-Spline Solution of a Singularly Perturbed Boundary Value Problem Arising in Biology. Chaos, Solitons & Fractals, 42, 2934-2948. http://dx.doi.org/10.1016/j.chaos.2009.04.036Caglar, H.N., Caglar, S.H. and Twizell, E.N. (1999) The Numerical Solution of Fifth-Order Boundary Value Problems with Sixth-Degree B-Spline Functions. Applied Mathematics Letters, 12, 25-30. http://dx.doi.org/10.1016/S0893-9659(99)00052-XSiddiqi, S.S. and Akram, G. (2007) Sextic Spline Solutions of Fifth-Order Boundary Value Problems. Applied Mathematics Letters, 20, 591-597. http://dx.doi.org/10.1016/j.aml.2006.06.012Khan, M. (1994) Finite Difference Solutions of Fifth-Order Boundary Value Problems. Ph.D. Thesis, Brunel University, England.Noor, M.A. and Mohyud-Din, S.T. (2009) A New Approach to Fifth-Order Boundary Value Problems. International Journal of Nonlinear Science, 7, 143-148.Zhang, J. (2009) The Numerical Solution of Fifth-Order Boundary Value Problems by the Variational Iteration Method. Computers & Mathematics with Applications, 58, 2347-2350. http://dx.doi.org/10.1016/j.camwa.2009.03.073