# On the Fractional Nagumo Equation with Nonlinear Diffusion and Convection

Authors: Abdon Atangana; Suares Clovis Oukouomi Noutchie
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

We presented the Nagumo equation using the concept of fractional calculus. With the help of two analytical techniques including the homotopy decomposition method (HDM) and the new development of variational iteration method (NDVIM), we derived an approximate solution. Both methods use a basic idea of integral transform and are very simple to be used.


---

## Body

## 1. Introduction

The Nagumo equation with linear or nonlinear diffusion and convection has broadly been useful to population dynamics, ecology, neurophysiology, chemical reactions, and flame propagation [1–3]. In particular, the case where the equation involves degenerate nonlinear diffusion is of considerable interest [4–8]. In this case, a travelling wave front solution of sharp type is known to exist for exactly one value of the wave speed. Such wave fronts, for instance, represent collective motion of populations in particular collective cell spreading, invasion in ecology, and concentration in chemical reactions [9–15]. However, it has been showing that the real world problems describe via fractional order derivative gives better prediction [16–20]. It is, therefore, important; further extend the nonlinear Nagumo equation using the concept of fractional derivative order.It is something very difficult to obtain the exact solution of nonlinear equation with fractional order derivative. Many scholars sometimes, to avoid this difficulty, solve this class of problem numerically. However, even with numerical scheme it is also difficult to provide a numerical solution for nonlinear equations. Thus, many scholars, to access the behaviour of the solution of the real problem under study, present an approximate solution of this type of equations.In the literature, there exist several analytical techniques [21–25] to deal with nonlinear equations including fractional type. The purpose of this work is to present an approximate solution for the generalized nonlinear Nagumo equation withnonlinear diffusion and convection via the well-known variational iteration method (VIM) and the homotopy decomposition method (HDM). The nonlinear fractional Nagumo equation considered here is given below as
(1)∂αu∂tα+βun∂u∂x=∂∂x[αun∂u∂x]+γu(1-um)(um-δ),0<α≤1,(2)u(x,0)=f(x),(3)u(0,t)=g(t),
where α,β,γ, and δ are constants, and for the sake of simplicity in this paper, we consider m=n=1=δ. The concept of fractional order is not well by some scholars; in order to accommodate those, we present in the next section the basic information regarding this concept.


## 2. Basic Information about Fractional Calculus

There exists a vast literature on different definitions of fractional derivatives. The most popular ones are the Riemann-Liouville and the Caputo derivatives. For Caputo we have the following.Definition 1 (see [26–30]).
A real functionf(x),x>0, is said to be in the space Cμ, μ∈R if there exists a real number p>μ, such that f(x)=xph(x), where h(x)∈C[0,∞), and it is said to be in space Cμm if f(m)∈Cμ, m∈N.Definition 2 (partial derivatives of fractional order [31–34]).
Assume now thatf(x) is a function of n variables xi,  i=1,…,n also of class C on D∈Rn. We define partial derivative of order α for f with respect to xi, the function
(4)a∂x_αf=1Γ(m-α)∫axi(xi-t)m-α-1∂ximf(xj)|xj=tdt,
where ∂xim is the usual partial derivative of integer order m.Definition 3 (see [26–30]).
The Riemann-Liouville fractional integral operator of orderα≥0, of a function f∈Cμ, μ≥-1, is defined as
(5)Jαf(x)=1Γ(α)∫0x(x-t)α-1f(t)dt,α>0,x>0,J0f(x)=f(x).
Properties of the operator can be found in [31–38]; we mention only the following.For(6)f∈Cμ,μ≥-1,α,β≥0,γ>-1,JαJβf(x)=Jα+βf(x),JαJβf(x)=JβJαf(x),Jαxγ=Γ(γ+1)Γ(α+γ+1)xα+γ.


## 3. Basic Information of the HDM and VIM

In this section, we shall present the basic information of the chosen analytical methods: the homotopy decomposition method and variational iteration method; we shall start with the homotopy decomposition method.


### 3.1. Information about the HDM

The interested reader can find the full detail of the methodology of the homotopy decomposition method in [39–42]. This relatively new method was recently used to solve some nonlinear fractional partial differential equations. However, since the new development of the variational iteration method using the Laplace transform was recently introduced, we shall present its methodology in the following subsection.


### 3.2. Some Information about the Variational Iteration Method

