Subscribe Now Subscribe Today
Research Article

Lacunary Spline Solutions of Fourth Order Initial Value Problem

K.H. Jwamer and R.G. Karem

In this study, we treat for the first time a Lacunary data (0, 1, 4) by constructing spline function of degree six in all subinterval of given partition which interpolates the Lacunary data (0, 1, 4) and the constructed spline function. We applied it to solve the fourth order initial value problem which is defined in section one. The numerical example showed that the presented spline function proved their effectiveness in solving the fourth order initial value problem. Also, we note that, the better error bounds are obtained for a small step size h.

Related Articles in ASCI
Similar Articles in this Journal
Search in Google Scholar
View Citation
Report Citation

  How to cite this article:

K.H. Jwamer and R.G. Karem, 2010. Lacunary Spline Solutions of Fourth Order Initial Value Problem. Asian Journal of Mathematics & Statistics, 3: 119-129.

DOI: 10.3923/ajms.2010.119.129



Initial value problems occur in a number of areas of applied mathematics among which are fluid mechanics, elasticity and quantum mechanics as well as science and engineering. Only small class of differential equations can be solved.

Several researchers have investigated some numerical methods for solving initial value problems, among which include cubic spline method, finite difference method, multi-derivative method and finite element method (Burden and Faires, 2001; Kanth and Vishnu, 2006).

The literature on the numerical solutions of initial value and boundary value problems by using lacunary spline functions is not too much. Saxena and Venturino (1994) used two-point boundary value problem by using lacunary spline function which interpolates the lacunary data (0, 2) and also Saki and Usmani (1983), Saxena and Venturino (1994), Siddiqi and Akram (2003, 2007), Siddiqi et al. (2007, 2008) are used qunitic spline functions for different boundary value problems. Saeed and Jwamer (2004) showed the existences, uniqueness and error bounds for lacunary interpolation by six degree spline, that spline has degree six only on the first and the final subintervals in given partitions and Jwamer (2005) showed theoretically only this spline which is found by Saeed and Jwamer (2004) approximate to the second order Cauchy problem but in the present work we showed the existence, uniqueness and error bounds for the same lacunary which is (0, 1, 4) but the degree of spline six at the all subintervals of given partition. Also, we showed theoretically and practically, that this spline approximate to the solution of fourth order initial value problem which defined in Eq. 1. Al-Bayati et al. (2009) and Sallam and Hussien (1984) are used different lacunary type of approximation solutions for different initial value problems.

The general fourth order initial value problem is considered of the form:


With the initial conditions:


By using that fεCn-1 ([0, 1]xR4) and that f is Lipschitz continuous in y, y', y'' and y'''.

The lacunary interpolation problem, which we have investigated in this study, consists in finding the six degree spline S (x) of deficiency four, interpolating data given on the function value and first and fourth order in the interval [0, 1]. Also, an extra initial condition is prescribed on the second derivative.


Here, we present six degree spline interpolation for one dimensional and given sufficiently smooth function f (x) defined on I = [0, 1] and Δn: 0 = x0<x1<x2<...<xn = 1 denote the uniform partition of I with knots xi = ih,where, I = 0, 1, 2,..., n. We denote by S4n,6 the class of six degree splines S (x) such that:


on the interval [x0, x1] where, a0,j, j = 3, 5 and 6 unknowns are to be determined.

Let us examine now intervals [xi, xi+t], i = 1, 2, … n-2. by taking into account the interpolating conditions, we can write the expression, for Si (x) in the following form:


where, ai,j, t = 1 (1) (n-1), j = 2, 3, 5 and 6 unknowns we need to be determine it.

On the last interval [xn-1, xn] we define Sn-1 (x) as follows:


where, an-1,j, j = 2, 3 5 and 6 unknowns are to be determined.


Here, the existence and uniqueness theorem for spline function of degree six which interpolate the lacunary data (0, 1, 4 ) is presented and examined.

