A Novel Method for Solutions of Fourth-Order Fractional Boundary Value Problems

Abstract

In this paper, we find the solutions of fourth order fractional boundary value problems by using the reproducing kernel Hilbert space method. Firstly, the reproducing kernel Hilbert space method is introduced and then the method is applied to this kind problems. The experiments are discussed and the approximate solutions are obtained to be more correct compared to the other obtained results in the literature.

1. Introduction

Boundary value problems come into view in many areas of science and engineering [1]. Many numerical methods have been presented for solving boundary value problems recently [2]. Wazwaz [3] has presented a modified decomposition method for investigating a special class of fourth order boundary value problems. The series solution of fourth order boundary value problems has been investigated by Adomian decomposition method in [4]. The differential transform method was implemented in [5] for series solution of fourth order boundary value problems. The approximate solution of two-point fourth order boundary value problems utilizing non-polynomial quintic spline functions have been given in [6,7]. Lodhi and Mishra [8] have enhanced a numerical method for solving fourth order boundary value problems utilizing quintic B-spline functions. A numerical technique depending on non-polynomial spline functions has been applied in [9] for the numerical solution of self-adjoint singularly perturbed fourth order boundary value problems. Akram and Amin [10] have applied quintic spline collocation method for approximate solution of singularly perturbed fourth order boundary value problems.

Fractional order boundary value problems have taken much interest by many investigators due to their implementations in many areas of science and engineering. Recently, much attentions have been given to investigate fractional order differential equations [11].

We take into consideration the following fourth order fractional boundary value problems arising in plate deflection theory by the reproducing kernel Hilbert space method.

$$u^{\left(4\right)}\left(z\right)+D^{\gamma }p\left(z\right)u\left(z\right)=g\left(z\right),z\in [a,b],$$

(1)

subject to the end conditions

$$u\left(a\right)-a_1=u\left(b\right)-a_2=u^{″}\left(a\right)-b_1=u^{″}\left(b\right)-b_2=0.$$

(2)

where $a_i^{\prime }s,b_i^{\prime }s,(i=1,2)$ are real constants and $0\le \gamma \le 1.$$p\left(z\right)$ and $g\left(z\right)$ are continuous on the $[a,b].$

Let H be a Hilbert space of functions on some set X, which every point consideration at points in X are continuous in the norm of H. Such space H, can be given by a positive definite kernel. This fact is first presented by Aronszajn [12]. These Hilbert spaces are called the reproducing kernel Hilbert spaces [13].

Arqub [14] has found the approximate solutions of DASs with nonclassical boundary conditions using a novel reproducing kernel algorithm. Arqub [15] has investigated the reproducing kernel algorithm for handling differential algebraic systems of ordinary differential equations. Arqub [16] has applied a fitted reproducing kernel Hilbert space method for the solutions of some certain classes of time-fractional partial differential equations subject to initial and Neumann boundary conditions. Azarnavid et al. [17] have investigated the Picard-reproducing kernel Hilbert space method for solving generalized singular nonlinear Lane–Emden type equations. Azarnavid et al. [18] have found multiplicity results by shooting reproducing kernel Hilbert space method for the catalytic reaction in a flat particle. Azarnavid et al. [19] have used an iterative kernel based method for fourth order nonlinear equation with nonlinear boundary condition.

We arrange the paper as: The main notions of fractional calculus are given in Section 2. The reproducing kernel Hilbert space method for fourth-order fractional order differential equations is presented in Section 3. The numerical results and discussions are presented in Section 4. The conclusion is presented in the last section.

2. Preliminaries and Notations

We give some definitions of fractional calculus in this section. There are many definitions of fractional derivatives. The definitions of Riemann–Liouville fractional integral and Caputo’s fractional derivative are presented as [17]:

Definition 1.

The Riemann–Liouville left and right fractional integral of order $\gamma >0$ is given as [11]:

$$I_{a^{+}}^{\gamma }u\left(z\right)=\frac{1}{\Gamma \left(\gamma \right)}{\int }_a^z{(z-\tau )}^{\gamma -1}u\left(\tau \right)\mathrm{d}\tau ,m-1<\gamma \le m,m\in N,$$

(3)