In its initial development, the essential nature of the method was to construct the following correction functional for (2) when α is a natural number:
(7)wn+1(x,t)=wn(x,t)+∫0tλ(t,τ)[-∂mw(x,τ)∂tm+L(wn(x,τ))mmmmmmmm∂mw(x,τ)∂tm+N(wn(x,τ))+k(x,τ)]dτ,
where λ(t,τ) is the so-called Lagrange multiplier [43] and wn(x,t) is the n-approximate solution. However, this development was not suitable for equations with fractional order derivative [43]. Therefore, in their work, they apply the new development of the VIM proposed in [44] to find the Lagrange multiplier. In this new VIM, the first step of the basic character of the method is to apply the Laplace transform on both sides of (2) to obtain
(8)smw(x,s)-sm-1w(x,0)-⋯wm-1(x,0)=L[L(w(x,t))+N(w(x,t))+k(x,t)].
The recursive formula of (8) can now be used to put forward the main recursive method connecting the Lagrange multiplier as
(9)wn+1(x,s)=wn(x,s)+λ(s)[smwn(x,s)-sm-1w(x,0)-⋯wm-1(x,0)mmmmmmm-L[L(wn(x,t))+N(wn(x,t))+k(x,t)]sm-1w(x,0)].
Now considering L[L(wn(x,t))+N(wn(x,t))+k(x,t)], the restricted term; the Lagrange multiplier can be obtained as [43]
(10)λ(s)=-1sm.
Now, applying the inverse Laplace transform on both sides of (9), we obtain the following iteration:
(11)wn+1(x,t)=wn(x,t)-L-1[1sm[-L[L(wn(x,t))+N(wn(x,t))+k(x,t)]]1smsmwn(x,s)-sm-1w(x,0)-⋯wm-1(x,0)hhhhhhhhhhh-L[1smL(wn(x,t))mmmminmhhhhhm+N(wn(x,t))+k(x,t)1sm]]].


## 3.1. Information about the HDM

The interested reader can find the full detail of the methodology of the homotopy decomposition method in [39–42]. This relatively new method was recently used to solve some nonlinear fractional partial differential equations. However, since the new development of the variational iteration method using the Laplace transform was recently introduced, we shall present its methodology in the following subsection.


## 3.2. Some Information about the Variational Iteration Method

In its initial development, the essential nature of the method was to construct the following correction functional for (2) when α is a natural number:
(7)wn+1(x,t)=wn(x,t)+∫0tλ(t,τ)[-∂mw(x,τ)∂tm+L(wn(x,τ))mmmmmmmm∂mw(x,τ)∂tm+N(wn(x,τ))+k(x,τ)]dτ,
where λ(t,τ) is the so-called Lagrange multiplier [43] and wn(x,t) is the n-approximate solution. However, this development was not suitable for equations with fractional order derivative [43]. Therefore, in their work, they apply the new development of the VIM proposed in [44] to find the Lagrange multiplier. In this new VIM, the first step of the basic character of the method is to apply the Laplace transform on both sides of (2) to obtain
(8)smw(x,s)-sm-1w(x,0)-⋯wm-1(x,0)=L[L(w(x,t))+N(w(x,t))+k(x,t)].
The recursive formula of (8) can now be used to put forward the main recursive method connecting the Lagrange multiplier as
(9)wn+1(x,s)=wn(x,s)+λ(s)[smwn(x,s)-sm-1w(x,0)-⋯wm-1(x,0)mmmmmmm-L[L(wn(x,t))+N(wn(x,t))+k(x,t)]sm-1w(x,0)].
Now considering L[L(wn(x,t))+N(wn(x,t))+k(x,t)], the restricted term; the Lagrange multiplier can be obtained as [43]
(10)λ(s)=-1sm.
Now, applying the inverse Laplace transform on both sides of (9), we obtain the following iteration:
(11)wn+1(x,t)=wn(x,t)-L-1[1sm[-L[L(wn(x,t))+N(wn(x,t))+k(x,t)]]1smsmwn(x,s)-sm-1w(x,0)-⋯wm-1(x,0)hhhhhhhhhhh-L[1smL(wn(x,t))mmmminmhhhhhm+N(wn(x,t))+k(x,t)1sm]]].


## 4. Application to the Fractional Nagumo Equation

In this section, we present the application of the homotopy decomposition method and the new development of the so-called variational iteration method to the nonlinear fractional Nagumo equation. We shall start with the HDM.


