A fuzzy solution of nonlinear partial differential equations

Author(s): Mawia Osman1, Zengtai Gong2, Altyeb Mohammed Mustafa1,3
1College of Mathematics and Statistics, Northwest Normal University, Lanzhou, Gansu, P.R. China.
2College of Mathematics and Statistics, Northwest Normal University, Lanzhou, Gansu, P.R. China
3Department of Applied Mathematics, Faculty of Mathematical Science, University of Khartoum, Khartoum, Sudan.
Copyright © Mawia Osman, Zengtai Gong, Altyeb Mohammed Mustafa. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

In this paper, the reduced differential transform method (RDTM) is applied to solve fuzzy nonlinear partial differential equations (PDEs). The solutions are considered as infinite series expansions which converge rapidly to the solutions. Some examples are solved to illustrate the proposed method.

Keywords: Fuzzy numbers, fuzzy-valued functions, gH-differentiability, reduced differential transform method, fuzzy partial differential equations.

1. Introduction

The fuzzy sets were introduced for the first time by Zadeh in [1]. Hundreds of examples have been supplied where the nature of uncertainty in the behavior of given system processes is fuzzy rather than stochastic nature. Recently, many authors showed interest in the study of the theoretical framework of fuzzy initial value problems. Chang and Zadeh [2] introduced the concept of fuzzy derivative. Dubosi and Prade [3] presented the extension principle. The differential and integral calculus for fuzzy-set-valued functions, shortly fuzzy-valued functions was developed in resent work, see [4,5,6,7,8,9].

It is known that phenomena of nature or physical systems can be modeled using partial differential equations (PDEs) such as wave equations, heat equations, Poisson’s equation and so on. Hence, studies of PDEs become one of the main topics of modern mathematical analysis and attracted much attention. Many authors developed different methods for solving different kinds of PDEs, see [10,11,12,13,14,15,16,17].

The differential transform method (DTM) was introduced by Zhou [18] and he applied it to solve initial value problems for electric circuit analysis. The DTM is based on Taylor,s series expansion and can be applied to solve both linear and nonlinear ordinary differential equations as well as PDEs. Keskin and Oturanc [19] proposed the RDTM, defining a set of transformation rules to overcome the complicated complex calculations of traditional DTM. Recently some researchers used RDTM for solving different equations, see [20,21,22,23,24,25,26].

This paper is structured as follows: In Section 2, we call some definition on a fuzzy number, fuzzy-valued function and strongly generalized Hukuhara differentiability. In Section 3, Taylor’s formula, one-dimensional DTM, and two-and three-dimensional RDTM is introduced. In Section 4, we provide some examples to show the efficiency and simplicity of RDTM. Finally, Section 5 consists of some brief conclusions.

2. Basic concepts

The fuzzy set u~E1 is called a fuzzy number if u~ is a normal, convex fuzzy set, upper semi-continuous and suppu={xR|u(x)>0} is compact. Here A¯ denotes the closure of A. We use E1 to denote the fuzzy number space [27,28].

For u~,v~E1, kR, the addition and scalar multiplication are defined by

[u~+v~]r=[u~]r+[v~]r,[ku~]r=k[u~]r, respectively, where [u~]r={x:u(x)r}=[u_r,u¯r] for any r[0,1].

We use the Hausdorff distance between fuzzy numbers [28] D:E1×E1[0,+) defined by

D(u~+v~)=supr[0,1]d([u~]r[v~]r)=supr[0,1]max{|u_rv_r|,|u¯rv¯r|}, where d is the Hausdorff metric. D(u~,v~) is called the distance between u~ and v~.

Definition 1. [29,30] Let u~ be a fuzzy number defined in F(R). The r-level set of u~, for any r[0,1] denoted by u~r is a crisp set that contains all elements in R, such that the membership value of u~ is greater or equal to r, that is u~r={xR|u~(x)r}. Whenever we represent the fuzzy number with r-level set, we mean that it is closed and bounded and it is denoted by [u_r,u¯r], where they represent the lower and upper bound r-level set of a fuzzy number.

The researchers [31,32] defined the parametrical representation of the fuzzy numbers as in the following definition:

Definition 2. [33] A fuzzy number u~ in parametric form is a pair [u_r,u¯r] of functions u_r and u¯r for any r[0,1], which satisfies the following requirements

  • u_r is a bounded non-decreasing left continuous function in (0,1];
  • u¯r is a bounded non-increasing left continuous function in (0,1];
  • u_ru¯r.
Some researchers classified the fuzzy numbers into several types of the fuzzy membership function and the triangular fuzzy membership function or also often referred to as triangular fuzzy number is the most widely used membership function.

In order to avoid the inconvenience, in the whole paper, the fuzzy numbers and fuzzy-valued functions are represented with a tilde sign at the top, while the real-value function and interval-valued functions are written directly.

Definition 3. [34] A fuzzy valued function f~ of two variable is a rule that assigns to each ordered pair of real numbers, (x,t), in a set D, a unique fuzzy number denoted by f~(x,t). The set D is the domain of f~ and its range is the set of values taken by f~, i.e., {f~(x,t)|(x,t)D}.

The parametric representation of the fuzzy valued function f:DE1 is expressed by f~(x,t;r)=[f_(x,t;r),f¯(x,t;r)], for all (x,t)D and α[0,1].

Definition 4. [34,35] A fuzzy valued function f:DE1 is said to be fuzzy continuous at (x0,t0)D if lim(x,t)(x0,t0)f~(x,t)=f(x0,t0). We say that f~ is fuzzy continuous on D if f~ is fuzzy continuous at every point (x0,t0) in D.

Definition 5. [36,37] The generalized Hukuhara difference of two fuzzy numbers u~,v~E1 is defined as follows:

u~gHv~=w~{ (i) u~=v~w~,or (ii) v~=u~(w~).
(1)
In terms of the rlevels, we get [u~gHv~]=[min{u_rv_r,u¯rv¯r},max{u_rv_r,u¯rv¯r}] and if the H-difference exists, then u~v~=u~gHv~; the conditions for existence of w~=u~gHv~E1 are
Case (i)   w_r=u_rv_r and w¯r=u¯rv¯r,r[0,1],with w_r increasing, w¯r decreasing, w_rw¯r.
(2)
Case (ii)  w_r=u¯rv¯r and w¯r=u_rv_r,r[0,1],with w_r increasing, w¯r decreasing, w_rw¯r.
(3)
It is easy to show that (i) and (ii) are both valid if and only if w~ is a crisp number. In this case, it is possible that the gH-difference of two fuzzy numbers does not exist. To address this shortcoming, a new difference of fuzzy numbers was presented in [37]

Definition 6. [38] Let u~(x,t):DE1 and (x0,t)D. We say that u~ is strongly generalized Hukuhara differentiable on (x0,t) (gH-differentiable for short) if there exists an element u~x|(x0,t)E1 such that

  • (i) &nbsp for all h>0 sufficiently small, u~(x0+h,t)gHu~(x0,t),u~(x0,t)gHu~(x0h,t) and the limits (in the metric D) limh0+u~(x0+h,t)gHu~(x0,t)h=limh0+=u~(x0,t)gHu~(x0h,t)h=u~xgH|(x0,t), or
  • &nbsp (ii) for all h>0 sufficiently small, u~(x0,t)gHu~(x0+h,t),u~(x0h,t)gHu~(x0,t) and the limits limh0+u~(x0,t)gHu~(x0+h,t)h=limh0+u~(x0h,t)gHu~(x0,t)h=u~xgH|(x0,t), or
  • (iii) &nbsp for all h>0 sufficiently small, u~(x0+h,t)gHu~(x0,t),u~(x0h,t)gHu~(x0,t) and the limits limh0+u~(x0+h,t)gHu~(x0,t)h=limh0+u~(x0h,t)gHu~(x0,t)h=u~xgH|(x0,t), or
  • (iv) &nbsp for all h>0 sufficiently small, u~(x0,t)gHu~(x0+h,t),u~(x0,t)gHu~(x0h,t) and the limits limh0+u~(x0,t)gHu~(x0+h,t)h=limh0+u~(x0,t)gHu~(x0h,t)h=u~xgH|(x0,t).

Lemma 1. [39] Let u~(x,t):DE1. Then the following statements hold:

  • (a) &nbsp If u~(x,t) is (i)-partial differentiable for x (i.e., u~ is partial differentiable for x under the meaning of Definition 5 (i), similarly to t), then
    [u~x]r=[u_(x,t)(r)x,u¯(x,t)(r)x],
    (4)
  • (b) &nbsp If u~(x,t) is (ii)-partial differentiable for x (i.e., u~ is partial differentiable for x under the meaning of Definition 5 (ii), similarly to t), then
    [u~x]r=[u¯(x,t)(r)x,u_(x,t)(r)x].
    (5)

Remark 1. For u~(x,t):DE1, the following results hold:

[2u~x2]r=[2u_(x,t)(r)x2,2u¯(x,t)(r)x2],
(6)
if (i,i),(ii,ii)-2u~x2 exist, and
[2u~x2]r=[2u¯(x,t)(r)x2,2u_(x,t)(r)x2],
(7)
if (i,ii),(ii,i)-2u~t2 exist.

3. Analysis of the method

In this section, we shall give some definitions and theorems of the Taylor series, one-dimensional DTM, and two-and three-dimensional RDTM.

Definition 7. [40] A Taylor series for the polynomial of degree n is defined as

Fn(x)=j=0n1j!(fj(c))(xc)j.
(8)

Theorem 1. If the function f(x) has (n+1) derivatives on an interval (cr,c+r) for some r>0, and limnRn(x)=0, for all x(cr,c+r), where Rn(x) is the error between Fn(x) and the polynomial function f(x) then the Taylor series expanded about x=c converges to f(x). Thus

f(x)=j=0n1j!(fj(c))(xc)j,for all,x(cr,c+r).
(9)

3.1. Differential transform method

We consider the following one-dimensional DTM:

Definition 8. [26] The differential transform F(j) of the function f(x) for jth order derivative is defined as

F(j)=1j!(djf(x)dxj)x=x0.
(10)

Definition 9. [26] The inverse differential transform of F(j) is defined as

f(x)=j=0(xx0)jF(j).
(11)
The Equation (9) is the Taylor series expansion of f(x) at x=x0. From Equations (10) and (11), the following basic operations of DTM can be deduced
  • &nbsp If f(x)=y1(x)±y2(x),  then  F(j)=Y1(j)±Y2(j).
  • &nbsp If f(x)=ay1(x), then  F(j)=aY1(j),where a is a constant.
  • &nbsp If f(x)=dy1(x)dx  then  F(j)=(j+1)Y1(j+1).
  • &nbsp If f(x)=d2y1(x)dx2  then  F(j)=(j+1)(j+2)Y1(j+2).
  • &nbsp If f(x)=dny1(x)dxn,  then  F(j)=(j+1)!j!Y1(j+1).
  • &nbsp If f(x)=y1(x)y2(x),  then  F(x)=i=0jY1(i)Y2(ji).
  • &nbsp If f(x)=xn,  then  F(j)=μ(jn)  where μ(jn)={1,   j=n,0,   jn.
  • &nbsp If f(x)=(1+x)n  then  F(j)=n(n1)(nj+1)j!.
  • &nbsp If f(x)=eax,  then  F(j)=ajj!,where a is a constants.
  • &nbsp If f(x)=sin(ωx+α),then F(j)=ωjj!sin(jπ2+α),  where ω and α are constants.
  • &nbsp If f(x)=cos(ωx+α),thenF(j)=ωjj!cos(jπ2+α),  where ω and α are constants.

3.2. Reduced differential transform method

We consider the following two-dimensional RDTM:

Definition 10. [41] If the function w(x,t) is analytical and differentiable continuously with respect to time t and space x in the domain of interest, then we get

Wj(x)=1j![jw(x,t)xj]t=0,
(12)
where the t-dimensional spectrum function Wj(x) is the transformed function of w(x,t). Here the lower case function w(x,t) represents the original function while the upper case Wj(x) stands for the transformed function.

Definition 11. [41] The inverse differential transform of Wj(x) is defined as

w(x,t)=j=0Wj(x)tj.
(13)
Thus combining (12) and (13), we can express the solution as
w(x,t)=j=0(1j![jw(x,t)xj]t=0)tj.
(14)
The basic concept of RDTM mainly comes from the power series expansion. For two-dimensional function
  • w(x,t)   then   Wj(x)=1j![jxju(x,t)]t=0.
  • u(x,t)=w(x,t)±v(x,y),   then   Uj(x)=Wj(x)±Vj(x).
  • u(x,t)=αw(x,t),   then   Uj(x)=αWj(x), &nbsp where α a is constant.
  • Misplaced &
  • u(x,t)=xmtnw(x,t),   then   Uj(x)=xmWjn(x).
  • u(x,t)=w(x,t)v(x,t),   then   Uj(x)=r=0jWr(x)Vjr(x)=r=0jVr(x)Wjr(x).
  • u(x,t)=rtrw(x,t),   then   Uj(x)=(j+r)!j!Wj+r(x).
  • u(x,t)=xw(x,t),   then   Uj(x)=xWj(x).
  • u(x,t)=2x2w(x,t),   then   Uj(x)=2x2Wj(x).
For three-dimensional function
  • w(x,y,t)   then   Wj(x,y)=1j![jxju(x,y,t)]t=0.
  • u(x,y,t)=w(x,y,t)±v(x,y,t)   then   Uj(x,y)=Wj(x,y)±Vj(x,y).
  • u(x,y,t)=λw(x,y,t)   then   Uj(x,y)=λWj(x,y),where λ a is constant.
  • u(x,y,t)=xmyntp   then   Uj(x)=xmynδ(jp),  where  δ(jp)={1       for   j=p,0       for  jp.
  • u(x,y,t)=xmyntpw(x,y,t)   then   Uj(x,y)=xmynWjp(x,y).
  • u(x,y,t)=w(x,y,t)v(x,y,t)   then   Uj(x,y)=r=0jWr(x,y)Vjr(x,y)=r=0jVr(x,y)Wjr(x,y).
  • u(x,y,t)=rtrw(x,y,t)   then    Uj(x,y)=(j+r)!j!Wj+1(x,y).
  • u(x,y,t)=xw(x,y,t)   then   Uj(x,y)=xWj(x,y).
  • u(x,y,t)=2x2w(x,y,t)   then   Uj(x,y)=2x2Wj(x,y).
  • u(x,y,t)=2y2w(x,y,t)   then   Uj(x,y)=2y2Wj(x,y).

4. Examples

In this section, we demonstrate how RDTM can be easily applied to obtain the exact solutions of the fuzzy partial differential equations.

Example 1. Consider the following one-dimensional initial value problem describing fuzzy heat-like equations

u~t=(x25)2u~x2,     0<x0,
(15)
subject to the initial condition
u~(x,0)=[rn,(2r)n](x2+x),
(16)
where (n=1,2,3,). Now
k~(s)={ s,              s[0,1], 2s          s(1,2], 0               s[0,2],
(17)
[k~n_](r)=rn and [k~n¯](α)=(2r)n. The parametric form of (15) is
u_t=(x2+5)2u_x2,     0<x0,
(18)
u¯t=(x2+5)2u¯x2,     0<x0,
(19)
for r[0,1], where u_ stands for u_(x,t;r), and u¯ stands for u¯(x,t;r). Applying the RDTM on Equations (18) and (19), we get the recurrence relation as
{(j+1)U_j+1(x,t;r)=(x2+5)2x2U_j(x;r),(j+1)U¯j+1(x,t;r)=(x2+5)2x2U¯j(x;r),
(20)
where U~j(x;r)=[U_j(x;r),U¯j(x;r)] is the transform function. From the initial condition (16), we get
U_0(x;r)=rn(x2+x),
(21)
U¯0(x;r)=(2r)n(x2+x).
(22)
Substituting U~0(x;r)=[U_0(x;r),U¯0(x;r)] into the recurrence relation (20), we get the following U~j(x;r) values successively U_1(x;r)=2(x2+5)rnU_2(x;r)=22(x2+5)rn2!U¯3(x;r)=23(x2+5)rn3!U¯4(x;r)=24(x2+5)rn4!} and U_1(x;r)=2(x2+5)(2r)nU_2(x;r)=22(x2+5)(2r)n2!U¯3(x;r)=23(x2+5)(2r)n3!U¯4(x;r)=24(x2+5)(2r)n4!} The inverse differential transform of U~j(x;r) is obtained from the relations u_(x,t;r)=j=0U_j(x;r)tj,u¯(x,t;r)=j=0U¯j(x;r)tj, and the exact solution is u~(x,t;r)=[rn,(2r)n]((x2+5)e2t+x5),0r1.

Example 2. Consider the following two-dimensional initial value problem describing fuzzy heat-like equations

u~t=12(x22u~x2y22u~y2),     0<x,y0,
(23)
subject to the initial condition
u~(x,y,0)=[(0.2+0.2r)n,(0.60.2r)n](x2+y2),
(24)
where (n=1,2,3,). Now
k~(s)={ 5(s0.2),        s[0.2,0.4], 5(0.6s),        s(0.4,0.6], 0                     s[0.2,0.6],
(25)
[k~n_](r)=(0.2+0.2r)n and [k~n¯](α)=(0.60.2r)n. The parametric form of (23) is
u_t=12(x22u_x2+y22u_y2),     0<x0,
(26)
u¯t=12(x22u¯x2+y22u¯y2),     0<x0,
(27)
for r[0,1], where u_ stands for u_(x,t)(r) and u¯ stands for u¯(x,t)(r),. Applying the RDTM on Equations (26) and (27), we get the recurrence relation as
{(j+1)U_j+1(x,y;r)=12(x22x2U_j(x,y;r)+y22y2U_j(x,y;r)),(j+1)U¯j+1(x,y;r)=12(x22x2U¯j(x,y;r)+y22y2U¯j(x,y;r)),
(28)
where U~j(x,y;r) is the transform function. From the initial condition (24), we get U_0(x,y;r)=(0.2+0.2r)n+(x2+y2),U¯0(x,y;r)=(0.60.2r)n+(x2+y2). Substituting U~0(x,y;r) into the recurrence relation (28), we get the following U~j(x,y;r) values successively U_1(x,y;r)=(0.2+0.2r)n+(x2+y2)U_2(x,y;r)=(0.2+0.2r)n+(x2+y2)2!U_3(x,y;r)=(0.2+0.2r)n+(x2+y2)3!U_4(x,y;r)=(0.2+0.2r)n+(x2+y2)4!} and U¯1(x,y;r)=(0.60.2r)n+(x2+y2)U¯2(x,y;r)=(0.60.2r)n+(x2+y2)2!U¯3(x,y;r)=(0.60.2r)n+(x2+y2)3!U¯4(x,y;r)=(0.60.2r)n+(x2+y2)4!} The solution for u~(x,t;r) is u_(x,y,t;r)=j=0U_j(x,y;r)tj=(0.2+0.2r)n+[(x2+y2)(1+t+t22!+t33!+t44!+)],u¯(x,y,t;r)=j=0U¯j(x,y;r)tj=(0.60.2r)n+[(x2+y2)(1+t+t22!+t33!+t44!+)], and the exact solution is u~(x,y,t;r)=[(0.2+0.2r)n,(0.60.2r)n]((x2+y2)exp(t)),0r1.

Example 3. We consider following two-dimensional initial value problem describing heat-like equations

u~t=ν~(x,y,t)14[x22u~x2y22u~y2],0<x,y0,
(29)
subject to the initial condition
u(x,y,0)=0~,
(30)
where ν~(x,y,t;r)=(1,0,1)(xy)2=[(r1)n,(1r)n](xy)2,0r1,(n=1,2,3,),0~E1. The parametric form of (29) is u_t=(r1)n(xy)2+14[x22u_x2+y22u_y2],0<x,y0,u¯t=(1r)n(xy)2+14[x22u¯x2+y22u¯y2],0<x,y0. Applying the RDTM, we get the recurrence relation as
{(j+1)U_j+1(x,y;r)=(r1)nx2y2δ(j)+14[x22x2U_j(x,y;r)+y22y2U_j(x,y;r)],(j+1)U¯j+1(x,y;r)=(1r)nx2y2δ(j)+14[x22x2U¯j(x,y;r)+y22y2U¯j(x,y;r)],
(31)
where δ(j)=1 when j=0, and δ(j)=0 when j0. Moreover U~j(x,y;r)=[U_j(x,y;r),U¯j(x,y;r)] is the transform function. From the initial conditions, we obtain
U_0(x,y;r)=0~,
(32)
U¯0(x,y;r)=0~.
(33)
Substituting U~0(x,y;r) into the recurrence relation (31), we get the following U~j(x,y;r) values successively U_1(x,y;r)=(r1)nx2y2U_2(x,y;r)=(r1)nx2y22!U_3(x,y;r)=(r1)nx2y23!U_4(x,y;r)=(r1)nx2y24!} and U¯1(x,y;r)=(1r)nx2y2U¯2(x,y;r)=(1r)nx2y22!U¯3(x,y;r)=(1r)nx2y23!U¯4(x,y;r)=(1r)nx2y24!} The solution for u~(x,t;r) is u_(x,y,t;r)=j=0U_j(x,y,t;r)tj=(r1)n[(x2y2)(1+t+t22!+t33!+t44!+)],u¯(x,y,t;r)=j=0U¯j(x,y,t;r)tj=(1r)n[(x2y2)(1+t+t22!+t33!+t44!+)], and the exact solution is u~(x,y,t;r)=[(r1)n,(1r)n]((x2y2)exp(t)),0r1.

Example 4. We consider the following fuzzy partial differential equation

u~t=2u~x2gHu~2,
(34)
subject to the initial condition
u~(x,0)=[(0.5+0.5r)n,(1.50.5r)n]1,
(35)
where (n = 1,2,3,…). Now
k~(s)={ 2(s0.5),s[0.5,1], 2(1.5s),s(1,1.5], 0,s[0.5,1.5],
(36)
[k~n_](r)=(0.5+0.5r)n and [k~n¯](α)=(1.50.5r)n. The parametric form of (34) is
u_t=2u_x2u_2,
(37)
u¯t=2u¯x2u¯2,
(38)
for r[0,1], where u_ stands for u_(x,t)(r) and u¯ stands for u¯(x,t)(r). Applying the RDTM on Equations (37) and (38), we get the recurrence relation as
{(k+1)U_j+1(x;r)=2x2[U_j(x;r)]i=1j[U_i(x;r)U_ji(x;r)],(k+1)U¯j+1(x;r)=2x2[U¯j(x;r)]i=1j[U¯i(x;r)U¯ji(x;r)].
(39)
The transformed initial condition (35) becomes
U_0(x;r)=(0.5+0.5r)+1,
(40)
U¯0(x;r)=(1.50.5r)+1.
(41)
For different values of j, we get the following results U_1(x;r)=(0.5+0.5r)1U_2(x;r)=(0.5+0.5r)+1U_3(x;r)=(0.5+0.5r)1U_4(x;r)=(0.5+0.5r)+1} and U¯1(x;r)=(1.50.5r)1U¯2(x;r)=(1.50.5r)+1U¯3(x;r)=(1.50.5r)1U¯4(x;r)=(1.50.5r)+1} The solution for u~(x,t;r) is u_(x,t;r)=j=0U_j(x;r)tj=(U_0(x;r)+U_1(x;r)t+U_2(x;r)t2+U_3(x;r)t3+)=(0.5+0.5r)+(1t+t2t3+),u¯(x,t;r)=j=0U¯j(x;r)tj=(U¯0(x;r)+U¯1(x;r)t+U¯2(x;r)t2+U¯3(x;r)t3+)=(1.50.5r)+(1t+t2t3+), and the exact solution is u~(x,t;r)=[(0.5+0.5r)n,(1.50.5r)n][11+t],0r1. Figure 1 illustrate that the left-hand functions of the r-level set of u~ (u lower) are always increasing functions of r and the right-hand functions of the r-level set of u~ (u upper) are always decreasing functions of r in the above examples.

5. Conclusion

In this paper, the reduced differential transform method (RDTM) has been successfully applied for solving fuzzy nonlinear partial differential equations under gH-differentiability. The solutions are considered as infinite series expansions that converge rapidly to the exact solutions. We solved some examples to illustrate the proposed method. The results reveal that the proposed method is a powerful and efficient technique for solving fuzzy nonlinear partial differential equations.

Author Contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References:

  1. Zaden, L. A. (1965). Fuzzy sets. lnformation and Control, 8, 338-353. [Google Scholor]
  2. Chang, S. S., & Zadeh, L. A. (1996). On fuzzy mapping and control. In Fuzzy Sets, Fuzzy Logic, and Fuzzy Systems: selected papers by Lotfi A Zadeh (pp. 180-184). [Google Scholor]
  3. Dubois, D., & Prade, H. (1982). Towards fuzzy differential calculus part 3: Differentiation. Fuzzy sets and Systems, 8(3), 225-233. [Google Scholor]
  4. Dubois, D., & Prade, H. (1982). Towards fuzzy differential calculus part 1: Integration of fuzzy mappings. Fuzzy Sets and Systems, 8(1), 1-17. [Google Scholor]
  5. Dubois, D., & Prade, H. (1982). Towards fuzzy differential calculus part 2: Integration on fuzzy intervals. Fuzzy Sets and Systems, 8(2), 105-116. [Google Scholor]
  6. Dubois, D., & Prade, H. (1982). Towards fuzzy differential calculus part 3: Differentiation. Fuzzy Sets and Systems, 8(3), 225-233. [Google Scholor]
  7. Puri, M. L., & Ralescu, D. A. (1983). Differentials of fuzzy functions. Journal of Mathematical Analysis and Applications, 91(2), 552-558. [Google Scholor]
  8. Puri, M. L., & Ralescu, D. (1983). Differential of fuzzy functions. Journal of Mathematical Analysis and Applications, 91, 321-325. [Google Scholor]
  9. Goetschel Jr, R., & Voxman, W. (1986). Elementary fuzzy calculus. Fuzzy Sets and Systems, 18(1), 31-43. [Google Scholor]
  10. Kharab, A., & Kharab, R. (1997). Spreadsheet solution of hyperbolic partial differential equations [for EM field calculations]. IEEE Transactions on Education, 40(1), 103-110. [Google Scholor]
  11. Kincaid, D., Kincaid, D. R., & Cheney, E. W. (2009). Numerical Analysis: Mathematics of Scientific Computing,(Vol. 2). American Mathematical Socity. [Google Scholor]
  12. Biro, O., & Preis, K. (1990). Finite element analysis of 3-D eddy currents. IEEE Transactions on Magnetics, 26(2), 418-423. [Google Scholor]
  13. McDonald, B. H., & Wexler, A. (1972). Finite-element solution of unbounded field problems. IEEE Transactions on Microwave Theory and Techniques, 20(12), 841-847.[Google Scholor]
  14. Gratkowski, S., Pichon, L., & Razek, A. (1996). New infinite elements for a finite element analysis of 2D scattering problems. IEEE Transactions on Magnetics, 32(3), 882-885. [Google Scholor]
  15. Osman, M., Gong, Z. T., & Mohammed, A. (2021). Differential transform method for solving fuzzy fractional wave equation. Journal of Computational Analysis and Applications, 29(3), 431-453.[Google Scholor]
  16. Osman, M., Gong, Z., & Mustafa, A. M. (2020). Comparison of fuzzy Adomian decomposition method with fuzzy VIM for solving fuzzy heat-like and wave-like equations with variable coefficients. Advances in Difference Equations, 2020, Article No 327. [Google Scholor]
  17. Osman, M., Gong, Z., & Mustafa, A. M. (2020). Applcation to the fuzzy fracrtional diffusion equation by using fuzzy fractional variational Homotopy pertubation iteration method. Advance Research Journal of Multidisciplinary Discoveries, 51(3), 15-26. [Google Scholor]
  18. Zhou, J. K. (1986). Differential Transformation and its Application for Electcuits Circuits. Huarjung University press, Wuhan, China. [Google Scholor]
  19. Keskin, Y., Servi, S., & Oturanç, G. (2011). Reduced differential transform method for solving Klein Gordon equations. In Proceedings of the World Congress on Engineering (Vol. 1, p. 2011). [Google Scholor]
  20. Arshad, M., Lu, D., & Wang, J. (2017). (N+ 1)-dimensional fractional reduced differential transform method for fractional order partial differential equations. Communications in Nonlinear Science and Numerical Simulation, 48, 509-519. [Google Scholor]
  21. Yu, J., Jing, J., Sun, Y., & Wu, S. (2016). (n+1)Dimensional reduced differential transform method for solving partial differential equations. Applied Mathematics and Computation, 273, 697-705. [Google Scholor]
  22. Keskin, Y., & Oturanç, G. (2010). Application of reduced differential transformation method for solving gas dynamics equation. International Journal of Contemporary Mathematical Sciences, 22(22), 1091-1096. [Google Scholor]
  23. Cenesiz, Y., Keskin, Y., & Kurnaz, A. (2010). The solution of the nonlinear dispersive K(m,n) equations by RDT method. International Journal of Nonlinear Science, 9,(4), 461-467. [Google Scholor]
  24. Raslan, K. R., Biswas, A., & Sheer, Z. F. A. (2012). Differential transform method for solving partial differential equations with variable coefficients. International Journal of Physical Sciences, 7(9), 1412-1419. [Google Scholor]
  25. Taha, B. A., & Abdul-Wahab, R. D. (2012). Numarical solutions of two-dimensional Burgers’ equations using reduced differential tranform method. Journal of Thi-Qar Science, 3(3), 252-264. [Google Scholor]
  26. Neog, B. C. (2015). Solutions of some system of non-linear PDEs using reduced differential transform method. IOSR Journal of Mathematics, 11(5), 37-44. [Google Scholor]
  27. Gong, Z., & Wang, L. (2012). The Henstock–Stieltjes integral for fuzzy-number-valued functions. Information Sciences, 188, 276-297. [Google Scholor]
  28. Negoita, C. V., & Ralescu, D. A. (1975). Applications of fuzzy sets to systems analysis (p. 187). Basel, Switzerland: Birkhäuser. Applications+of+fuzzy+sets+to+systems+analysis&btnG=” target=”_blank”>[Google Scholor]
  29. Kaleva, O. (2006). A note on fuzzy differential equations. Nonlinear Analysis: Theory, Methods & Applications, 64(5), 895-900. [Google Scholor]
  30. Kaleva, O. (1987). Fuzzy differential equations. Fuzzy Sets and Systems, 24(3), 301-317. [Google Scholor]
  31. Ma, M., Friedman, M., & Kandel, A. (1999). Numerical solutions of fuzzy differential equations. Fuzzy Sets and Systems, 105(1), 133-138. [Google Scholor]
  32. Friedman, M., Ma, M., & Kandel, A. (1999). Numerical solutions of fuzzy differential and integral equations. Fuzzy Sets and Systems, 106(1), 35-48. [Google Scholor]
  33. Rahman, N. A. A., & Ahmad, M. Z. (2017). Solving fuzzy fractional differential equations using fuzzy Sumudu transform. Journal of Nonlinear Sciences and Applications, 10(5), 2620-2632.[Google Scholor]
  34. Allahviranloo, T., Gouyandeh, Z., Armand, A., & Hasanoglu, A. (2015). On fuzzy solutions for heat equation based on generalized Hukuhara differentiability. Fuzzy Sets and Systems, 265, 1-23. [Google Scholor]
  35. Guang-Quan, Z. (1991). Fuzzy continuous function and its properties. Fuzzy Sets and Systems, 43(2), 159-171. [Google Scholor]
  36. Yang, H., & Gong, Z. (2019). Ill-posedness for fuzzy Fredholm integral equations of the first kind and regularization methods. Fuzzy Sets and Systems, 358, 132-149. [Google Scholor]
  37. Bede, B., & Stefanini, L. (2013). Generalized differentiability of fuzzy-valued functions. Fuzzy Sets and Systems, 230, 119-141. [Google Scholor]
  38. Bede, B., & Gal, S. G. (2005). Generalizations of the differentiability of fuzzy-number-valued functions with applications to fuzzy differential equations. Fuzzy Sets and Systems, 151(3), 581-599. [Google Scholor]
  39. Chalco-Cano, Y., & Roman-Flores, H. (2008). On new solutions of fuzzy differential equations. Chaos, Solitons & Fractals, 38(1), 112-119. [Google Scholor]
  40. Neog, B. C. (2015). Solutions of some system of non-linear PDEs using reduced differential transform method. IOSR Journal of Mathematics, 11(5), 37-44. [Google Scholor]
  41. Keskin, Y., & Oturanc, G. (2009). Reduced differential transform method for partial differential equations. International Journal of Nonlinear Sciences and Numerical Simulation, 10(6), 741-750. [Google Scholor]
  42. Mustafa, A. M., Gong, Z., & Osman, M. (2020). The solution of fuzzy variational problem and fuzzy optimal control problem under granular differentiability concept. International Journal of Computer Mathematics, https://doi.org/10.1080/00207160.2020.1823974.[Google Scholor]