and

$$I_{b^{-}}^{\gamma }u\left(z\right)=-\frac{1}{\Gamma \left(\gamma \right)}{\int }_b^z{(z-\tau )}^{\gamma -1}u\left(\tau \right)\mathrm{d}\tau ,m-1<\gamma \le m,m\in N.$$

(4)

where Γ defines the Gamma function.

Definition 2.

The Riemann–Liouville fractional derivative of order $\gamma >0$ is presented as [11]:

$$D_a^{\gamma }u\left(z\right)=\frac{1}{\Gamma (m-\gamma )}\frac{d^m}{d{z}^m}{\int }_a^z{(z-\tau )}^{m-\gamma -1}u\left(\tau \right)\mathrm{d}\tau ,m-1<\gamma \le m,m\in N.$$

(5)

Definition 3.

The Caputo’s fractional derivative of order $\gamma >0$ is presented as [11]:

$$D_a^{\gamma }u\left(z\right)=\frac{1}{\Gamma (m-\gamma )}{\int }_a^z{(z-\tau )}^{m-\alpha -1}\frac{d^{m}u\left(\tau \right)}{d{\tau }^m}\mathrm{d}\tau ,m-1<\alpha \le m,m\in N.$$

(6)

Lemma 1.

Assume that f and g are analytic functions on $(a-h,a+h)$ and let $0<\gamma \le 1.$ Then we have [11]:

$$\begin{array}{ccc}\hfill D_a^{\gamma }\left[f\left(z\right)g\left(z\right)\right]& =& \frac{{(z-a)}^{-\gamma }}{\Gamma (1-\gamma )}g\left(a\right)(f\left(z\right)-f\left(a\right))+\left(D_a^{\alpha }g\left(z\right)\right)f\left(z\right)\hfill \\ & +& \sum _{l=1}^{\infty }\left(\genfrac{}{}{0pt}{}{\alpha }{l}\right)\left(I_a^{1-\alpha }g\left(z\right)\right)D_a^{l}f\left(z\right).\hfill \end{array}$$

(7)

3. Reproducing Kernel Hilbert Space Method

We will construct reproducing kernel Hilbert spaces to solve the problem.

Definition 4.

$V_2^1[0,1]$ is the first reproducing kernel Hilbert space that we need.

$$V_2^1[0,1]=\{r\in AC[0,1]:r^{\prime }\in L^2[0,1]\}.$$

The inner product and norm for this reproducing kernel Hilbert space are defined by:

$${〈r,q〉}_{V_2^1}=r\left(0\right)q\left(0\right)+{\int }_0^1{r}^{\prime }\left(z\right)q^{\prime }\left(z\right)\mathrm{d}z,r,q\in V_2^1[0,1]$$

and

$${∥r∥}_{V_2^1}=\sqrt{{〈r,r〉}_{V_2^1}},r\in V_2^1[0,1].$$

Lemma 2.

The reproducing kernel function $U_z$ of $V_2^1[0,1]$ is obtained as:

$$U_t\left(z\right)=\left\{\begin{array}{c}1+z,z\le t,\hfill \\ 1+t,z>t.\hfill \end{array}\right.$$

Definition 5.

We construct the reproducing kernel Hilbert space $V_2^5[0,1]$ as:

$$\begin{array}{c}\hfill V_2^5[0,1]=\{r\in AC[0,1]:r^{\prime },r^{″},r^{‴},r^{\left(4\right)}\in AC[0,1],r^{\left(5\right)}\in L^2[0,1],\\ \hfill r\left(0\right)=r\left(1\right)=r^{″}\left(0\right)=r^{″}\left(1\right)=0\}.\end{array}$$

The inner product and the norm for this special Hilbert space is defined as:

$${〈r,q〉}_{V_2^5}=\sum _{i=0}^4{r}^{\left(i\right)}\left(0\right)q^{\left(i\right)}\left(0\right)+{\int }_0^1{r}^{\left(5\right)}\left(z\right)q^{\left(5\right)}\left(z\right)\mathrm{d}z,r,q\in V_2^5[0,1]$$

and

$${∥r∥}_{V_2^5}=\sqrt{{〈r,r〉}_{V_2^5}},r\in V_2^5[0,1].$$

Theorem 1.

We obtain the reproducing kernel function for the reproducing kernel Hilbert space $V_2^5[0,1]$ by:

$$\begin{array}{c}\hfill A_y\left(x\right)=\left\{\begin{array}{cc}{a}_y\left(x\right),\hfill & x\le y,\hfill \\ b_y\left(x\right),\hfill & x>y.\hfill \end{array}\right.\end{array}$$

where,

$$\begin{array}{ccc}\hfill a_y\left(x\right)& =\hfill & \frac{x^9}{362880}+\frac{635131x^6{y}^4}{28582381440}+\frac{11x^6{y}^9}{142911907200}-\frac{11x^6{y}^8}{15879100800}+\frac{x^7{y}^2}{10080}\hfill \\ \hfill & \hfill & +\frac{231x^{6}y}{8270365}-\frac{105619x^7{y}^7}{3001150051200}+\frac{91723x^7{y}^6}{428735721600}+\frac{23x^5{y}^9}{31758201600}\hfill \\ \hfill & \hfill & -\frac{88249x^7{y}^5}{1442911907200}-\frac{88249x^7{y}^4}{28582381440}-\frac{91723x^7{y}^3}{3572797680}-\frac{193x^7{y}^9}{1000383350400}\hfill \\ \hfill & \hfill & +\frac{193x^7{y}^8}{111153705600}-\frac{579x^{7}y}{8270365}-\frac{2969x^{8}y}{66692223360}+\frac{13x^8{y}^9}{190549209600}\hfill \\ \hfill & \hfill & -\frac{13x^8{y}^8}{21172134400}+\frac{193x^8{y}^7}{111153705600}-\frac{11x^8{y}^6}{15879100800}-\frac{69x^8{y}^5}{10586067200}\hfill \\ \hfill & \hfill & -\frac{69x^8{y}^4}{2117213440}+\frac{11x^8{y}^3}{132325840}-\frac{13x^9{y}^9}{1714942886400}+\frac{13x^9{y}^8}{190549209600}\hfill \\ \hfill & \hfill & -\frac{193x^9{y}^7}{1000383350400}+\frac{11x^9{y}^6}{142911907200}+\frac{23x^9{y}^5}{31758201600}-\frac{635131x^3{y}^4}{238186512}\hfill \\ \hfill & \hfill & +\frac{23x^9{y}^4}{6351640320}-\frac{11x^9{y}^3}{1190932560}-\frac{91x^{9}y}{33081460}+\frac{2969xy}{1654073}-\frac{127013x^6{y}^3}{2381865120}\hfill \\ \hfill & \hfill & -\frac{91x{y}^9}{33081460}+\frac{819x{y}^8}{33081460}-\frac{579x{y}^7}{8270365}+\frac{231x{y}^6}{8270365}+\frac{635131x^6{y}^5}{142911907200}\hfill \\ \hfill & \hfill & +\frac{4347x{y}^5}{16540730}+\frac{4347x{y}^4}{3308146}-\frac{5544x{y}^3}{1654073}+\frac{127013x^3{y}^3}{19848876}-\frac{636517x^6{y}^6}{428735721600}\hfill \\ \hfill & \hfill & -\frac{91723x^3{y}^7}{3572797680}+\frac{636517x^3{y}^6}{3572797680}-\frac{635131x^3{y}^5}{1190932560}-\frac{635131x^5{y}^3}{1190932560}\hfill \\ \hfill & \hfill & -\frac{11x^3{y}^9}{1190932560}+\frac{11x^3{y}^8}{132325840}-\frac{5544x^{3}y}{1654073}-\frac{222389x^4{y}^4}{158791008}+\frac{91723x^6{y}^7}{428735721600}\hfill \\ \hfill & \hfill & -\frac{88249x^4{y}^7}{28582381440}+\frac{635131x^4{y}^6}{28582381440}-\frac{319739x^4{y}^5}{4763730240}-\frac{635131x^4{y}^3}{238186512}\hfill \\ \hfill & \hfill & +\frac{23x^4{y}^9}{6351640320}-\frac{69x^4{y}^8}{2117213440}+\frac{4347x^{4}y}{3308146}+\frac{222389x^5{y}^4}{793955040}+\frac{4347x^{5}y}{16540730}\hfill \\ \hfill & \hfill & -\frac{88249x^5{y}^7}{142911907200}+\frac{635131x^5{y}^6}{142911907200}-\frac{319739x^5{y}^5}{23818651200}-\frac{69x^5{y}^8}{10586067200}.\hfill \end{array}$$

The reproducing kernel function $A_y\left(x\right)$ is symmetric. Therefore, when exchange x and y in $a_y\left(x\right),$ we will obtain $b_y\left(x\right).$

Proof. 

We have

$${〈r,A_y〉}_{V_2^5}=\sum _{i=0}^4{A}_y^{\left(i\right)}\left(0\right)r^{\left(i\right)}\left(0\right)+{\int }_0^1{A}_y^{\left(5\right)}\left(z\right)r^{\left(5\right)}\left(z\right)\mathrm{d}z,$$

by Definition 5. We obtain

$$\begin{array}{ccc}\hfill {〈r,A_y〉}_{V_2^5}& =\hfill & A_y\left(0\right)r\left(0\right)+A_y^{\prime }\left(0\right)r^{\prime }\left(0\right)+A_y^{″}\left(0\right)r^{″}\left(0\right)\hfill \\ \hfill & \hfill & +A_y^{\left(3\right)}\left(0\right)r^{‴}\left(0\right)+A_y^{\left(4\right)}\left(0\right)r^{\left(4\right)}\left(0\right)\hfill \\ \hfill & \hfill & +A_y^{\left(5\right)}\left(1\right)r^{\left(4\right)}\left(1\right)-A_y^{\left(5\right)}\left(0\right)r^{\left(4\right)}\left(0\right)\hfill \\ \hfill & \hfill & -A_y^{\left(6\right)}\left(1\right)r^{\left(3\right)}\left(1\right)+A_y^{\left(6\right)}\left(0\right)r^{\left(3\right)}\left(0\right)\hfill \\ \hfill & \hfill & +A_y^{\left(7\right)}\left(1\right)r^{″}\left(1\right)-A_y^{\left(7\right)}\left(0\right)r^{″}\left(0\right)\hfill \\ \hfill & \hfill & -A_y^{\left(8\right)}\left(1\right)r^{\prime }\left(1\right)+A_y^{\left(8\right)}\left(0\right)r^{\prime }\left(0\right)+{\int }_0^1{A}_y^{\left(9\right)}\left(z\right)r^{\prime }\left(z\right)\mathrm{d}z,\hfill \end{array}$$

by integration by parts. Since $r\left(0\right)=r\left(1\right)=r^{″}\left(0\right)=r^{″}\left(1\right)=0,$ we get

$$\begin{array}{ccc}\hfill {〈r,A_y〉}_{V_2^5}& =\hfill & A_y^{\prime }\left(0\right)r^{\prime }\left(0\right)+A_y^{\left(3\right)}\left(0\right)r^{\left(3\right)}\left(0\right)+A_y^{\left(4\right)}\left(0\right)r^{\left(4\right)}\left(0\right)\hfill \\ \hfill & \hfill & +A_y^{\left(5\right)}\left(1\right)r^{\left(4\right)}\left(1\right)-A_y^{\left(5\right)}\left(0\right)r^{\left(4\right)}\left(0\right)\hfill \\ \hfill & \hfill & -A_y^{\left(6\right)}\left(1\right)r^{\left(3\right)}\left(1\right)+A_y^{\left(6\right)}\left(0\right)r^{\left(3\right)}\left(0\right)\hfill \\ \hfill & \hfill & -A_y^{\left(8\right)}\left(1\right)r^{\prime }\left(1\right)+A_y^{\left(8\right)}\left(0\right)r^{\prime }\left(0\right)+{\int }_0^1{A}_y^{\left(9\right)}\left(z\right)r^{\prime }\left(z\right)\mathrm{d}z.\hfill \end{array}$$

Since

$$A_y^{\left(5\right)}\left(1\right)=A_y^{\left(6\right)}\left(1\right)=A_y^{\left(8\right)}\left(1\right)=0,$$

we get

$$\begin{array}{ccc}\hfill {〈r,A_y〉}_{V_2^5}& =\hfill & A_y^{\prime }\left(0\right)r^{\prime }\left(0\right)+A_y^{\left(3\right)}\left(0\right)r^{\left(3\right)}\left(0\right)+A_y^{\left(4\right)}\left(0\right)r^{\left(4\right)}\left(0\right)\hfill \\ \hfill & \hfill & -A_y^{\left(5\right)}\left(0\right)r^{\left(4\right)}\left(0\right)+A_y^{\left(6\right)}\left(0\right)r^{\left(3\right)}\left(0\right)\hfill \\ \hfill & \hfill & +A_y^{\left(8\right)}\left(0\right)r^{\prime }\left(0\right)+{\int }_0^1{A}_y^{\left(9\right)}\left(z\right)r^{\prime }\left(z\right)\mathrm{d}z.\hfill \end{array}$$

We have

$$\begin{array}{ccc}\hfill A_y\left(0\right)& =& \frac{2969y}{1654073}+\frac{4347y^4}{3308146}-\frac{91y^9}{33081460}+\frac{819y^8}{33081460}\hfill \\ \hfill & +& \frac{4347y^5}{16540730}-\frac{5544y^3}{1654073}-\frac{579y^7}{8270365}+\frac{231y^6}{8270365},\hfill \end{array}$$

$$\begin{array}{ccc}\hfill A_y^{\left(3\right)}\left(0\right)& =& -\frac{33264y}{1654073}-\frac{635131y^4}{39697752}-\frac{11y^9}{198488760}+\frac{33y^8}{66162920}\hfill \\ \hfill & -& \frac{635131y^5}{198488760}+\frac{127013y^3}{3308146}-\frac{91723y^7}{595466280}+\frac{636517y^6}{595466280},\hfill \end{array}$$

$$\begin{array}{ccc}\hfill A_y^{\left(4\right)}\left(0\right)& =& \frac{52164y}{1654073}+\frac{222389y^4}{6616292}+\frac{23y^9}{264651680}-\frac{207y^8}{264651680}\hfill \\ \hfill & -& \frac{319739y^5}{198488760}-\frac{635131y^3}{9924438}-\frac{88249y^7}{1190932560}+\frac{635131y^6}{1190932560},\hfill \end{array}$$

$$\begin{array}{ccc}\hfill A_y^{\left(5\right)}\left(0\right)& =& \frac{52164y}{1654073}+\frac{222389y^4}{6616292}+\frac{23y^9}{264651680}-\frac{207y^8}{264651680}\hfill \\ \hfill & -& \frac{319739y^5}{198488760}-\frac{635131y^3}{9924438}-\frac{88249y^7}{1190932560}+\frac{635131y^6}{1190932560},\hfill \end{array}$$

$$\begin{array}{ccc}\hfill A_y^{\left(6\right)}\left(0\right)& =& \frac{33264y}{1654073}+\frac{635131y^4}{39697752}+\frac{11y^9}{198488760}-\frac{33y^8}{66162920}\hfill \\ \hfill & +& \frac{635131y^5}{198488760}-\frac{127013y^3}{3308146}+\frac{91723y^7}{595466280}-\frac{636517y^6}{595466280},\hfill \end{array}$$

$$\begin{array}{ccc}\hfill A_y^{\left(8\right)}\left(0\right)& =& -\frac{2969y}{1654073}-\frac{4347y^4}{3308146}+\frac{91y^9}{33081460}-\frac{819y^8}{33081460}\hfill \\ \hfill & -& \frac{4347y^5}{16540730}+\frac{5544y^3}{1654073}+\frac{579y^7}{8270365}-\frac{231y^6}{8270365}.\hfill \end{array}$$

Therefore, we reach

$$\begin{array}{c}\hfill {〈r,A_y〉}_{V_2^5}={\int }_0^y{A}_y^{\left(9\right)}\left(z\right)r^{\prime }\left(z\right)\mathrm{d}z+{\int }_y^1{A}_y^{\left(9\right)}\left(z\right)r^{\prime }\left(z\right)\mathrm{d}z.\end{array}$$

We have

$$\begin{array}{c}\hfill A_y^{\left(9\right)}\left(z\right)=\left\{\begin{array}{cc}1+k\left(y\right),\hfill & z<y,\hfill \\ k\left(y\right),\hfill & z>y.\hfill \end{array}\right.\end{array}$$

where

$$\begin{array}{ccc}k\left(y\right)\hfill & =\hfill & \frac{819y^8}{33081460}-\frac{1651104y}{1654073}+\frac{4347y^5}{165540730}+\frac{4347y^4}{3308146}\hfill \\ \hfill & \hfill & -\frac{91y^9}{33081460}-\frac{5544y^3}{1654073}-\frac{579y^7}{8270365}+\frac{231y^6}{8270365}.\hfill \end{array}$$

Therefore, we find

$$\begin{array}{c}\hfill {〈r,A_y〉}_{V_2^5}={\int }_0^y(1+k\left(y\right))r^{\prime }\left(z\right)\mathrm{d}z+{\int }_y^1\left(k\left(y\right)\right)r^{\prime }\left(z\right)\mathrm{d}z,\end{array}$$

$$\begin{array}{c}\hfill {〈r,A_y〉}_{V_2^5}=(1+k\left(y\right))(r\left(y\right)-r\left(0\right))+k\left(y\right)(r\left(1\right)-r\left(y\right)),\end{array}$$

$$\begin{array}{c}\hfill {〈r,A_y〉}_{V_2^5}=r\left(y\right).\end{array}$$

Thus, the proof is completed. □

We investigate the solutions of the problem (1) and (2) in the reproducing kernel Hilbert space $V_2^5[0,1].$ We need to homogenize the conditions to apply the reproducing kernel Hilbert space method. We use the following transformation to homogenize these conditions.

$$\begin{array}{ccc}\hfill v\left(z\right)& =& u\left(z\right)-f\left(z\right)\hfill \end{array}$$

where

$$\begin{array}{ccc}\hfill f\left(z\right)& =& \frac{a^{2}b{b}_1+2a^{2}b{b}_2-a^{2}z{b}_1-2a^{2}z{b}_2-2a{b}^2{b}_1}{6(a-b)}\hfill \\ \hfill & -& \frac{a{b}^2{b}_2+2abz{b}_1-2az{b}_2+3a{z}^2{b}_2+2b^{2}z{b}_1}{6(a-b)}\hfill \\ \hfill & +& \frac{z{b}^2{b}_2-3b{z}^2{b}_1+z^3{b}_1-z^3{b}_2+6a{a}_2}{6(a-b)}-\frac{6b{a}_1+6z{a}_1-6z{a}_2}{6(a-b)}.\hfill \end{array}$$

After using the above transformation we obtain the following problem:

$$Hv=v^{\left(4\right)}\left(z\right)+D^{\gamma }p\left(z\right)v\left(z\right)=M\left(z\right),z\in [a,b],$$

(8)

subject to the end conditions

$$v\left(a\right)=v\left(b\right)=v^{″}\left(a\right)=v^{″}\left(b\right)=0$$

(9)

where

$$\begin{array}{c}\hfill M\left(z\right)=g\left(z\right)-D^{\gamma }p\left(z\right)f\left(z\right).\end{array}$$

Lemma 3.

H is a bounded linear operator.

Proof. 

We will show

$${∥Hv∥}_{V_2^1[0,1]}^2\le P{∥v∥}_{V_2^5[0,1]}^2,$$

where P is a positive constant. We know

$${∥Hv∥}_{V_2^1[0,1]}^2={〈Hv,Hv〉}_{V_2^1[0,1]}={\left[Hv\left(0\right)\right]}^2+{\int }_0^1{\left[H{v}^{\prime }\left(z\right)\right]}^{2}dz.$$

We obtain

$$v\left(y\right)={〈v(·),A_y(·)〉}_{V_2^5[0,1]}$$

and

$$Hv\left(y\right)={〈v(·),H{A}_y(·)〉}_{V_2^5[0,1]}.$$

Therefore, we reach

$$\left|Hv\right|\le {∥v∥}_{V_2^5[0,1]}{∥H{A}_y∥}_{V_2^5[0,1]}=P_1{∥v∥}_{V_2^5[0,1]}.$$

Thus, we acquire

$${\left[Hv\left(0\right)\right]}^2\le P_1^2{∥v∥}_{V_2^5[0,1]}^2.$$

Since

$${\left(Hv\right)}^{\prime }\left(z\right)={〈v(·),{\left(H{A}_y\right)}^{\prime }(·)〉}_{V_2^5[0,1]},$$

we obtain

$$\left|{\left(Hv\right)}^{\prime }\right|\le {∥v∥}_{E_2^5[0,1]}{∥{\left(H{A}_y\right)}^{\prime }∥}_{E_2^5[0,1]}=P_2{∥v∥}_{V_2^5[0,1]}.$$

Then, we find

$${\left[Hv\right]}^2\le P_2^2{∥v∥}_{V_2^5[0,1]}^2.$$

At last we obtain

$$\begin{array}{ccc}\hfill {∥Hv∥}_{V_2^1[0,1]}^2& =& {\left[Hv\left(0\right)\right]}^2+{\int }_0^1{\left[{\left(Hv\right)}^{\prime }\left(z\right)\right]}^{2}dz\le (P_1^2+P_2^2){∥v∥}_{V_2^5[0,1]}^2,\hfill \end{array}$$

where $P=P_1^2+P_2^2$ is a positive constant. This completes the proof. □

We construct $a_i\left(z\right)=U_{z_i}\left(z\right)$ and ${\psi }_i\left(z\right)=H^{*}{a}_i\left(z\right)$, where $H^{*}$ is conjugate operator of H. The orthonormal system ${\left\{{\widehat{\psi }}_i\left(z\right)\right\}}_{i=1}^{\infty }$ of $E_2^5[0,1]$ can be obtained by Gram–Schmidt orthogonalization operation of ${\left\{{\psi }_i\left(z\right)\right\}}_{i=1}^{\infty }$,

$${\widehat{\psi }}_i\left(z\right)=\sum _{k=1}^i{\beta }_{ik}{\psi }_k\left(z\right),({\beta }_{ii}>0,i=1,2,\dots ).$$

(10)

Theorem 2.

If $v\left(z\right)$ is the exact solution of (8) and (9), then we obtain

$$v\left(z\right)=\sum _{i=1}^{\infty }\sum _{k=1}^i{\beta }_{ik}M\left(z_k\right){\widehat{\psi }}_i\left(z\right).$$

(11)

where ${\left\{z_i\right\}}_{i=1}^{\infty }$ is dense in $[0,1]$.

Proof. 

We prove this theorem by using the reproducing property, the features of adjoint operator and the complete system as:

$$\begin{array}{ccc}\hfill v\left(z\right)& =& \sum _{i=1}^{\infty }{〈v\left(z\right),{\widehat{\psi }}_i\left(z\right)〉}_{V_2^5[0,1]}{\widehat{\psi }}_i\left(z\right)\hfill \\ \hfill & =& \sum _{i=1}^{\infty }\sum _{k=1}^i{\beta }_{ik}{〈v\left(z\right),{\psi }_k\left(z\right)〉}_{V_2^5[0,1]}{\widehat{\psi }}_i\left(z\right)\hfill \\ \hfill & =& \sum _{i=1}^{\infty }\sum _{k=1}^i{\beta }_{ik}{〈v\left(z\right),H^{*}{a}_k\left(z\right)〉}_{V_2^5[0,1]}{\widehat{\psi }}_i\left(z\right)\hfill \\ \hfill & =& \sum _{i=1}^{\infty }\sum _{k=1}^i{\beta }_{ik}{〈Hv\left(z\right),a_k\left(z\right)〉}_{V_2^1[0,1]}{\widehat{\psi }}_i\left(z\right)\hfill \\ \hfill & =& \sum _{i=1}^{\infty }\sum _{k=1}^i{\beta }_{ik}{〈M\left(z\right),U_{z_k}〉}_{V_2^1[0,1]}{\widehat{\psi }}_i\left(z\right)\hfill \\ \hfill & =& \sum _{i=1}^{\infty }\sum _{k=1}^i{\beta }_{ik}U\left(z_k\right){\widehat{\psi }}_i\left(z\right).\hfill \end{array}$$

This completes the proof. □

The approximate solution $v_n\left(z\right)$ can be obtained as:

$$v_n\left(z\right)=\sum _{i=1}^n\sum _{k=1}^i{\beta }_{ik}M\left(z_k\right){\widehat{\psi }}_i\left(z\right).$$

(12)

4. Numerical Experiments

Two test examples have been taken into consideration to demonstrate the accuracy of the reproducing kernel Hilbert space method in this section. The numerical simulation has been executed in MAPLE 18.

Example 1.

We consider the fourth order fractional boundary value problem as [17]:

$$u^{\left(4\right)}\left(z\right)+0.05D^{\gamma }u\left(z\right)=g\left(z\right),z\in [0,1],$$

(13)

with the boundary conditions

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

(14)

We have the exact solution of the problem as $z^4(z-1)$.

The absolute errors for different values of γ applying the reproducing kernel Hilbert space method are given in the

Table 1. Absolute errors for Example 1.

x	$\mathit{\gamma }=0.99$	$\mathit{\gamma }=1.0$
$0.0$	$0.000000000000000$	$0.000000000000000$
$0.1$	$1.172251330\times {10}^{-4}$	$1.079370100\times {10}^{-4}$
$0.2$	$2.050389600\times {10}^{-4}$	$1.892279000\times {10}^{-4}$
$0.3$	$2.501585000\times {10}^{-4}$	$2.301822740\times {10}^{-4}$
$0.4$	$2.502148300\times {10}^{-4}$	$2.286283000\times {10}^{-4}$
$0.5$	$2.134363900\times {10}^{-4}$	$1.922819800\times {10}^{-4}$
$0.6$	$1.537822000\times {10}^{-4}$	$1.353536000\times {10}^{-4}$
$0.7$	$8.948229000\times {10}^{-5}$	$7.512914000\times {10}^{-5}$
$0.8$	$3.720321000\times {10}^{-5}$	$2.799920000\times {10}^{-5}$
$0.9$	$7.854290000\times {10}^{-6}$	$3.736680000\times {10}^{-6}$
$1.0$	$1.591755364\times {10}^{-8}$	$4.209367700\times {10}^{-7}$

. Figure 1, Figure 2 and Figure 3 demonstrate the approximate and exact solution for $\gamma =1.0$ and $\gamma =0.99.$

Example 2.

We take into consideration the following fourth order fractional boundary value problem [17]:

$$u^{\left(4\right)}\left(z\right)+D^{\gamma }zu\left(z\right)=g\left(z\right),z\in [0,1],$$

(15)

with the boundary conditions

$$u\left(0\right)=0,u\left(1\right)=0,u^{″}\left(0\right)=0,u^{″}\left(1\right)=26(\gamma -1).$$

(16)

We have the exact solution of the problem as $z^6(z^{\gamma }-z^{2-\gamma }).$ The absolute errors applying the reproducing kernel Hilbert space method are presented in the

Table 2. Absolute errors for Example 2.

x	$\mathit{\gamma }=0.3$
$0.0$	$1.542329002\times {10}^{-11}$
$0.1$	$1.509810086\times {10}^{-4}$
$0.2$	$2.928833663\times {10}^{-4}$
$0.3$	$4.165814337\times {10}^{-4}$
$0.4$	$5.106740380\times {10}^{-4}$
$0.5$	$5.655431300\times {10}^{-4}$
$0.6$	$5.714676400\times {10}^{-4}$
$0.7$	$5.202103100\times {10}^{-4}$
$0.8$	$4.058811200\times {10}^{-4}$
$0.9$	$2.300090700\times {10}^{-4}$
$1.0$	$2.085180231\times {10}^{-7}$

. Figure 4 and Figure 5 show numerical solutions for $\gamma =0.3.$

5. Conclusions

We studied the approximate solution of the fourth order fractional boundary value problems in this paper. We applied the reproducing kernel Hilbert space method to our problem. We demonstrated our results by tables and figures. We proved the accuracy of the reproducing kernel Hilbert space method for solutions of fourth order fractional differential equations.