### 4.1. Application of the HDM

Let us consider the fractional nonlinear Nagumo Equation with the following initial condition:(12)u(x,0)=x2,∂αu∂tα+βu∂u∂x=∂∂x[au∂u∂x]+γu(1-u)(u-1),mmmmmmmmmmmmmmmmii0<α≤1.
Using the steps involved in the HDM we arrive at the following:
(13)u(x,t)=u(x,0)+1Γ[1-α]∫0t(t-τ)1-α[-βu∂u∂x+∂∂x[au∂u∂x]mmmmmmmmmmmmnmmmm[au∂u∂x]+γu(1-u)(u-1)]dτ.
Now, assume the solution of the above equation can be expressed in series form as follows:
(14)u(x,t)=∑n=0∞pnun(x,t).
Replacing this in (13) and after comparing the term of the same power of p, we obtain the following recursive formulas:
(15)u0(x,t)=u(x,t),u1(x,t)=1Γ[α]∫0t(t-τ)α-1[-βu0∂u0∂x+∂∂x[au0∂u0∂x]mmmmmmmmmmmmm+γu0(1-u0)(u0-1)[au0∂u0∂x]]dτ,u2(x,t)=1Γ[α]∫0t(t-τ)α-1[-βH21+aH22+aH23mmmmmmmmnmmm+γH24+γH25-u1]dτ.
Here,
(16)H21=u0(x,t)∂xu1(x,t)+u1(x,t)∂xu0(x,t),H22=2∂xu0(x,t)∂xu1(x,t),H23=u0(x,t)∂xxu1(x,t)+u1(x,t)∂xxu0(x,t),H24=u02(x,t)u1(x,t)+u13+u12(x,t)u0(x,t),H25=2u0(x,t)u1(x,t).
The general recursive formula for p≥3 is given as
(17)un(x,t)=1Γ[α]∫0t(t-τ)α-1[-βHn1+aHn2+aHn3mmmmmmmmmmmmm+γHn4+γHn5-un-1]dτ,
with
(18)Hn1(x,t)=∑j=0n-1uj(x,t)∂xun-j(x,t),Hn2(x,t)=∑j=0n-1∂xuj(x,t)∂xun-j(x,t),Hn3(x,t)=∑j=0n-1uj(x,t)∂xxun-j(x,t),Hn4(x,t)=∑j=0n-1∑k=0juj(x,t)uj-k(x,t)un-j-1(x,t),Hn5(x,t)=∑j=0n-1uj(x,t)un-j(x,t),
so that, integrating the above set of integral equations, we obtain the following:
(19)u0(x,t)=x2,u1(x,t)=-tαx2(-6a+2xβ+(-1+x2)2γ)Γ(1+α),u2(x,t)=t2αx2(t2α(-6a+2xβ+(-1+x2)2γ)3Γ(1+3α)Γ3(1+α)Γ(1+4α)(48a2+γmmmmmmm-2a(3+24xβ+4γ-22x2γmmmmmmmmnim+x(10xβ2hhhhhhhhhhhhhhlhhhh+xγ(-2+x2-2γ+3x2γ-x6γ)hhhhhhhhhhhhhhlhhhh+2β(1+(2-8x2+3x4)γ))))mmmmmnn×1Γ(1+2α)mmmmmin+tαx4γ(-6a+2xβ+(-1+x2)2γ)2Γ(1+2α)Γ2(1+α)Γ(1+3α)mmmmmin-t2α(-6a+2xβ+(-1+x2)2γ)3Γ(1+3α)Γ3(1+α)Γ(1+4α)).
Here, we have computed only three terms in the series solution. However, using the recursive formula, we can compute the remaining terms, and the approximate solution is given as follows:
(20)u(x,t)=u0(x,t)+u1(x,t)+u2(x,t)+⋯+⋯.


### 4.2. New Development of the Variational Iteration Method

