Fractional variational iteration method (FVIM) is performed to give an approximate analytical solution of nonlinear fractional Riccati differential equation. Fractional derivatives are described in the Riemann-Liouville derivative. A new application of fractional variational iteration method (FVIM) was extended to derive analytical solutions in the form of a series for these equations. The behavior of the solutions and the effects of different values of fractional order α are indicated graphically. The results obtained by the FVIM reveal that the method is very reliable, convenient, and effective method for nonlinear differential equations with modified Riemann-Liouville derivative
1. Introduction
In recent years, fractional calculus used in many areas such as electrical networks, control theory of dynamical systems, probability and statistics, electrochemistry of corrosion, chemical physics, optics, engineering, acoustics, viscoelasticity, material science and signal processing can be successfully modelled by linear or nonlinear fractional order differential equations [1–10]. As it is well known, Riccati differential equations concerned with applications in pattern formation in dynamic games, linear systems with Markovian jumps, river flows, econometric models, stochastic control, theory, diffusion problems, and invariant embedding [11–17]. Many studies have been conducted on solutions of the Riccati differential equations. Some of them, the approximate solution of ordinary Riccati differential equation obtained from homotopy perturbation method (HPM) [18–20], homotopy analysis method (HAM) [21], and variational iteration method proposed by He [22]. The He’s homotopy perturbation method proposed by He [23–25] the variational iteration method [26] and Adomian decomposition method (ADM) [27] to solve quadratic Riccati differential equation of fractional order.
The variational iteration method (VIM), which proposed by He [28, 29], was successfully applied to autonomous ordinary and partial differential equations and other fields. He [30] was the first to apply the variational iteration method to fractional differential equations. In recent years, a new modified Riemann-Liouville left derivative is suggested by Jumarie [31–35]. Recently, the fractional Riccati differential equation is solved with help of new homotopy perturbation method (HPM) [23].
In this paper, we extend the application of the VIM in order to derive analytical approximate solutions to nonlinear fractional Riccati differential equation:D*αy(x)=A(x)+B(x)y(x)+C(x)y2(x),x∈R,0<α≤1,t>0,
subject to the initial conditionsy(k)(0)=dk,k=0,1,2,…,n-1,
where α is fractional derivative order, n is an integer, A(x), B(x), and C(x) are known real functions, and dk is a constant.
The goal of this paper is to extend the application of the variational iteration method to solve fractional nonlinear Riccati differential equations with modified Riemann-Liouville derivative.
The paper is organized as follows: In Section 2, we give definitions related to the fractional calculus theory briefly. In Section 3, we define the solution procedure of the fractional variational iteration method to show inefficiency of this method, we present the application of the FVIM for the fractional nonlinear Riccati differential equations with modified Riemann-Liouville derivative and numerical results in Section 4. The conclusions are then given in the final Section 5.
2. Basic Definitions
Here, some basic definitions and properties of the fractional calculus theory which can be found in [31–35].
Definition 2.1.
Assume f:R→R,x→f(x) denote a continuous (but not necessarily differentiable) function, and let the partition h>0 in the interval [0,1]. Jumarie’s derivative is defined through the fractional difference [34]:
Δ(α)=(FW-1)αf(x)=∑k=0∞(-1)k(αk)f[x+(α-k)h],
where FWf(x)=f(x+h). Fractional derivative is defined as the following limit form [1, 7]:
f(α)=limh→0Δα[f(x)-f(0)]hα.
This definition is close to the standard definition of derivatives (calculus for beginners), and as a direct result, the αth derivative of a constant, 0<α<1, is zero.
Definition 2.2.
The left-sided Riemann-Liouville fractional integral operator of order α≥0, of a function f∈Cμ,μ≥-1 is defined as
Jaαf(x)=1Γ(α)∫ax(x-τ)α-1f(τ)dτ,forα>0,x>0,Ja0f(x)=f(x).
The properties of the operator Jα can be found in [1, 7, 36].
Definition 2.3.
The modified Riemann-Liouville derivative [33, 34] is defined as
Dαx0f(x)=1Γ(n-α)dndxn∫x0(x-τ)n-α(f(τ)-f(0))dτ,
where x∈[0,1],n-1≤α<n,andn≥1.
In addition, we want to give as in the following some properties of the fractional modified Riemann-Liouville derivative.
Fractional the integration of part:
Iαba(u(α)v)=(uv)|ab-Iαba(uv(α)).
Definition 2.4.
Fractional derivative of compounded functions [33, 34] is defined as
dαf≅Γ(1+α)df,0<α<1.
Definition 2.5.
The integral with respect to (dx)α [33, 34] is defined as the solution of the fractional differential equation:
dy≅f(x)(dx)α,x≥0,y(0)=0,0<α<1.
Lemma 2.6.
Let f(x) denote a continuous function [33, 34] then the solution of (2.5) is defined as
y=∫0xf(τ)(dτ)α=α∫0x(x-τ)αf(τ)dτ,0<α<1.
For example, f(x)=xβ in (2.10) one obtains
∫0xτβ(dτ)α=Γ(β+1)Γ(α+1)Γ(α+β+1)xβ+α,0<α<1.
Definition 2.7.
Assume that the continuous function f:R→R,x→f(x) has a fractional derivative of order kα, for any positive integer k and any α, 0<α≤1; then the following equality holds, which is
f(x+h)=∑k=0∞hαkαk!fαk(x),0<α≤1.
On making the substitution h→x and x→0, we obtain the fractional Mc-Laurin series:
f(x)=∑k=0∞xαkαk!fαk(0),0<α≤1.
3. Fractional Variational Iteration Method
To describe the solution procedure of the fractional variational iteration method [31–35], we consider the following fractional Riccati differential equation:D*αy(x)=A(x)+B(x)y(x)+C(x)y2(x),x∈R,0<α≤1,t>0.
According to the VIM, we can build a correct functional for (3.1) as follows:yn+1(x)=yn(x)+Iα[λ(τ)(dαyn(τ)dτα-(A(τ)+B(τ)y(τ)+C(τ)y2(τ)))]yn+1(x)=yn(x)+1Γ(α)∫0x(x-τ)α-1λ(τ)(dαyn(τ)dτα-(A(τ)+B(τ)y(τ)+C(τ)y2(τ)))dτ.
Using (2.3), we obtain a new correction functional:yn+1(x)=yn(x)+1Γ(α+1)∫0xλ(τ)(dαyn(τ)dτα-(A(τ)+B(τ)y(τ)+C(τ)y2(τ)))(dτ)α.
It is obvious that the sequential approximations yk, k≥0 can be established by determining λ, a general Lagrange’s multiplier, which can be identified optimally with the variational theory. The function ỹn is a restricted variation which means δỹn=0. Therefore, we first designate the Lagrange multiplier λ that will be identified optimally via integration by parts. The successive approximations yn+1(x),n≥0 of the solution y(x) will be readily obtained upon using the obtained Lagrange multiplier and by using any selective function y0. The initial values are usually used for choosing the zeroth approximation y0. With λ determined, then several approximations yk, k≥0 follow immediately [37]. Consequently, the exact solution may be procured by usingy(x)=limn→∞yn(x).
4. Applications
In this section, we present the solution of two examples of the Riccati differential equations as the applicability of FVIM.
Example 4.1.
Let usconsider the fractional Riccati differential equation, we get
dαydxα=-y2(x)+1,0<α≤1,
with initial conditions:
y(0)=0.
Construction the following functional:
yn+1(x)=yn(x)+1Γ(α+1)∫0xλ(τ){dαyndτα+yn2(τ)-1}(dτ)α,
we have
δyn+1(x)=δyn(x)+1Γ(α+1)δ∫0xλ(τ){dαyndτα+yn2(τ)-1}(dτ)α=δyn+λδyn|τ=x-1Γ(α+1)∫0xdαλ(τ)dταδyn(τ)(dτ)α.
Similarly, we can get the coefficients of δyn to zero:
1+λ(τ)|τ=x=0,dαλ(τ)dτα=0.
The generalized Lagrange multiplier can be identified by the above equations:
λ(x)=-1.
substituting (4.6) into (4.3) produces the iteration formulation as follows:
yn+1(x)=yn(x)-1Γ(α+1)∫0x{dαydτα+y2(τ)-1}(dτ)α.
Taking the initial value y0(x)=0, we can derive
y1(x)=y0(x)-1Γ(α+1)∫0x{dαy0dτα+y02(τ)-1}(dτ)α=xαΓ(α+1),y2(x)=y1(x)-1Γ(α+1)∫0x{dαy1dτα+y12(τ)-1}(dτ)α=xαΓ(α+1)-1Γ(α+1)∫0x{1+x2αΓ2(α+1)-1}(dτ)α=xαΓ(α+1)-Γ(2α+1)x3αΓ2(α+1)Γ(3α+1),y3(x)=xαΓ(α+1)-Γ(2α+1)x3αΓ2(α+1)Γ(3α+1)+2Γ(2α+1)Γ(4α+1)x5αΓ3(α+1)Γ(3α+1)Γ(5α+1)-Γ2(2α+1)Γ(6α+1)x7αΓ4(α+1)Γ2(3α+1)Γ(7α+1),y4(x)=xαΓ(α+1)-Γ(2α+1)x3αΓ2(α+1)Γ(3α+1)+2Γ(2α+1)Γ(4α+1)x5αΓ3(α+1)Γ(3α+1)Γ(5α+1)-Γ2(2α+1)Γ(6α+1)x7αΓ4(α+1)Γ2(3α+1)Γ(7α+1)-4Γ(2α+1)Γ(4α+1)Γ(6α+1)x7αΓ4(α+1)Γ(3α+1)Γ(5α+1)Γ(7α+1)+2Γ(2α+1)Γ(6α+1)Γ(8α+1)x9αΓ5(α+1)Γ2(3α+1)Γ(7α+1)Γ(9α+1)+4Γ2(2α+1)Γ(4α+1)Γ(8α+1)x9αΓ5(α+1)Γ2(3α+1)Γ(5α+1)Γ(9α+1)-2Γ3(2α+1)Γ(6α+1)Γ(10α+1)x11αΓ6(α+1)Γ3(3α+1)Γ(7α+1)Γ(11α+1)-4Γ2(2α+1)Γ2(4α+1)Γ(10α+1)x11αΓ6(α+1)Γ2(3α+1)Γ2(5α+1)Γ(11α+1)+4Γ3(2α+1)Γ(4α+1)Γ(6α+1)Γ(12α+1)x13αΓ7(α+1)Γ3(3α+1)Γ(5α+1)Γ(7α+1)Γ(13α+1)-Γ4(2α+1)Γ2(6α+1)Γ(14α+1)x15αΓ8(α+1)Γ4(3α+1)Γ2(7α+1)Γ(15α+1)⋮
Then, the approximate solutions in a series form are
y(x)=Limn→∞yn(x)=xαΓ(α+1)-Γ(2α+1)x3αΓ2(α+1)Γ(3α+1)+2Γ(2α+1)Γ(4α+1)x5αΓ3(α+1)Γ(3α+1)Γ(5α+1)-Γ2(2α+1)Γ(6α+1)x7αΓ4(α+1)Γ2(3α+1)Γ(7α+1)-4Γ(2α+1)Γ(4α+1)Γ(6α+1)x7αΓ4(α+1)Γ(3α+1)Γ(5α+1)Γ(7α+1)+⋯.
As α=1 is
y(x)=Limn→∞yn(x)=x-x33+2x515-17x763+⋯.
The exact solution of (4.1) is y(x)=(e2x-1)/(e2x+1), when α=1.
Figure 1 indicates the solution obtained using FVIM versus the exact solution when α=1. Figure 2 is plotted for approximate solution of time-fractional Riccati differential equation for α=0.7, 0.8, 0.9, and 1. Equation (4.1) is solved by using the homotopy perturbation method (HPM) [24]. FVIM solutions indicate that the present algorithm performs extreme efficiency, simplicity, and reliability. The results obtained from FVIM are fully compatible with those of the HPM.
Table 1 shows the approximate solutions for (4.1) obtained for different values of α using the variational iteration method and HPM [24]. From the numerical results in Table 1, it is clear that the approximate solutions are in high agreement with the exact solutions, when α=1, and the solution continuously depends on the time-fractional derivative. Example 4.1 has been solved using HAM [21], ADM [27], VIM [26], and HPM [23–25].
Approximate solutions for (4.1).
t
Ref. [24]α=0.5
y3(x)α=0.5
[24]α=0.75
y3(x)α=0.75
[24]α=1
y3(x)α=1
y(x)exact
0
0
0
0
0
0
0
0
0.1
0.273875
0.086513
0.184795
0.190102
0.099668
0.099667
0.099667
0.2
0.454125
0.161584
0.313795
0.310033
0.197375
0.197375
0.197375
0.3
0.573932
0.238256
0.414562
0.405062
0.291312
0.291320
0.291312
0.4
0.644422
0.321523
0.492889
0.483479
0.379944
0.380005
0.379948
0.5
0.674137
0.413682
0.462117
0.550470
0.462078
0.462375
0.462117
0.6
0.671987
0.515445
0.597393
0.610344
0.536857
0.537923
0.537049
0.7
0.648003
0.626403
0.631772
0.666961
0.603631
0.606768
0.604367
0.8
0.613306
0.745278
0.660412
0.723760
0.661706
0.669695
0.664036
0.9
0.579641
0.870074
0.687960
0.783638
0.709919
0.728139
0.716297
1.0
0.558557
0.998176
0.718260
0.848783
0.746032
0.784126
0.761594
The graph indicates the solution y(x) for (4.1), when α=1.
Plots of approx. solution y3(x) for different values of α.
Example 4.2.
Let usconsider the fractional Riccati differential equation, we get
dαydxα=2y(x)-y2(x)+1,0<α≤1,
with initial conditions
y(0)=0.
Construction the following functional:
yn+1(x)=yn(x)+1Γ(α+1)∫0xλ(τ){dαyndτα-2yn(τ)+yn2(τ)-1}(dτ)α.
we have
δyn+1(x)=δyn(x)+1Γ(α+1)δ∫0xλ(τ){dαyndτα-2yn(τ)+yn2(τ)-1}(dτ)α=δyn+λδyn|τ=x-1Γ(α+1)∫0xdαλ(τ)dταδyn(τ)(dτ)α.
Similarly, we can get the coefficients of δyn to zero:
1+λ(τ)|τ=x=0,dαλ(τ)dτα=0.
The generalized Lagrange multiplier can be identified by the above equations:
λ(x)=-1,
substituting (4.16) into (4.13) produces the iteration formulation as follows:
yn+1(x)=yn(x)-1Γ(α+1)∫0x{dαyndτα-2yn(τ)+yn2(τ)-1}(dτ)α.
Taking the initial value y0(x)=0, we can derive
y1(x)=y0(x)-1Γ(α+1)∫0x{dαy0dτα-2y0(τ)+y02(τ)-1}(dτ)α=xαΓ(α+1),y2(x)=y1(x)-1Γ(α+1)∫0x{dαy1dτα+y12(τ)-1}(dτ)α=xαΓ(α+1)-1Γ(α+1)∫0x{1-2ταΓ(α+1)+τ2αΓ2(α+1)-1}(dτ)α=xαΓ(α+1)+2x2αΓ(2α+1)-Γ(2α+1)x3αΓ2(α+1)Γ(3α+1),y3(x)=xαΓ(α+1)+2x2αΓ(2α+1)-Γ(2α+1)x3αΓ2(α+1)Γ(3α+1)+4x3αΓ(3α+1)-4Γ(3α+1)x4αΓ(α+1)Γ(2α+1)Γ(4α+1)-2Γ(2α+1)x4αΓ2(α+1)Γ(4α+1)-4Γ(4α+1)x5αΓ2(2α+1)Γ(5α+1)+2Γ(2α+1)Γ(4α+1)x5αΓ3(α+1)Γ(3α+1)Γ(5α+1)-Γ2(2α+1)Γ(6α+1)x7αΓ4(α+1)Γ2(3α+1)Γ(7α+1)⋮
Then, the approximate solutions in a series form are
y(x)=Limn→∞yn(x)=xαΓ(α+1)+2x2αΓ(2α+1)-Γ(2α+1)x3αΓ2(α+1)Γ(3α+1)+4x3αΓ(3α+1)-4Γ(3α+1)x4αΓ(α+1)Γ(2α+1)Γ(4α+1)-2Γ(2α+1)x4αΓ2(α+1)Γ(4α+1)-4Γ(4α+1)x5αΓ2(2α+1)Γ(5α+1)+2Γ(2α+1)Γ(4α+1)x5αΓ3(α+1)Γ(3α+1)Γ(5α+1)-Γ2(2α+1)Γ(6α+1)x7αΓ4(α+1)Γ2(3α+1)Γ(7α+1)+⋯.
As α=1 is
y(x)=Limn→∞yn(x)=x-x33+2x515-x763+⋯.
The exact solution of (4.11) is y(x)=1+2tanh(2t+(1/2)log((2-1)/(2+1))), when α=1.
Figure 3 is plotted for approximate solution of time-fractional Riccati differential equation found in Example 4.2. In Figure 4, we have shown the graphic of approximate solution of (4.11) for α=0.7,0.8,0.9,and1. Figures 2 and 4 show that a decrease in the fractional order α causes an increase in the function.
Table 2 indicates the approximate solutions for (4.11) obtained for different values of α using the variational iteration method and HPM [24]. From the numerical results in Table 2, it is clear that the approximate solutions are in high agreement with the exact solutions, when α=1, and the solution continuously depends on the time-fractional derivative.
Approximate solutions for (4.11).
t
Ref. [24]α=0.5
y3(x)α=0.5
[24]α=0.75
y3(x)α=0.75
[24]α=1
y3(x)α=1
y(x)exact
0
0
0
0
0
0
0
0
0.1
0.321730
0.577431
0.216866
0.244460
0.110294
0.110266
0.110295
0.2
0.629666
0.912654
0.428892
0.469709
0.241965
0.241585
0.241976
0.3
0.940941
1.166253
0.654614
0.698718
0.395106
0.393515
0.395104
0.4
1.250737
1.353549
0.891404
0.924319
0.568115
0.564013
0.567812
0.5
1.549439
1.482633
1.132763
1.137952
0.757564
0.749528
0.756014
0.6
1.825456
1.559656
1.370240
1.331462
0.958259
0.945155
0.953566
0.7
2.066523
1.589984
1.594278
1.497600
1.163459
1.144826
1.152946
0.8
2.260633
1.578559
1.794879
1.630234
1.365240
1.341552
1.346363
0.9
2.396839
1.530028
1.962239
1.724439
1.554960
1.527690
1.526911
1.0
2.466004
1.448805
2.087384
1.776542
1.723810
1.695238
1.689498
The graph indicates the solution y(x) for (4.11), when α=1.
Plots of approx. solution y3(x) for different values of α.
Example 4.3.
Let usconsider the fractional Riccati differential equation, we get
dαydxα=x2+y2(x),0<α≤1,
with initial conditions:
y(0)=1.
Construction the following functional:
yn+1(x)=yn(x)+1Γ(α+1)∫0xλ(τ){dαyndτα-τ2-yn2(τ)}(dτ)α.
we have
δyn+1(x)=δyn(x)+1Γ(α+1)δ∫0xλ(τ){dαyndτα-τ2-yn2(τ)}(dτ)α=δyn+λδyn|τ=x-1Γ(α+1)∫0xdαλ(τ)dταδyn(τ)(dτ)α.
Similarly, we can get the coefficients of δyn to zero:
1+λ(τ)|τ=x=0,dαλ(τ)dτα=0.
The generalized Lagrange multiplier can be identified by the above equations:
λ(x)=-1.
substituting (4.26) into (4.23) produces the iteration formulation as follows:
yn+1(x)=yn(x)-1Γ(α+1)∫0x{dαyndτα-τ2-yn2(τ)}(dτ)α.
Taking the initial value y0(x)=1, we can derive
y1(x)=y0(x)-1Γ(α+1)∫0x{dαy0dτα-τ2-y02(τ)}(dτ)α=1+xαΓ(α+1)+2x2+αΓ(α+3),y2(x)=y1(x)-1Γ(α+1)∫0x{dαy1dτα-y12(τ)-τ2}(dτ)α=1+xαΓ(α+1)+2x2+αΓ(α+3)+2x2αΓ(2α+1)+4x2+2αΓ(2α+3)+Γ(2α+1)x3αΓ2(α+1)Γ(3α+1)+4Γ(2α+3)x2+3αΓ(α+1)Γ(α+3)Γ(3α+3)+4Γ(2α+5)x4+3αΓ2(α+3)Γ(3α+5),y3(x)=1+xαΓ(α+1)+2x2+αΓ(α+3)+2x2αΓ(2α+1)+4x2+2αΓ(2α+3)+Γ(2α+1)x3αΓ2(α+1)Γ(3α+1)+4Γ(2α+3)x2+3αΓ(α+1)Γ(α+3)Γ(3α+3)+4Γ(2α+5)x4+3αΓ2(α+3)Γ(3α+5)+8x2+3αΓ(3α+3)+16Γ(5α+3)Γ(2α+3)x2+6αΓ(2α+1)Γ(α+1)Γ(α+3)Γ(3α+3)Γ(6α+3)+16Γ(3α+5)x4+4αΓ(α+3)Γ(2α+3)Γ(4α+5)+4Γ(4α+1)x5αΓ2(2α+1)Γ(5α+1)+2Γ(2α+1)x4αΓ2(α+1)Γ(4α+1)+2Γ(2α+1)Γ(4α+1)x5αΓ3(α+1)Γ(3α+1)Γ(5α+1)+8Γ(3α+3)x2+4αΓ(α+1)Γ(α+3)Γ(3α+4)+8Γ(3α+3)x2+4αΓ(α+3)Γ(2α+1)Γ(4α+3)+4x3αΓ(3α+1)+Γ2(2α+1)Γ(4α+1)x5αΓ4(α+1)Γ2(3α+1)Γ(7α+1)+4Γ(5α+1)x6αΓ(3α+1)Γ2(α+1)Γ(6α+1)+8Γ(2α+5)x4α+4Γ2(α+3)Γ(4α+5)+4Γ(3α+1)x4αΓ(α+1)Γ(2α+1)Γ(4α+1)+16Γ(4α+5)x4+5αΓ2(2α+3)Γ(5α+5)+8Γ(2α+3)x2+4αΓ(α+1)Γ(α+3)Γ(4α+3)+8Γ(2α+3)Γ(4α+3)x5α+2Γ2(α+1)Γ(α+3)Γ(3α+3)Γ(5α+3)+8Γ(2α+5)Γ(4α+5)x6α+2Γ(α+1)Γ2(α+3)Γ(3α+5)Γ(5α+5)+4Γ(2α+1)Γ(4α+3)x5α+2Γ(3α+1)Γ(α+3)Γ(3α+1)Γ(5α+3)+16Γ(2α+3)Γ(4α+5)x5α+4Γ(α+1)Γ2(α+3)Γ(3α+3)Γ(5α+5)+8Γ(2α+1)Γ(5α+3)x6α+2Γ(3α+1)Γ2(α+1)Γ(2α+3)Γ(6α+3)+16Γ(2α+5)Γ(5α+5)x6α+4Γ(2α+1)Γ2(α+3)Γ(3α+5)Γ(6α+5)+32Γ(5α+5)x6α+4Γ(α+1)Γ(α+3)Γ(3α+3)Γ(6α+5)+32Γ(2α+5)Γ(5α+7)x6α+6Γ(2α+3)Γ2(α+3)Γ(3α+5)Γ(6α+7)+16Γ2(2α+3)Γ(6α+5)x7α+4Γ2(α+1)Γ2(α+3)Γ2(3α+3)Γ(7α+5)+16Γ(4α+3)x2+5αΓ(2α+1)Γ(2α+3)Γ(5α+3)+16Γ(2α+5)Γ(4α+7)x5α+6Γ3(α+3)Γ(3α+5)Γ(5α+7)+16Γ2(2α+5)Γ(6α+9)x8+7αΓ4(α+3)Γ(7α+9)Γ2(3α+5)+8Γ(2α+1)Γ(2α+3)Γ(6α+3)x7α+2Γ3(α+1)Γ(3α+1)Γ(3α+3)Γ(α+3)Γ(7α+3)+8Γ(2α+1)Γ(2α+5)Γ(6α+5)x7α+4Γ2(α+1)Γ(3α+1)Γ(3α+5)Γ2(α+3)Γ(7α+5)+32Γ(2α+3)Γ(2α+5)Γ(6α+7)x7α+6Γ(α+1)Γ(3α+3)Γ(3α+5)Γ3(α+3)Γ(7α+7)⋮
Then, the approximate solutions in a series form are
y(x)=Limn→∞yn(x)=1+xαΓ(α+1)+2x2+αΓ(α+3)+2x2αΓ(α+1)Γ(2α+1)+4x2+2αΓ(2α+3)+Γ(2α+1)x3αΓ2(α+1)Γ(3α+1)+4Γ(2α+3)x2+3αΓ(α+1)Γ(α+3)Γ(3α+3)+4Γ(2α+5)x4+3αΓ2(α+3)Γ(3α+5)+⋯.
As α=1 is
y(x)=Limn→∞yn(x)=1+x+x2+5x36+5x46+8x515+29x690+47x7315+41t8360+299t911340+4t10525+184t1151975+t122268+4t1312285+t1559535+⋯.
The exact solution of (4.21) is
y(x)=t(J-3/4(t2/2)Γ(1/4)+2J3/4(t2/2)Γ(3/4))J1/4(t2/2)Γ(1/4)-2J-1/4(t2/2)Γ(3/4),
where Jv(t) is the Bessel function of first kind, when α=1.
Figure 5 is plotted for approximate solution of time-fractional Riccati differential equation found in Example 4.3. In Figure 6, we have shown the graphic of approximate solution of (4.21) for α=0.5,0.65,0.75,and1. Figures 2, 4, and 6 show that a decrease in the fractional order α causes an increase in the function.
Table 3 indicates the approximate solutions for (4.21) obtained for different values of α using the HPM [23]. From the numerical results in Table 3, it is clear that the approximate solutions are in substantial agreement with the exact solutions, when α=1, and the solution continuously depends on the time-fractional derivative.
Approximate solutions for (4.21).
t
Ref. [23]α=0.50
y3(x)α=0.5
[23]α=0.75
y3(x)α=0.75
[23]α=1
y3(x)α=1
y(x)exact
0
1
1
1
1
1
1
1
0.1
1.412871
1.671055
1.219154
1.249863
1.111463
1.110922
1.111450
0.2
1.693287
2.342131
1.426528
1.513463
1.252997
1.248193
1.252533
0.3
2.035965
3.241501
1.684735
1.852218
1.439293
1.420818
1.435450
0.4
2.500705
4.468508
2.031414
2.306699
1.692749
1.641742
1.675200
0.5
3.161671
6.144389
2.517316
2.928996
2.047355
1.929432
1.989583
0.6
4.121732
8.429666
3.215023
3.789643
2.553642
2.310101
2.399200
0.7
5.523246
11.54066
4.227199
4.985354
3.284790
2.820861
2.927450
0.8
7.558999
15.77047
5.695700
6.649324
4.343940
3.514150
3.600533
0.9
10.48397
21.51700
7.811804
8.965376
5.872787
4.463992
4.447450
1.0
14.62818
29.32098
10.82770
12.18753
8.061507
5.774817
5.500000
The graph indicates the solution y(x) for (4.21), when α=1.
Plots of approx. solution y3(x) for different values of α.
5. Conclusions
In this paper, variational iteration method having integral w.r.t. (dτ)α has been successfully implemented to finding approximate analytical solution of fractional Riccati differential equations. Variational iteration method known as very powerful and an effective method for solving fractional Riccati differential equation. It is also a promising method to solve other nonlinear equations. In this paper, we have discussed modified variational iteration method having integral w.r.t. (dτ)α used for the first time by Jumarie. The obtained results indicate that this method is powerful and meaningful for solving the nonlinear fractional differential equations. Three examples indicate that the results of variational iteration method having integral w.r.t. (dτ)α are in excellent agreement with those obtained by HPM, ADM, and HAM, which is available in the literature.
OldhamK. B.SpanierJ.1974New York, NY, USAAcademic PressPodlubnyI.1999New York, NY, USAAcademic PressKilbasA. A.SrivastavaH. M.TrujilloJ. J.2006Amsterdam, The NetherlandsElsevierPodlubnyI.1999San Diego, Calif, USAAcademic PressCaputoM.Linear models of dissipation whose Q is almost frequency independent, part II196713552953910.1111/j.1365-246X.1967.tb02303.xKilbasA. A.SrivastavaH. H.TrujilloJ. J.2006Amsterdam, The NetherlandsElsevierMillerK. S.RossB.1993New York, NY, USAWileySamkoS. G.KilbasA. A.MarichevO. I.1993Yverdon, SwitzerlandGordon and BreachZaslavskyG. M.2005Oxford University PressMerdanM.YıldırımA.GökdoğanA.Numerical solution of time-fraction Modified Equal Width Wave EquationEngineering Computations. In pressEinickeG. A.WhiteL. B.BitmeadR. R.The use of fake algebraic Riccati equations for co-channel demodulation2003519228822932-s2.0-004291290510.1109/TSP.2003.815376GerberM.HasselblattB.KeesingD.The riccati equation: pinching of forcing and solutions20031221291342-s2.0-0344927291KalmanR. E.HoY. C.NarendraK. S.Controllability of linear dynamical systems19631189213BittantiS.ColaneriP.De NicolaoG.The periodic Riccati equation1991Berlin, GermanySpringer127162Communications and Control EngineeringBittantiS.ColaneriP.GuardabassiG. O.Periodic solutions of periodic Riccati equations19842976656672-s2.0-0021455225AndersonB. D. O.MooreJ. B.1979Englewood Cliffs, NJ, USAPrentice-HallReidW. T.197286New York, NY, USAAcademic PressAminikhahH.HemmatnezhadM.An efficient method for quadratic Riccati differential equation20101548358392-s2.0-7035035497210.1016/j.cnsns.2009.05.009AbbasbandyS.Homotopy perturbation method for quadratic Riccati differential equation and comparison with Adomian's decomposition method200617214854902-s2.0-3114447933010.1016/j.amc.2005.02.014AbbasbandyS.Iterated He's homotopy perturbation method for quadratic Riccati differential equation200617515815892-s2.0-3364488577010.1016/j.amc.2005.07.035TanY.AbbasbandyS.Homotopy analysis method for quadratic Riccati differential equation20081335395462-s2.0-3484888079810.1016/j.cnsns.2006.06.006AbbasbandyS.A new application of He's variational iteration method for quadratic Riccati differential equation by using Adomian's polynomials2007207159632-s2.0-3425061549010.1016/j.cam.2006.07.012KhanN. A.AraA.JamilM.An efficient approach for solving the Riccati equation with fractional orders2011619268326892-s2.0-7995547435410.1016/j.camwa.2011.03.017OdibatZ.MomaniS.Modified homotopy perturbation method: application to quadratic Riccati differential equation of fractional order20083611671742-s2.0-3534886986110.1016/j.chaos.2006.06.041CangJ.TanY.XuH.LiaoS. J.Series solutions of non-linear Riccati differential equations with fractional order2009401192-s2.0-6454914882810.1016/j.chaos.2007.04.018JafariH.TajadodiH.He’s variational iteration method for solving fractional Riccati differential equation20102010876473810.1155/2010/764738MomaniS.ShawagfehN.Decomposition method for solving fractional Riccati differential equations20061822108310922-s2.0-3375151221510.1016/j.amc.2006.05.008HeJ. H.Variational iteration method—a kind of non-linear analytical technique: some examples19993446997082-s2.0-0000092673HeJ. H.WuX. H.Variational iteration method: new development and applications2007547-88818942-s2.0-3474887067710.1016/j.camwa.2006.12.083HeJ. H.Some applications of nonlinear fractional differential equations and their approximations199915286902-s2.0-0033102037JumarieG.Stochastic differential equations with fractional Brownian motion input19932461113113110.1080/00207729308949547JumarieG.New stochastic fractional models for Malthusian growth, the Poissonian birth process and optimal management of populations2006443-42312542-s2.0-3374508222810.1016/j.mcm.2005.10.003JumarieG.Laplace's transform of fractional order via the Mittag-Leffler function and modified Riemann-Liouville derivative20092211165916642-s2.0-7034921207210.1016/j.aml.2009.05.011JumarieG.Table of some basic fractional calculus formulae derived from a modified Riemann-Liouville derivative for non-differentiable functions20092233783852-s2.0-5704918653810.1016/j.aml.2008.06.003JumarieG.On the solution of the stochastic differential equation of exponential growth driven by fractional Brownian motion20051878178262-s2.0-1814437844710.1016/j.aml.2004.09.012JangM.-J.ChenC.-L.LiuY.-C.Two-dimensional differential transform for partial differential equations20011212-32612702-s2.0-003587725110.1016/S0096-3003(99)00293-3FarazN.KhanY.JafariH.YildirimA.MadaniM.Fractional variational iteration method via modified Riemann-Liouville derivative20102344134172-s2.0-7795664186410.1016/j.jksus.2010.07.025