Theorem 1: Existence and Uniqueness Spline Function
Given the real numbers y(r), i = 0, 1, 2,…, n and r = 0, 1, 4 and r = 0, 1, 4 and y'' (x0) and y'' (xn) then, there exist a unique spline of degree six as given in the Eq. 3-5 such that:


We defined the spline function S(x) as follows:

where, the coefficients of these polynomials are to be determined by the following conditions:




To find, uniquely, the coefficients in S0 (x) of Eq. 3 by using the condition 7 where, i = 0, we obtain the following:





From the boundary condition 8 we have:


Solving Eq. 9-11 to obtain the following:




Substituting these values of a0,3, a0,5 and a0,6 we get:


We shall find the coefficients of Si (x) on the interval [xi,xi+1] for I = 1, 2, 3,..., n-2.

Here we have:






From Eq. 9 we have:

substitute it in each of Eq. 17, 18 and 20 and we obtain the following linear system:




Clearly, the above linear system Eq. 22 has unique solutions and after solving it we obtain the following relations:





Note that the values ai,2, ai,3, ai,5, and ai,6 for i = 1, 2,…, n-2 has unique and to find each of them only the first step means for i = 1 use Eq. 16 and after that it is easy that we show the above linear has unique solutions.

Finally, for finding the coefficients of Sn-1 (x), we have:






Solving Eq. 27-29, we see that the coefficients an-1,i; i = 2, 3, 5 and 6 are uniquely determined. Hence the proof of Theorem 1 is completed.


The error bound of the spline function S (x) which is a solution of the problem 6 is obtained for the uniform partition I by the following theorem:

Theorem 2
Let y∈C6 [0, 1] and S (x) be a unique spline function of degree six which is a solution of the problem 6. Then for x∈ [xi, xi+1]; i = 0, 1, 2,…, n-1:

And for i = 1, 2, 3,…, n-2:

And for i = n-1:

where, W6 (h) denotes the modules of continuity of y(6), defined by W6 (h) = max {|y(6) (x)-y(6) (y)|; where |x-y|<h and x, y∈[0, 1]}

To prove this theorem we need the following lemma:

Lemma 1
Let y∈C6 [0, 1]. Then,

for i = 0, 1, …, n-1



And W6 (h) denotes the modules of continuity of y(6).

Proof of Lemma 1
If y∈C6 [0, 1] then using Taylor’s expansion formula, we have:

where, xii<xi+1 and similar expressions for the derivatives of y (x) can be used.

Now from Eq. 23 and using Eq. 31 we obtain:


where, xis,i<xi+1 for i = 1, 2,..., n; s = 1, 2, 3, 4 and:


where, xi1,0, θ2,0, θ3,0, θ4,0<x1

We see that the system of Eq. 32 and 33 is unknown ei,2, i = 1, 2, …,n-1 has the unique solution:


It is clear that:


This completes the proof of the lemma 1.

Proof of Theorem 2
Let xε [xi, xi+1] where, i = 1, 2, …, n-2.

We have from Eq. 4 with apply Taylor’s expansion formula we have


Using Eq. 26 and 34, lemma 1 and Taylor’s series, we have


From Eq. 4 we have:


From Eq. 36 we obtain:


From Eq. 25, 35 and 37 and lemma 1 we get:




From Eq. 6, we have Si(4) (xi)-y(4) (xi) = 0, from which and using Eq. 38 and 39 we obtain:


To find |Si(3) (x)-y(3) (x)|, from Eq. 4 we have:


Using Eq. 24, 35, 38 and Taylor’s formula in Eq. 41 we obtain


To find |Si(2) (x)-y(2) (x)|, from Eq. 4 we have:


Using Eq. 25, 35, 38 and 42, lemma 1 and Taylor’s formula in Eq. 43 we obtain




From Eq. 6, we have Si' (xi)-y' (xi) = 0, then using Eq. 44 and 45 we obtain:




And since from Eq. 6, we have Si (xi)-y (xi) = 0, then put it in above Eq. 47 and using Eq. 46 we obtain:

This proves Theorem 2 for x∈ [xi, xi+1], I = 1, 2,..., n-1.

For x∈ [x0, x1], we have from Eq. 3:


Using Eq. 15 and Taylor’s series in Eq. 48 we obtain:

Carrying on similar steps as for the case x∈ [xi, xi+1], i = 1, 2, …, n-1, we find the following inequalities:

This proves Theorem 2 for x∈ [x0, x1]. And for x∈ [xn-1, xn] the same manner in above is used to get an-1,2, an-1,3, an-1,5 and an-1,6, the proof of Theorem 2 is complete.


This study presents numerical result to demonstrate the convergence of the spline function of degree six which is constructed in the section three to the fourth order initial value problem.

We consider that the fourth order initial value problem y(4)-y = 0 where, x∈ [0, 1], y (0) = y'(0) = y''(0) = y'''(0) = 1.

From Eq. 3 it’s easy to verify that:

Table 1: Absolute maximum error for S (x) and it’s derivative

Also it is easy from Eq. 3 and 4 to verify that:

From Eq. 6 we have:

From Eq. 3 and 4, with using the values of ai,j, i = 0, 1,…,n-1 and j = 2, 3, 5, 6 given in the Eq. 13-16 and 23-26, we get:

It turns out that the six degree spline which is presented in this study, yield approximate solution that is O(h6) as stated in Theorem 1. The results are shown in the Table 1 for different step sizes h in the algorithm similar to Al-Bayati et al. (2009).


In this study, we get to the conclusion that, this spline which is defined in Eq. 3-5 approximate to the solution of Eq. 1 and also if we change the step size h by taking small value we obtain the best approximate solution to the exact solution of Eq. 1.

Al-Bayati, A.Y., R.K. Saeed and F.K. Hama-Salh, 2009. The existence, uniquenes and errors bounds for approximation splines interpolation for solving second-order initial value problems. Am. J. Math. Statistics, 5: 123-129.
Direct Link  |  

Burden, R.L. and J.D. Faires, 2001. Numerical Analysis. 7th Edn., Brooks Cole Pacific, Grove, ISBN: 0534382169.

Jwamer, K.H., 2005. Solution of cauchys problem by using spline interpolations. J. Al-Nahrain Univ. Iraq, 8: 97-99.

Kanth, A.S V. and B. Vishnu, 2006. Cubic spline for a class on non-linear singular boundary value problems arising in physiology. Applied Math. Comput., 174: 768-774.
Direct Link  |  

Saeed, R.K. and K.H. Jwamer, 2004. Lacunary interpolation by splines function (0, 1, 4) case. J. Dohuk Univ., 4: 193-196.

Saki, M. and R.Usmani, 1983. Spline solutions for nonlinear fourth order two-point boundary value problem. Publ. Res. Math. Sci., 19: 135-144.

Sallam, S. and M.A. Hussien, 1984. Deficient spline function approximation to second-order differential equations. Applied Math. Modell., 6: 408-412.
CrossRef  |  

Saxena, A. and E. Venturino, 1994. Solving two-point boundary value problems by means of deficient quartic splines. Applied Math. Comput., 66: 25-40.
CrossRef  |  Direct Link  |  

Siddiqi, S.S. and G. Akram, 2003. Quintic spline solutions of forth order boundary value problem. ArXiv; Math. Nal 0306357, pp: 12.

Siddiqi, S.S. and G. Akram, 2007. Solutions of tenth-order boundary value problems using eleven degree splines. Applied Math. Comput., 185: 115-127.
Direct Link  |  

Siddiqi, S.S., G. Akram and A. Elahi, 2008. Quintic spline solution of linear fifth-Order boundary value problems. Applied Math. 196: 214-220.

Siddiqi, S.S., G. Akram and S. Nazeem, 2007. Quintic spline solution of linear sixth-order boundary value problems. Applied Math. Comput., 189: 887-892.
CrossRef  |  

©  2020 Science Alert. All Rights Reserved