In this subsection, we test the efficiency of the new development of the variation iteration method by solving the nonlinear fractional equation (1). Therefore, following the methodology of NDVIM, we are at the following.The Lagrange multiplier is(21)λ(s)=-1sα
and the recursive formula is given as
(22)un+1(x,t)=un(x,t)-L-1[1sα[+γun(1-un)(un-1)∂un∂x]L[-βun∂un∂xmmmmmmmmmmmm+∂∂x[aun∂un∂x]mmmmmmmmmmmm+γun(1-un)(un-1)∂un∂x]]]
with the initial term
(23)u0(x,t)=x2,
so that using the iteration formulas we obtain
(24)u1(x,t)=x2-tαΓ(1+α)×(6ax2-2x3β+x2(1-x2)(x2-1)γ),(25)u2(x,t)=(x2(t3α(-144a3mmmmmnmm+12a2mmmmmnmm×(18xβ+4γ-22x2γ+21x4γ)mmmmmmnm-4a(12xβγ-54x3βγmmmmmmmmmmn+48x5βγ+γ2mmmmmmmmnmm+43x4γ2-48x6γ2+18x8γ2mmmmmmmmmnm+2x2(12β2-7γ2))mmmmm+x(2βγ2+52x4βγ2-52x6βγ2mmmmimmm+18x8βγ2+26x7γ3-14x9γ3mmmimmmm+3x11γ3+4x2(3β3-5βγ2)mmimmmmm+6x5(5β2γ-4γ3)+2x(5β2γ-γ3)mmimmmmm+x3(-36β2γ+11γ3)))mmmmm×Γ(1+α)Γ2(1+2α)Γ(1+4α)mmmmm+t2α(180a2+6xβγ-20x3βγ+14x5βγmmmmmmimm+γ2+12x4γ2-10x6γ2+3x8γ2mimmmmmmm-4a(16x4γ23xβ+9γmimmmmmmmmmm-25x2γ+16x4γ)mimmmmmmm+2x2(5β2-3γ2))mmmmm×Γ(1+3α)Γ3(1+α)Γ(1+4α)mmmmm+Γ(1+2α)Γ(1+3α)mmmmm×(t4αx4γmmmmmmm×(-6a+2xβ+γ-2x2γ+x4γ)3mmmnnmnm×Γ(1+3α)mmmnnmnm-(2tα(16a-2xβ-γ+2x2γ-x4γ)mmnmmmnnmm-Γ(1+α)x4γ)Γ(1+α)mmmnnmnm×Γ3(1+α)Γ(1+4α))))×(Γ3(1+α)Γ(1+2α)mmmmm×Γ(1+3α)Γ(1+4α)Γ3)-1.
Using the recursive formula, the remaining term can be obtained but here, due to the length of this term, we computed only three terms and the approximate solution case given as
(26)u(x,t)=u2(x,t).
In the following section, we compare the approximate solution via HDM and NDVIM.


## 4.1. Application of the HDM

Let us consider the fractional nonlinear Nagumo Equation with the following initial condition:(12)u(x,0)=x2,∂αu∂tα+βu∂u∂x=∂∂x[au∂u∂x]+γu(1-u)(u-1),mmmmmmmmmmmmmmmmii0<α≤1.
Using the steps involved in the HDM we arrive at the following:
(13)u(x,t)=u(x,0)+1Γ[1-α]∫0t(t-τ)1-α[-βu∂u∂x+∂∂x[au∂u∂x]mmmmmmmmmmmmnmmmm[au∂u∂x]+γu(1-u)(u-1)]dτ.
Now, assume the solution of the above equation can be expressed in series form as follows:
(14)u(x,t)=∑n=0∞pnun(x,t).
Replacing this in (13) and after comparing the term of the same power of p, we obtain the following recursive formulas:
(15)u0(x,t)=u(x,t),u1(x,t)=1Γ[α]∫0t(t-τ)α-1[-βu0∂u0∂x+∂∂x[au0∂u0∂x]mmmmmmmmmmmmm+γu0(1-u0)(u0-1)[au0∂u0∂x]]dτ,u2(x,t)=1Γ[α]∫0t(t-τ)α-1[-βH21+aH22+aH23mmmmmmmmnmmm+γH24+γH25-u1]dτ.
Here,
(16)H21=u0(x,t)∂xu1(x,t)+u1(x,t)∂xu0(x,t),H22=2∂xu0(x,t)∂xu1(x,t),H23=u0(x,t)∂xxu1(x,t)+u1(x,t)∂xxu0(x,t),H24=u02(x,t)u1(x,t)+u13+u12(x,t)u0(x,t),H25=2u0(x,t)u1(x,t).
The general recursive formula for p≥3 is given as
(17)un(x,t)=1Γ[α]∫0t(t-τ)α-1[-βHn1+aHn2+aHn3mmmmmmmmmmmmm+γHn4+γHn5-un-1]dτ,
with
(18)Hn1(x,t)=∑j=0n-1uj(x,t)∂xun-j(x,t),Hn2(x,t)=∑j=0n-1∂xuj(x,t)∂xun-j(x,t),Hn3(x,t)=∑j=0n-1uj(x,t)∂xxun-j(x,t),Hn4(x,t)=∑j=0n-1∑k=0juj(x,t)uj-k(x,t)un-j-1(x,t),Hn5(x,t)=∑j=0n-1uj(x,t)un-j(x,t),
so that, integrating the above set of integral equations, we obtain the following:
(19)u0(x,t)=x2,u1(x,t)=-tαx2(-6a+2xβ+(-1+x2)2γ)Γ(1+α),u2(x,t)=t2αx2(t2α(-6a+2xβ+(-1+x2)2γ)3Γ(1+3α)Γ3(1+α)Γ(1+4α)(48a2+γmmmmmmm-2a(3+24xβ+4γ-22x2γmmmmmmmmnim+x(10xβ2hhhhhhhhhhhhhhlhhhh+xγ(-2+x2-2γ+3x2γ-x6γ)hhhhhhhhhhhhhhlhhhh+2β(1+(2-8x2+3x4)γ))))mmmmmnn×1Γ(1+2α)mmmmmin+tαx4γ(-6a+2xβ+(-1+x2)2γ)2Γ(1+2α)Γ2(1+α)Γ(1+3α)mmmmmin-t2α(-6a+2xβ+(-1+x2)2γ)3Γ(1+3α)Γ3(1+α)Γ(1+4α)).
Here, we have computed only three terms in the series solution. However, using the recursive formula, we can compute the remaining terms, and the approximate solution is given as follows:
(20)u(x,t)=u0(x,t)+u1(x,t)+u2(x,t)+⋯+⋯.


## 4.2. New Development of the Variational Iteration Method

In this subsection, we test the efficiency of the new development of the variation iteration method by solving the nonlinear fractional equation (1). Therefore, following the methodology of NDVIM, we are at the following.The Lagrange multiplier is(21)λ(s)=-1sα
and the recursive formula is given as
(22)un+1(x,t)=un(x,t)-L-1[1sα[+γun(1-un)(un-1)∂un∂x]L[-βun∂un∂xmmmmmmmmmmmm+∂∂x[aun∂un∂x]mmmmmmmmmmmm+γun(1-un)(un-1)∂un∂x]]]
with the initial term
(23)u0(x,t)=x2,
so that using the iteration formulas we obtain
(24)u1(x,t)=x2-tαΓ(1+α)×(6ax2-2x3β+x2(1-x2)(x2-1)γ),(25)u2(x,t)=(x2(t3α(-144a3mmmmmnmm+12a2mmmmmnmm×(18xβ+4γ-22x2γ+21x4γ)mmmmmmnm-4a(12xβγ-54x3βγmmmmmmmmmmn+48x5βγ+γ2mmmmmmmmnmm+43x4γ2-48x6γ2+18x8γ2mmmmmmmmmnm+2x2(12β2-7γ2))mmmmm+x(2βγ2+52x4βγ2-52x6βγ2mmmmimmm+18x8βγ2+26x7γ3-14x9γ3mmmimmmm+3x11γ3+4x2(3β3-5βγ2)mmimmmmm+6x5(5β2γ-4γ3)+2x(5β2γ-γ3)mmimmmmm+x3(-36β2γ+11γ3)))mmmmm×Γ(1+α)Γ2(1+2α)Γ(1+4α)mmmmm+t2α(180a2+6xβγ-20x3βγ+14x5βγmmmmmmimm+γ2+12x4γ2-10x6γ2+3x8γ2mimmmmmmm-4a(16x4γ23xβ+9γmimmmmmmmmmm-25x2γ+16x4γ)mimmmmmmm+2x2(5β2-3γ2))mmmmm×Γ(1+3α)Γ3(1+α)Γ(1+4α)mmmmm+Γ(1+2α)Γ(1+3α)mmmmm×(t4αx4γmmmmmmm×(-6a+2xβ+γ-2x2γ+x4γ)3mmmnnmnm×Γ(1+3α)mmmnnmnm-(2tα(16a-2xβ-γ+2x2γ-x4γ)mmnmmmnnmm-Γ(1+α)x4γ)Γ(1+α)mmmnnmnm×Γ3(1+α)Γ(1+4α))))×(Γ3(1+α)Γ(1+2α)mmmmm×Γ(1+3α)Γ(1+4α)Γ3)-1.
Using the recursive formula, the remaining term can be obtained but here, due to the length of this term, we computed only three terms and the approximate solution case given as
(26)u(x,t)=u2(x,t).
In the following section, we compare the approximate solution via HDM and NDVIM.


## 5. Numerical Results

We devote this section to the comparison of the numerical solutions obtained via the HDM and the NDVIM for different values of the fractional order derivative. In this case, we choseα=1.5, γ=1, and a=4. The following figures show the numerical solution of the time fractional nonlinear Nagumo equation. Figure 1 show, the approximate solution obtained via the HDM for α=0.9, Figure 2 shows the approximate solution via the NDVIM. Figure 3 Show the approximate solution obtained via the HDM for α=0.25 and Figure 4 shows the approximate solution via the NDVIM. Table 1 shows the comparison of the numerical values of the solution obtained via the HDM and the NDVIM, respectively, for different values of alpha.Table 1
Comparison of numerical values for the approximate solutions via HDM and NDVIM.
x
t
HDMα=0.25
NDVIMα=0.25
HDMα=0.9
NDVIMα=0.9
−10
0
100
100
100
100
1
-
1.14275
×
10
18
-
1.14319
×
10
18
-
3.24381
×
10
17
-
3.24545
×
10
17
2
-
2.28554
×
10
18
-
2.28628
×
10
18
-
3.93359
×
10
18
-
3.93465
×
10
18
6
-
6.85676
×
10
18
-
6.85845
×
10
18
-
2.05326
×
10
20
-
2.05347
×
10
20
10
-
1.1428
×
10
1
9
-
1.14305
×
10
1
9
-
1.29153
×
10
21
-
1.29161
×
10
2
1
−5
0
25
25
25
25
1
-
2.98096
×
10
12
-
2.99989
×
10
12
-
8.45746
×
10
1
1
-
8.52756
×
10
1
1
2
-
5.9636
×
10
12
-
5.99543
×
10
12
-
1.02663
×
10
1
3
-
1.03118
×
10
1
3
6
-
1.78972
×
10
1
3
-
1.79697
×
10
1
3
-
5.36272
×
10
1
4
-
5.37156
×
10
1
4
10
-
2.98327
×
10
1
3
-
2.99391
×
10
1
3
-
3.37376
×
10
1
5
-
3.37727
×
10
1
5
5
0
25
25
25
25
1
-
3.50932
×
10
1
2
-
3.53141
×
10
1
2
-
9.95678
×
10
1
1
-
1.00386
×
10
1
2
2
-
7.02052
×
10
1
2
-
7.05766
×
10
1
2
-
1.20856
×
10
1
3
-
1.21387
×
10
1
3
6
-
2.10687
×
10
1
3
-
2.11533
×
10
1
3
-
6.31281
×
10
1
4
-
6.32312
×
10
1
4Figure 1
Approximate solution of the time-fractional nonlinear Nagumo equation via the HDM for the value of alpha equal to 0.9.Figure 2
Approximate solution of the time-fractional nonlinear Nagumo equation via the NDVIM for the value of alpha equal to 0.9.Figure 3
Approximate solution of the time-fractional nonlinear Nagumo equation via the HDM for the value of alpha equal to 0.25.Figure 4
We present in Table1 the numerical values of the approximate solutions obtained via both methods for different values of alpha.Both methods used the idea of iteration; the initial components are obtained as the Taylor series of the exact solution. On one hand, the new development of variational iteration method makes use of the Laplace transform, the Lagrange multiplier, and finally the inverse Laplace transform. On the other hand, the HDM uses just a simple integral and the perturbation technique. Both techniques are simple to implement and are very accurate.


## 6. Conclusion

The Nagumo equation is a very complex equation, for which the exact solution does not exist. The Nagumo equation was extended to the concept of fractional order derivative. The resulting equation was further analyzed within the framework of the homotopy decomposition method and the new development of variational iteration method. Both methods use a simple idea of integral transform. The numerical results are presented to test the efficiency and the accuracy of both methods. From their iteration formulas, one can conclude that these two methods are simple to be used and are powerful weapons to handle fractional nonlinear equation type.



---

*Source: 963985-2014-08-05.xml*