Next Article in Journal
What Can We Learn from Multi-Objective Meta-Optimization of Evolutionary Algorithms in Continuous Domains?
Previous Article in Journal
Methods for Assessing Human–Machine Performance under Fuzzy Conditions
Previous Article in Special Issue
Consensus-Based Multi-Person Decision Making with Incomplete Fuzzy Preference Relations Using Product Transitivity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Switching Point Solution of Second-Order Fuzzy Differential Equations Using Differential Transformation Method

by
Nadeem Salamat
1,*,†,‡,
Muhammad Mustahsan
1,‡ and
Malik M. Saad Missen
2
1
Mathematics Department, Khwaja Fareed University of Engineering and Information Technology, Rahimyar khan-64200, Pakistan
2
Department of computer Science, Islamia University Bahawalpur, Bahawalpur 63100, Pakistan
*
Author to whom correspondence should be addressed.
Current address: Department of Computer Science and Information Technology, The Islamia University, Bahawalpur, Bahawalpur 63100, Pakistan.
These authors contributed equally to this work.
Mathematics 2019, 7(3), 231; https://doi.org/10.3390/math7030231
Submission received: 31 January 2019 / Revised: 24 February 2019 / Accepted: 25 February 2019 / Published: 1 March 2019
(This article belongs to the Special Issue Nonlinear Analysis Using Fuzzy Mathematics)

Abstract

:
The first-order fuzzy differential equation has two possible solutions depending on the definition of differentiability. The definition of differentiability changes as the product of the function and its first derivative changes its sign. This switching of the derivative’s definition is handled with the application of min, max operators. In this paper, a numerical technique for solving fuzzy initial value problems is extended to solving higher-order fuzzy differential equations. Fuzzy Taylor series is used to develop the fuzzy differential transformation method for solving this problem. This leads to a single solution for higher-order differential equations.

1. Introduction

The fuzzy set theory was pioneered by Lotfi Zadeh [1]. Fuzzy sets are widely used in a variety of mathematical modeling, and fuzzy derivative was introduced in [2,3,4,5,6,7,8,9] and is an extension of ordinary derivatives. This provides the foundation of fuzzy differential equations, and uncertain systems are modeled using fuzzy differential equations.
Fuzzy differential equations are solved using many methods, such as differential or fuzzy inclusion [10,11,12], the extension principle and characterization theorem [13,14,15], and artificial neural networks [16]. Numerical methods such as the differential transformation method [17,18], extension of numerical solution methods for ordinary differential equations (ODEs) [19,20,21,22,23,24], interactive derivatives [25], multistep methods [26], and Runge–Kutta (RK) methods [27] are also applied to solving these equations. The metric properties of fuzzy functions are studied in [28,29]. The existence of solutions for second-order fuzzy differential equations is proved in [21,30]. As a result, higher-order fuzzy differential equations can be solved by using generalized fuzzy derivatives. A system of differential equations and second-order differential equations are solved in [6,31].
A first-order differential equation with an initial value is called an initial value problem, and this sort of equation has a unique solution. If the initial values are fuzzified, then the problem is called a fuzzy initial value problem (FIVP). The differential equation with fuzzy initial values has two solutions for every fuzzification parameter. These solutions depend upon the definition of fuzzy differentiability, and the definition changes as the sign of the product of function and its first derivative change [32].
Fuzzy differential equations are interpreted as a family of differential inclusions. Fuzzy interval arithmetic is very simple to compute because it is applicable to the endpoints of an interval, and this leads to an estimate of the confidence interval. These fuzzy solutions are constructed from the solution of the crisp problem generated by the fuzzy initial value problem. Fuzzy solutions are developed using the parametric representation of fuzzy functions. In this representation, fuzzy sets are represented by lower and upper approximations: the lower approximation is increasing and the upper approximation is a decreasing function. The function’s property of increasing and decreasing is utilized with the definition of switching points. The fuzzy differentiability of type-(i) switches to that of type-(ii) at the point where the function changes its increasing nature to decreasing or vice versa. This property is also used for the first-order derivative function, and similar logic is implemented for developing the solution of second-order differential equations.
The differential transformation method (DTM) [33] is a way of finding the numerical solution of linear and nonlinear differential equations [34]. The DTM is also used for solving fuzzy differential equations (FDEs) with FIVPs [10,35]. In this paper, a solution of a second-order fuzzy differential equation is proposed. The second-order differential equation is converted to an equal system of first-order differential equations, and initial values are fuzzified. The switching points are verified for the solution and first-order derivative functions. The DTM is used to solve the system of equations, and the solution leads to the approximation of a fuzzy solution of the problem.
Two-point initial value problems are considered in [14,36,37,38], and a fuzzy boundary value problem is solved in [39]. In these approaches, the definition of differentiability in the given interval is chosen initially without switching. In this paper, the FIVP is approached by using FDTM, and the idea is extended to second-order fuzzy differential equations with the adaptability of min and max operators to handle the switching points. The solution set produces the upper and lower approximations for the bounded solution of the second-order differential equation. The first-order FIVP represents a single solution with the application of two definitions of differentiability, and for second-order FIVP, the four definitions of differentiability are combined with the application of min, max operators, where the definition of differentiability depends upon the interval of the validity of the solution and the switching point within that interval. This leads to a unique solution of the differential equation. The rest of the paper is organized as follows: First, we start with preliminaries in Section 2. The switching points and fuzzy differentiability for second-order derivatives are explained in Section 3. In Section 4, we describe the DTM, and its extension is derived in Section 4.1. In Section 5, numerical examples are discussed, and finally, Section 7 consists of the conclusion with a summary of our contribution, as well as a discussion of future directions of investigation.

2. Preliminaries

Numerous techniques and terms are used to define a fuzzy number and fuzzy membership functions [40]. In this paper, the following definitions are used.
Definition 1
(Fuzzy Set). An ordered pair of a set and a membership function is called a fuzzy set, and a membership function is defined as η : R [ 0 , 1 ] . The fuzzy set has the following properties:
1. 
A set of all the points where the membership function has a nonzero value is called the support of the fuzzy set.
2. 
The set of points where η ( x ) = 1 is called the core of a fuzzy set η.
3. 
A fuzzy set η is normal if at least one x 0 R such that η ( x 0 ) = 1 .
4. 
A fuzzy set η is convex if η ( t x + ( 1 t ) y ) min ( η ( x ) , η ( y ) ) for all x , y R and t [ 0 , 1 ] .
5. 
A fuzzy set is upper semi-continuous in R
6. 
Compactness of a fuzzy set is defined as | x R , η ( x ) > 0 | ¯ , where . ¯ stands for a closure operator.
The fuzzy real numbers, numbers with a triangular sendograph, are represented by E 1 such that R E . These fuzzy number are represented as three points, a , b , c R , with a < b < c ; then, the α -level fuzzy set is represented by h ̲ ( γ ) = a + ( b a ) γ and h ¯ ( γ ) = c ( c b ) γ . The γ cuts for all 0 < γ 1 are defined as [ h ] γ = { h ( x ) γ x R } . These are called the level sets, and the zero-level compact set is represented by [ h ] 0 = { h ( x ) > 0 for   all   x R } ¯ . The level sets are nested sets, and the γ cuts of any fuzzy number are always closed and bounded intervals.
Definition 2.
Fuzzy numbers are represented as an ordered pair of functions ( h ̲ ( γ ) , h ¯ ( γ ) ) , where γ [ 0 , 1 ] , and h ̲ ( γ ) and h ¯ ( γ ) are called the lower and upper boundaries of h ( γ ) . The h ¯ ( γ ) is bounded, left continuous, and decreasing, and h ̲ ( γ ) is bounded, right continuous, and an increasing function on [ 0 , 1 ] , and h ̲ ( γ ) h ¯ ( γ ) for all γ [ 0 , 1 ] . A crisp number is represented by an equivalent fuzzy notation with an equal approximation of lower and upper h ( γ ) . Let h and g be two fuzzy numbers defined in γ-level notation, where h = ( h ̲ , h ¯ ) and g = ( g ̲ , g ¯ ) ; then, the algebra of fuzzy numbers is defined as
Addition: 
Two fuzzy numbers h and g are added as follows: ( h + g ) ̲ ( γ ) = h ̲ ( γ ) + g ̲ ( γ ) and ( h + g ) ¯ ( γ ) = h ¯ ( γ ) + g ¯ ( γ ) , and its γ-level addition will be [ h + g ] γ = [ ( h + g ) ̲ ( γ ) , ( h + g ) ¯ ( γ ) ]
Subtraction: 
The subtraction of fuzzy triangular numbers is computed as follows: ( h g ) ̲ ( γ ) = h ̲ ( γ ) g ̲ ( γ ) and ( h g ) ¯ ( γ ) = h ¯ ( γ ) g ¯ ( γ ) . The subtraction of fuzzy numbers h g is defined as addition, and the parameterized subtraction is defined as [ h g ] γ = [ ( h g ) ̲ ( γ ) , ( h g ) ¯ ( γ ) ]
Scalar Multiplication: 
The scalar multiplication of a scalar a with a fuzzy number is computed as follows: a h ̲ ( γ ) = a h ̲ ( γ ) , a h ¯ ( γ ) = a h ¯ ( γ ) when a 0 and a h ̲ ( γ ) = a h ¯ ( γ ) , a h ¯ ( γ ) = a h ̲ ( γ ) when a < 0 , and the scalar multiplication will become a h γ = [ a h ̲ ( γ ) , a h ¯ ( γ ) ]
Multiplication: 
Let h , g be two fuzzy numbers; then, multiplication of fuzzy numbers is computed as follows: h g = [ h g ̲ , h g ¯ ] , where ( h g ) ̲ ( γ ) = min [ h ̲ ( γ ) g ̲ ( γ ) , h ¯ ( γ ) g ̲ ( γ ) , h ̲ ( γ ) g ¯ ( γ ) , h ¯ ( γ ) g ¯ ( γ ) ] and ( h g ) ¯ ( γ ) = max [ h ̲ ( γ ) g ̲ ( γ ) , h ¯ ( γ ) g ̲ ( γ ) , h ̲ ( γ ) g ¯ ( γ ) , h ¯ ( γ ) g ¯ ( γ ) ]
Division: 
Division between two fuzzy numbers is defined conditionally. If 0 [ g ̲ 0 , g ¯ 0 ] , then u v = [ ( h g ) ̲ , ( h g ) ¯ ] is defined as ( h g ) ̲ ( γ ) = min [ h ̲ ( γ ) g ̲ ( γ ) , h ¯ ( γ ) g ̲ ( γ ) , h ̲ ( γ ) g ¯ ( γ ) , h ¯ ( γ ) g ¯ ( γ ) ] and ( h g ) ¯ ( γ ) = max [ h ̲ ( γ ) g ̲ ( γ ) , h ¯ ( γ ) g ̲ ( γ ) , h ̲ ( γ ) g ¯ ( γ ) , h ¯ ( γ ) g ¯ ( γ ) ]
Definition 3
(H-Difference). The set of real numbers is defined by fuzzy membership function E, and h , g , f are three fuzzy numbers in E such that h = g + f h g = f . Then, the H-difference of two fuzzy real numbers is denoted as h g . Its γ-level representation will become
[ h g ] γ = [ h ̲ γ g ¯ γ , h ¯ γ g ̲ γ ] .
Definition 4
(g-Difference). The set of real numbers is defined by fuzzy membership function E, and h , g , f are three fuzzy numbers in E such that h = g + f h g = f . Then, the g-difference of two fuzzy real numbers is defined as x g y , and its γ cuts are
[ u g v ] γ = β γ [ h ̲ β g h g ̲ β , h ¯ β g h g ¯ β ] .
Definition 5
(Hausdorff Distance). Let F : R f × R f R + { 0 } be a function. Then, the Hausdorff distance is defined as
D ( u , v ) = S u p γ [ 0 , 1 ] max { h ̲ γ v ̲ γ , h ¯ γ v ¯ γ } = [ u ] γ g [ v ] γ ,
where a complete metric space is represented by D. The g-difference of intervals always exists. This definition is equivalent to the usual definitions for metric spaces of fuzzy numbers.
Definition 6
(Generalized Hukuhara Difference). The generalized Hukuhara difference between two fuzzy numbers is defined as a and b as follows:
[ a g H b ] γ = [ c ̲ γ , c ¯ γ ]
where c ̲ γ = min ( a ̲ γ a ̲ γ , a ¯ γ a ¯ γ ) and c ¯ γ = max ( a ̲ γ a ̲ γ , a ¯ γ a ¯ γ ) with the conditions that c ̲ γ is a nondecreasing function and c ¯ γ is a nonincreasing function such that c ̲ γ < c ¯ γ .
Definition 7
(Generalized Hukuhara Differentiability [41]). Let f be a fuzzy function and f : [ a , b ] K F ( R n ) , and let t be any point in the domain. Then, f is a strongly generalized differentiable function f at t and will be defined by the following.
for all h > 0 If h is sufficiently small and h > 0 , then
  • Differentiable type-(i): There exist f ( t + h ) f ( t ) and f ( t ) f ( t h ) ; then,
    f ( t ) = lim h 0 f ( t + h ) f ( t ) h = lim h 0 f ( t ) f ( t h ) h
    or
  • Differentiable type-(ii): If f ( t ) f ( t + h ) and f ( t h ) f ( t ) exist and
    f ( t ) = lim h 0 f ( t + h ) f ( t ) h = lim h 0 f ( t ) f ( t h ) h
for all h < 0 If h is sufficiently small and h < 0 , then
  • Differentiable type-(i): There exist f ( t + h ) f ( t ) and f ( t ) f ( t h ) ; then,
    f ( t ) = lim h 0 f ( t h ) f ( t ) h = lim h 0 f ( t ) f ( t + h ) h
  • Differentiable type-(ii): There exist f ( t ) f ( t + h ) and f ( t h ) f ( t ) ; then,
    f ( t ) = lim h 0 f ( t h ) f ( t ) h = lim h 0 f ( t ) f ( t + h ) h
Reasoning as above, there exists the Hukuhara difference f ( x 0 ) H f ( x 0 h ) . Now, if we suppose f < 0 , then we see that we cannot use the above kind of reasoning to prove that the H-differences f ( x 0 + h ) H f ( x 0 ) and f ( x 0 ) H f ( x 0 h ) and the derivative f ( x 0 ) exist.
Definition 8
(Second-Order Generalized Hukuhara Differentiability [41]). Let f be a strongly differentiable function in the generalized sense on the domain, and let there exist f ( t ) R f such that, for all h > 0 , we have
(I) 
f ( t + h ) f ( t ) and f ( t ) f ( t h ) and
f ( t ) = lim h 0 f ( t + h ) f ( t ) h = lim h 0 f ( t ) f ( t h ) h
(II) 
f ( t ) f ( t + h ) and f ( t h ) f ( t ) and
f ( t ) = lim h 0 f ( t h ) f ( t ) h = lim h 0 f ( t ) f ( t + h ) h
If f has differentiability of type-(i) and f exists, then there are two possibilities for the differentiability of function f . Similarly, if f has differentiability of type-(ii), then f has two possibilities. Thus f has four definitions of differentiability.
Definition 9
(Generalized Differentiability ([42])). Let x 0 [ a , b ] and k be such that x 0 + k [ a , b ] ; then, the g_derivative of the function f : [ a , b ] K F ( R n ) at x 0 is defined as
f g ( x 0 ) = lim k 0 1 k [ f ( x 0 + k ) g f ( x 0 ) ] .
If f g ( x 0 ) exists, then f is called differentiable in the generalized sense.
Theorem 1.
Let F : [ a , b ] K F ( R n ) be such that
[ F ( x ) ] γ = [ f ( x ) γ , g ( x ) γ ]
where ( f ( x ) γ and g ( x ) γ are differentiable real-valued functions with respect to x. Then, F ( x ) is g_ differentiable, and mathematically,
[ F ( x ) ] γ = [ min ( ( f ( x ) ) γ , ( g ( x ) ) γ ) , max ( ( f ( x ) ) γ , ( g ( x ) ) γ ) ]
The interval is divided into subintervals. The functional value at the end point of the previous interval becomes the initial value of the next subinterval and at each boundary point of the subinterval. The conditions are verified, and if there is a switching point, then the definition of fuzzy derivative changes from type-(i) to type-(ii) and vice versa. The use of min and max operators handles the switching points, which are defined in the next subsection. (Proof of Theorem 1 follows from [42].)
Proposition 1.
([43]) Let F : ( t 0 + a , T ) K F ( R n ) be Hukuhara differentiable 2 times and
[ ( F ) γ ( t ) ) ] = [ ( f 1 ) γ ( t ) ) , ( f 2 ) γ ( t ) ) ] γ [ 0 , 1 ]
such that ( f 1 ) γ ( t ) ) = min ( ( f ( x ) ) γ , ( g ( x ) ) γ ) and ( f 2 ) γ ( t ) ) = max ( ( f ( x ) ) γ , ( g ( x ) ) γ ) . If f , g are differentiable, then we can write F ( 2 ) in a similar fashion:
[ F γ ( 2 ) ( t ) ) ] = [ ( f 1 ( 2 ) ) γ ( t ) ) , ( f 2 ( 2 ) ) γ ( t ) ) ] γ [ 0 , 1 ]
In fact, the second-order derivative of a real-valued function represents the four definitions of fuzzy derivatives of fuzzy functions. The definition of the derivative changes at the switching point.
Theorem 2
(Theorem 3.6 in [41]). Let F : R E be a function and denote F ( t ) γ = [ f ( t ) γ , g ( t ) γ ] for each γ [ 0 , 1 ] . Then,
  • If F is (I)-differentiable, then f and g are differentiable functions and D 1 1 F ( t ) γ = F ( t ) γ = [ f ( t ) γ , g ( t ) γ ]
  • If f is (II)-differentiable, then f and g are differentiable functions and D 2 1 F ( t ) γ = F ( t ) γ = [ g ( t ) γ , f ( t ) γ ] .
Proof. 
For details, see [41]. □
Theorem 3
(Theorem 3.9 in [41]). Let [ f ( t ) ] = [ f ( t ) , g ( t ) ] be an interval-valued function, where f is an increasing function and g is a decreasing function, and let the system D 1 ( 1 ) f : I R f and D 2 ( 1 ) f : I R f be two fuzzy functions. Then,
1. 
Let f and g be two differentiable functions. If F and D 1 ( 1 ) F are (I)-differentiable, then D 1 , 1 ( 2 ) [ F ( t ) ] γ = [ f ( t ) γ , g ( t ) γ ] .
2. 
Let f and g be two differentiable functions. If F is (I)-differentiable and D 1 ( 1 ) F is (II)-differentiable, then D 1 , 2 ( 2 ) [ F ( t ) ] γ = [ g ( t ) γ , f ( t ) γ ] .
3. 
Let f and g be two differentiable functions. If F is (II)-differentiable and D 1 ( 1 ) F is (I)-differentiable, then D 2 , 1 ( 2 ) [ F ( t ) ] γ = [ g ( t ) γ , f ( t ) γ ] .
4. 
Let f and g be two differentiable functions. If F and D 1 ( 1 ) F are (II)-differentiable, then D 2 , 2 ( 2 ) [ F ( t ) ] γ = [ f ( t ) γ , g ( t ) γ ] .
Proof. 
For details, see [41]. □
Theorem 4
(Theorem 3.1 in [44]). Let t 0 [ a , b ] and assume that f : [ a , b ] × R F R F is continuous. A mapping x : [ a , b ] R F is a solution to the initial value problem x = f ( t ; x ; x ) and x ( t 0 ) = k 1 ; x ( t 0 ) = k 2 if and only if x and x are continuous and satisfy one of the following conditions:
  • x ( t ) = k 2 ( t t 0 ) + t 0 t ( t 0 t f ( s , x ( s ) , x ( s ) ) d s ) d s + k 1 , where x and x are (i)-differentials,
  • x ( t ) = k 1 ( 1 ) ( k 2 ( t t 0 ) ( 1 ) t 0 t ( t 0 t f ( s , x ( s ) , x ( s ) ) d s ) d s , where x and x are (ii)-differentials,
  • x ( t ) = k 1 ( 1 ) ( k 2 ( t t 0 ) + t 0 t ( t 0 t f ( s , x ( s ) , x ( s ) ) d s ) d s , where x is (i)-differentiable and x is a (ii)-differential, or
  • x ( t ) = k 2 ( t t 0 ) k 1 ( 1 ) t 0 t ( t 0 t f ( s , x ( s ) , x ( s ) ) d s ) d s + k 1 , where x is (ii)-differentiable and x is a (i)-differential.
Remark 1.
(i) This theorem shows that four systems of fuzzy differentiability represent only two systems, and two pairs are equivalent. According to the conditions of the triangular fuzzy function, f is an increasing function and g is a decreasing function. If F follows the first definition, then F = [ f , g ] . This shows that f is an increasing function and g is a decreasing function. In the other case, F = [ g , f ] , which shows that f is now decreasing and g is an increasing function. (ii) Theorem 4 is for H-differentiability, and its mode-generalized form is g H differentiability, so it is also g H differentiable.

3. Fuzzy Differentiability and Switching Points

Definition 10
(Switching Point). The points in an interval where fuzzy differentiability of type-(i) changes to type-(ii) and vice versa are called switching points.
For the first-order initial value differential equations, the above-cited switching points apply. This logic can be derived from any one of the two functions f and g, and the definition of fuzzy differentiability is associated with any of the above functions. In Table 1, f is increasing. The definition changes from type-(i) to type-(ii) if f changes its increasing nature to decreasing and vice versa.
The possibility for f is described in the Table 1. If both the function and its derivative are increasing functions (similarly for g decreasing) then, type-(i) is applicable, if f changes its nature and is a decreasing function then differentiability of type-(ii) is applicable. This behavior of function is verified over the interval. We extend the same logic to the second derivative, and we get the representations listed in Table 2.
The above table shows that the fuzzy function follows the same differentiability if there is no switch in increasing or decreasing nature of function or there are two switches. The pair with one switching point also have similar behavior. We transform the second-order differential equation into equal system of first-order differential equations.
y ( t ) = f ( t , y ( t ) , y ( t ) )
This second-order differential equation is decomposed into the system of first-order differential equation as
y 1 ( t ) = y 2 ( t ) y 2 ( t ) = f ( t , y 1 ( t ) , y 2 ( t ) )
where transformation y 1 ( t ) = y ( t ) and y 2 ( t ) = y ( t ) is used to describe the equivalent system of first-order differential equations.
The solution system described in [36] will become
  • I-I system:
    y 1 ̲ ( t , γ ) = y 2 ̲ ( t , γ ) , y 2 ̲ ( t , γ ) = f ̲ ( t , y 1 ̲ ( t , γ ) , y 1 ¯ ( t , γ ) , y 2 ̲ ( t , γ ) , y 2 ¯ ( t , γ ) ) y 1 ¯ ( t , γ ) = y 2 ¯ ( t , γ ) , y 2 ¯ ( t , γ ) = f ¯ ( t , y 1 ̲ ( t , γ ) , y 1 ¯ ( t , γ ) , y 2 ̲ ( t , γ ) , y 2 ¯ ( t , γ ) ) y 1 ̲ ( 0 , γ ) = y 1 ̲ 0 ( γ ) , y 1 ¯ ( 0 , γ ) = y 1 ¯ 0 ( γ ) y 2 ̲ ( 0 , γ ) = y 2 ̲ 0 ( γ ) , y 2 ¯ ( 0 , γ ) = y 2 ¯ 0 ( γ )
  • I-II system:
    y 1 ̲ ( t , γ ) = y 2 ̲ ( t , γ ) , y 2 ̲ ( t , γ ) = f ¯ ( t , y 1 ̲ ( t , γ ) , y 1 ¯ ( t , γ ) , y 2 ̲ ( t , γ ) , y 2 ¯ ( t , γ ) ) y 1 ¯ ( t , γ ) = y 2 ¯ ( t , γ ) , y 2 ¯ ( t , γ ) = f ̲ ( t , y 1 ̲ ( t , γ ) , y 1 ¯ ( t , γ ) , y 2 ̲ ( t , γ ) , y 2 ¯ ( t , γ ) ) y 1 ̲ ( 0 , γ ) = y 1 ̲ 0 ( γ ) , y 1 ¯ ( 0 , γ ) = y 1 ¯ 0 ( γ ) y 2 ̲ ( 0 , γ ) = y 2 ̲ 0 ( γ ) , y 2 ¯ ( 0 , γ ) = y 2 ¯ 0 ( γ )
  • II-I system:
    y 1 ̲ ( t , γ ) = y 2 ¯ ( t , γ ) , y 2 ̲ ( t , γ ) = f ̲ ( t , y 1 ̲ ( t , γ ) , y 1 ¯ ( t , γ ) , y 2 ̲ ( t , γ ) , y 2 ¯ ( t , γ ) ) y 1 ¯ ( t , γ ) = y 2 ̲ ( t , γ ) , y 2 ¯ ( t , γ ) = f ¯ ( t , y 1 ̲ ( t , γ ) , y 1 ¯ ( t , γ ) , y 2 ̲ ( t , γ ) , y 2 ¯ ( t , γ ) ) y 1 ̲ ( 0 , γ ) = y 1 ̲ 0 ( γ ) , y 1 ¯ ( 0 , γ ) = y 1 ¯ 0 ( γ ) y 2 ̲ ( 0 , γ ) = y 2 ¯ 0 ( γ ) , y 2 ¯ ( 0 , γ ) = y 2 ̲ 0 ( γ )
  • II-II system:
    y 1 ̲ ( t , γ ) = y 2 ¯ ( t , γ ) , y 2 ̲ ( t , γ ) = f ¯ ( t , y 1 ̲ ( t , γ ) , y 1 ¯ ( t , γ ) , y 2 ̲ ( t , γ ) , y 2 ¯ ( t , γ ) ) y 1 ¯ ( t , γ ) = y 2 ̲ ( t , γ ) , y 2 ¯ ( t , γ ) = f ̲ ( t , y 1 ̲ ( t , γ ) , y 1 ¯ ( t , γ ) , y 2 ̲ ( t , γ ) , y 2 ¯ ( t , γ ) ) y 1 ̲ ( 0 , γ ) = y 1 ̲ 0 ( γ ) , y 1 ¯ ( 0 , γ ) = y 1 ¯ 0 ( γ ) y 2 ̲ ( 0 , γ ) = y 2 ¯ 0 ( γ ) , y 2 ¯ ( 0 , γ ) = y 2 ̲ 0 ( γ )
This system is used to initiate the fuzzy solution of the differential equations, and switching points are determined by using the fact that f has differentiability type-(i) when f . f > 0 ; otherwise, it has type-(ii). This shows that there is only one possible solution at one point, and the other solution is possible at another point. In between these two points, the definition of differentiability changes from type-(i) to type-(ii) and vice versa. Similarly, if f . f > 0 , then f has type-(i) differentiability [41]. The point where f . f > 0 changes to f . f < 0 is called the switching point. This concept is similarly applied to second derivatives.
Definition 11.
(Fuzzy Power Series ([45])). Let u and u 0 be two numbers of the fuzzy number space L ( R ) , and their difference u u 0 is a nonnegative fuzzy number. Let ( v n ) be a sequence of nonnegative fuzzy numbers. Now, define [ u u 0 ] γ = [ w ] γ = [ w ̲ γ , w ¯ γ ] with w ̲ γ = h ̲ ( γ ) u 0 ̲ ( γ ) and w ¯ γ = h ¯ ( γ ) u 0 ¯ ( γ ) and v γ = [ g ̲ ( γ ) , g ¯ ( γ ) ] . Then, the power series of fuzzy numbers with the coefficients v n is given by n = 1 v n ( u u 0 ) n , which can be expressed in terms of γ-level sets as follows:
n = 1 v n [ w ] γ n = [ n = 1 g ̲ n ( γ ) [ w ̲ ( γ ) ] n , n = 1 g ¯ n ( γ ) [ w ¯ ( γ ) ] n ]
Definition 12.
The differential transform of the n t h -order derivative of a function f ( x ) is defined as
F ( n ) = 1 n ! d n ( f ( x ) ) d x n x = x 0
Then, F ( n ) , the inverse transform of the above relation, is defined as
f ( x ) = n = 0 F ( n ) ( x x 0 ) n
Equation (9) represents the Taylor series expansion of function f ( x ) at x = x 0 .

4. Differential Transformation Method

First, the function and its differential transformation is provided in Table 3.
These basic functions are used for developing the differential transformation solution of the problem.

4.1. Extension of Differential Transformation Method to Fuzzy Differentiability

First, definitions of fuzzy differentiability of type-(i) and type-(ii) are given.
Definition 13.
([35]) If X ( t , r ) is strongly differentiable of order n in the domain T, then
  • fis (i)-differentiable:
    φ ¯ ( t , n , r ) = d n ( x ¯ ) ( t , r ) d t n t T X i ¯ ( N , r ) = φ ¯ ( t , n , r ) = d n ( x ¯ ) ( t , r ) d t n ] t = t i , n N φ ̲ ( t , n , r ) = d n ( x ̲ ) ( t , r ) d t n , t T X i ̲ ( N , r ) = φ ̲ ( t , n , r ) = d n ( x ̲ ) ( t , r ) d t n ] t = t i , n N
  • fis (ii)-differentiable
    φ ¯ ( t , n , r ) = d n ( x ̲ ) ( t , r ) d t n t T X i ¯ ( N , r ) = φ ̲ ( t , n , r ) = d n ( x ¯ ) ( t , r ) d t n ] t = t i , i f   n   i s   o d d φ ̲ ( t , k , r ) = d k ( x ¯ ) ( t , r ) d t n t T X i ̲ ( N , r ) = φ ¯ ( t , n , r ) = d n ( x ¯ ) ( t , r ) d t n ] t = t i , i f   n   i s   e v e n
    where X ¯ ( N , r ) and X ̲ ( N , r ) are the upper and lower spectrum of x ( t , r ) at t = t i in domain N. When the whole interval of validity is divided into subintervals, the switching points are verified at the end points of the subintervals. If there is a switching point, then the definition of fuzzy differentiability changes from type-(i) to type-(ii). If f is first differentiable, then x ( t , r ) is represented as
    x ¯ ( t , n , r ) = n = 0 ( t t i ) n n ! X ( n , r ) ¯ x ̲ ( t , n , r ) = n = 0 ( t t i ) n n ! X ( n , r ) ̲
The change in definition is handled with the application of switching points. So, only this definition is followed and the above equations are inverse transformation of X ( n ) , if X ( n ) is defined as
X i ̲ ( N , r ) = M ( n ) d n ( q ( t ) x ( t , r ) ̲ ) d t n ] t = t i , n = 0 , 1 , 2 , X i ¯ ( N , r ) = M ( n ) d n ( q ( t ) x ( t , r ) ¯ ) d t n ] t = t i , n = 0 , 1 , 2 ,
Then the function x ( t , r ) can be rewritten as
x ¯ ( t , n , r ) = 1 q ( t ) n = 0 ( t t i ) n n ! X ( N , r ) ¯ M ( n ) x ̲ ( t , n , r ) = 1 q ( t ) n = 0 ( t t i ) n n ! X ( N , r ) ̲ M ( n )
where M ( n ) and q ( t ) are positive numbers and M ( n ) is a weighting function and q ( t ) is a kernel function. For simplicity, put q ( t ) = 1 and M ( n ) = H n n ! and H is a time horizon. Thus the above equations can be written as
x ¯ ( t , n , r ) = n = 0 ( t t i H ) n X ( N , r ) ¯ x ̲ ( t , n , r ) = n = 0 ( t t i H ) n X ( N , r ) ̲
The differential transformation method is a numerical technique based on the Taylor series.
Definition 14.
The function is transformed by differential transform as,
[ Y ( n ) ] γ = 1 n ! [ d n ( y ( x ) ̲ γ ) d x n , d n ( y ( x ) ¯ γ ) d x n ] x = x 0
where Y ( n ) is a fuzzy function, and inverse of differential transform for the function [ Y ( n ) ] γ is defined as
[ y ( n ) ] γ = [ n = 0 Y γ ( n ) ̲ ( x x 0 ̲ ) γ n , n = 0 Y γ ( n ) ¯ ( x x 0 ¯ ) γ n ]
where Equation (11) is known as the fuzzy Taylor series expansion of fuzzy function y ( x ) at fuzzy point x = x 0 .
Let us consider the second-order initial value problem
y ( t ) = f ( t , y ( t ) , y ( t ) ) , y ( 0 ) = y 0 y ( 0 ) = y 0
This equation is converted into an equal problem of system of first-order differential equations, written as
y 1 ( t ) = y 2 ( t ) , y 1 ( 0 ) = y ( 0 ) y 2 ( t ) = f ( t , y 1 ( t ) , y 2 ( t ) ) , y 2 ( 0 ) = y ( 0 )
using the substitution y 1 ( t ) = y ( t ) and y 2 ( t ) = y ( t ) . This system of initial value problem is fuzzified and problem if solved using the differential transformation.
( k + 1 ) Y 1 ̲ ( k + 1 , r ) = Y 2 ̲ ( k , r ) ( k + 1 ) Y 1 ¯ ( k + 1 , r ) = Y 2 ¯ ( k , r ) ( k + 1 ) Y 2 ̲ ( k + 1 , r ) = [ F ( t ̲ , t ¯ , Y 1 ̲ ( k ) , Y 1 ¯ ( k , r ) , Y 2 ̲ ( k , r ) , Y 2 ¯ ( k , r ) ) ( k + 1 ) Y 2 ̲ ( k + 1 , r ) = [ F ( t ̲ , t ¯ , Y 1 ̲ ( k ) , Y 1 ¯ ( k , r ) , Y 2 ̲ ( k , r ) , Y 2 ¯ ( k , r ) )
where F denotes the transformed function. Now the whole domain is divided into N sub-domains and the switching points are verified at these points. From the initial conditions, we get Y 1 ̲ ( 0 , r ) = y 1 ̲ 0 ( r ) , Y 1 ¯ ( 0 , r ) = y 1 ¯ 0 ( r ) and Y 2 ̲ ( 0 , r ) = y 2 ̲ 0 ( r ) , Y 2 ¯ ( 0 , r ) = y 2 ¯ 0 ( r ) for the first sub-domain, the Taylor series can be written as
y 1 ̲ 0 ( t , r ) = Y 1 ̲ 0 ( 0 , r ) + Y 1 ̲ 0 ( 1 , r ) ( t a ) + Y 1 ̲ 0 ( 1 , r ) ( t a ) 2 + + Y 1 ̲ 0 ( n , r ) ( t a ) n y 1 ¯ 0 ( t , r ) = Y 1 ¯ 0 ( 0 , r ) + Y 1 ¯ 0 ( 1 , r ) ( t a ) + Y 1 ¯ 0 ( 1 , r ) ( t a ) 2 + + Y 1 ¯ 0 ( n , r ) ( t a ) n y 2 ̲ 0 ( t , r ) = Y 2 ̲ 0 ( 0 , r ) + Y 2 ̲ 0 ( 1 , r ) ( t a ) + Y 2 ̲ 0 ( 1 , r ) ( t a ) 2 + + Y 2 ̲ 0 ( n , r ) ( t a ) n y 2 ¯ 0 ( t , r ) = Y 2 ¯ 0 ( 0 , r ) + Y 2 ¯ 0 ( 1 , r ) ( t a ) + Y 2 ¯ 0 ( 1 , r ) ( t a ) 2 + + Y 2 ¯ 0 ( n , r ) ( t a ) n
In these equations, subscript 0 shows that the series is defined at the initial point. This series can be evaluated at t = t 1 , then we define
y 1 ̲ 0 ( t 1 , r ) = Y 1 ̲ 0 ( 0 , r ) + Y 1 ̲ 0 ( 1 , r ) ( t 1 a ) + Y 1 ̲ 0 ( 1 , r ) ( t 1 a ) 2 + + Y 1 ̲ 0 ( n , r ) ( t 1 a ) n = Y 1 ̲ 0 ( 0 , r ) + Y 1 ̲ 0 ( 1 , r ) h + Y 1 ̲ 0 ( 1 , r ) h 2 + + Y 1 ̲ 0 ( n , r ) h n = j = 0 n Y 1 ̲ 0 ( j , r ) h j y 1 ¯ 0 ( t 1 , r ) = Y 1 ¯ 0 ( 0 , r ) + Y 1 ¯ 0 ( 1 , r ) ( t 1 a ) + Y 1 ¯ 0 ( 1 , r ) ( t 1 a ) 2 + + Y 1 ¯ 0 ( n , r ) ( t 1 a ) n = Y 1 ¯ 0 ( 0 , r ) + Y 1 ¯ 0 ( 1 , r ) h + Y 1 ¯ 0 ( 1 , r ) h 2 + + Y 1 ¯ 0 ( n , r ) h n = j = 0 n Y 1 ¯ 0 ( j , r ) h j y 2 ̲ 0 ( t 1 , r ) = Y 2 ̲ 0 ( 0 , r ) + Y 2 ̲ 0 ( 1 , r ) ( t 1 a ) + Y 2 ̲ 0 ( 1 , r ) ( t 1 a ) 2 + + Y 2 ̲ 0 ( n , r ) ( t 1 a ) n = Y 2 ̲ 0 ( 0 , r ) + Y 2 ̲ 0 ( 1 , r ) h + Y 2 ̲ 0 ( 1 , r ) h 2 + + Y 2 ̲ 0 ( n , r ) h n = j = 0 n Y 2 ̲ 0 ( j , r ) h j y 2 ¯ 0 ( t 1 , r ) = Y 2 ¯ 0 ( 0 , r ) + Y 2 ¯ 0 ( 1 , r ) ( t 1 a ) + Y 2 ¯ 0 ( 1 , r ) ( t 1 a ) 2 + + Y 2 ¯ 0 ( n , r ) ( t 1 a ) n = Y 2 ¯ 0 ( 0 , r ) + Y 2 ¯ 0 ( 1 , r ) h + Y 2 ¯ 0 ( 1 , r ) h 2 + + Y 2 ¯ 0 ( n , r ) h n = j = 0 n Y 2 ¯ 0 ( j , r ) h j
The final value of y ( t 1 , r ) is the initial value for the next sub-domain. At this point switching is implement and verified by applying the min and max operators. Thus y 1 ̲ ( t 1 , r ) = min ( y 1 ̲ ( t 1 , r ) , y 1 ¯ ( t 1 , r ) ) and y 1 ¯ ( t 1 , r ) = max ( y 1 ̲ ( t 1 , r ) , y 1 ¯ ( t 1 , r ) ) . If y 1 ̲ ( t 1 , r ) = y 1 ̲ ( t 1 , r ) then there is no switch otherwise there is a switch and definition of differentiability changes from type-(i) to type-(ii) and vice versa. This is also verified for y 2 . As a result, the initial conditions for the series at t 1 are recomputed from the series of the previous subdomain.
y 1 ̲ 1 ( t 2 , r ) = Y 1 ̲ 1 ( 0 , r ) + Y 1 ̲ 1 ( 1 , r ) ( t 2 t 1 ) + Y 1 ̲ 1 ( 1 , r ) ( t 2 t 1 ) 2 + + Y 1 ̲ 1 ( n , r ) ( t 1 a ) n = Y 1 ̲ 1 ( 0 , r ) + Y 1 ̲ 1 ( 1 , r ) h + Y 1 ̲ 1 ( 1 , r ) h 2 + + Y 1 ̲ 1 ( n , r ) h n = j = 0 n Y 1 ̲ 1 ( j , r ) h j y 1 ¯ 1 ( t 2 , r ) = Y 1 ¯ 1 ( 0 , r ) + Y 1 ¯ 1 ( 1 , r ) ( t 2 t 1 ) + Y 1 ¯ 1 ( 1 , r ) ( t 2 t 1 ) 2 + + Y 1 ¯ 1 ( n , r ) ( t 1 a ) n = Y 1 ¯ 1 ( 0 , r ) + Y 1 ¯ 1 ( 1 , r ) h + Y 1 ¯ 1 ( 1 , r ) h 2 + + Y 1 ¯ 1 ( n , r ) h n = j = 0 n Y 1 ¯ 1 ( j , r ) h j y 2 ̲ 1 ( t 2 , r ) = Y 2 ̲ 1 ( 0 , r ) + Y 2 ̲ 1 ( 1 , r ) ( t 2 t 1 ) + Y 2 ̲ 1 ( 1 , r ) ( t 2 t 1 ) 2 + + Y 2 ̲ 1 ( n , r ) ( t 2 t 1 ) n = Y 2 ̲ 1 ( 0 , r ) + Y 2 ̲ 1 ( 1 , r ) h + Y 2 ̲ 1 ( 1 , r ) h 2 + + Y 2 ̲ 1 ( n , r ) h n = j = 0 n Y 2 ̲ 0 ( j , r ) h j y 2 ¯ 1 ( t 1 , r ) = Y 2 ¯ 1 ( 0 , r ) + Y 2 ¯ 1 ( 1 , r ) ( t 2 t 1 ) + Y 2 ¯ 1 ( 1 , r ) ( t 2 t 1 ) 2 + + Y 2 ¯ 1 ( n , r ) ( t 2 t 1 ) n = Y 2 ¯ 1 ( 0 , r ) + Y 2 ¯ 1 ( 1 , r ) h + Y 2 ¯ 1 ( 1 , r ) h 2 + + Y 2 ¯ 1 ( n , r ) h n = j = 0 n Y 2 ¯ 1 ( j , r ) h j
Continuing the same process and series at the point t i + 1 will be
y 1 ̲ i + 1 ( t i , r ) = Y 1 ̲ i ( 0 , r ) + Y 1 ̲ i ( 1 , r ) ( t i + 1 t i ) + Y 1 ̲ i ( 1 , r ) ( t i + 1 t i ) 2 + + Y 1 ̲ i ( n , r ) ( t i + 1 t i ) n = Y 1 ̲ i ( 0 , r ) + Y 1 ̲ i ( 1 , r ) h + Y 1 ̲ i ( 1 , r ) h 2 + + Y 1 ̲ i ( n , r ) h n = j = 0 n Y 1 ̲ i ( j , r ) h j y 1 ¯ i + 1 ( t i , r ) = Y 1 ¯ i ( 0 , r ) + Y 1 ¯ i ( 1 , r ) ( t i + 1 t i ) + Y 1 ¯ i ( 1 , r ) ( t i + 1 t i ) 2 + + Y 1 ¯ i ( n , r ) ( t i + 1 t i ) n = Y 1 ¯ i ( 0 , r ) + Y 1 ¯ i ( 1 , r ) h + Y 1 ¯ i ( 1 , r ) h 2 + + Y 1 ¯ i ( n , r ) h n = j = 0 n Y 1 ¯ i ( j , r ) h j y 2 ̲ i + 1 ( t i , r ) = Y 2 ̲ i ( 0 , r ) + Y 2 ̲ i ( 1 , r ) ( t i + 1 t i ) + Y 2 ̲ i ( 1 , r ) ( t i + 1 t i ) 2 + + Y 2 ̲ i ( n , r ) ( t i + 1 t i ) n = Y 2 ̲ i ( 0 , r ) + Y 2 ̲ i ( 1 , r ) h + Y 2 ̲ i ( 1 , r ) h 2 + + Y 2 ̲ i ( n , r ) h n = j = 0 n Y 2 ̲ i ( j , r ) h j y 2 ¯ i + 1 ( t i , r ) = Y 2 ¯ i ( 0 , r ) + Y 2 ¯ i ( 1 , r ) ( t i + 1 t i ) + Y 2 ¯ i ( 1 , r ) ( t i + 1 t i ) 2 + + Y 2 ¯ i ( n , r ) ( t i + 1 t i ) n = Y 2 ¯ i ( 0 , r ) + Y 2 ¯ i ( 1 , r ) h + Y 2 ¯ i ( 1 , r ) h 2 + + Y 2 ¯ i ( n , r ) h n = j = 0 n Y 2 ¯ i ( j , r ) h j

5. Numerical Examples

In this section, we consider the numerical examples of first and second order.
Example 1.
Consider the following FIVP:
y ( t ) + y ( t ) + y ( t ) = t y ( 0 ) = 0 , y ( 0 ) = 1
Solution:
The analytical solution of this problem is y ( t ) = exp ( t 2 ) ( cos ( 3 t 2 + 1 3 sin ( 3 t 2 ) + t 1 . Before fuzzifying the differential equation, we first analyze the switching points of the solution function as mentioned in Section 3. For this purpose, the graph of the solution function and its first and second derivatives in Figure 1.
The graphs show that f ( t ) . f ( t ) > 0 and f ( t ) . f ( t ) < 0 for t [ 0 , 1.204 ] . This shows that f has differentiability of type-(i), f has differentiability of type-(ii), and the second-order differential equation has the solution of the D 1 , 2 2 system. This solution is switched for t ( 1.204 , 5 ] , where f ( t ) . f ( t ) > 0 and f ( t ) . f ( t ) > 0 , and the equation has the solution of type D 1 , 1 2 within the interval. This shows that the solution switches from type D 1 , 2 2 to D 1 , 1 1 .
Now, the second-order differential equation is converted into a system of first-order differential equations. Differential Equation (13) is converted into a system of first-order differential equations, represented in Equation (14).
y 1 ( t ) = y 2 ( t ) y 2 ( t ) = t y 2 ( t ) y 1 ( t )
The initial conditions also change to y 1 ( 0 ) = 0 , y 2 ( 0 ) = 1 . Replacing t * = t t i , then above equation will become
y 1 ( t * ) = y 2 ( t * ) y 2 ( t * ) = t * + t i y 2 ( t * ) y 1 ( t * )
Every real number can be represented as the triangular fuzzy number. The initial values are fuzzified, and this problem becomes an FIVP. These values are represented as the interval-valued fuzzy set
y 1 ( 0 , γ ) = [ y 1 ̲ , y 1 ¯ ] γ = [ 1 + γ , 1 γ ]
Similarly, the second value is fuzzified as
y 2 ( 0 , γ ) = [ y 2 ̲ , y 2 ¯ ] γ = [ 1 + γ , 1 γ ]
The system of ODEs is converted into the equivalent system of ordinary differential equations with fuzzy initial value conditions. This process splits the system of ODEs into two equivalent systems of ODEs, and Equation (15) will become
y 1 ̲ ( t * , γ ) = y 2 ̲ ( t * , γ ) y 1 ¯ ( t * , γ ) = y 2 ¯ ( t * , γ ) y 2 ̲ ( t * , γ ) = t * + t i y 2 ¯ ( ( t * , γ ) ) y 1 ¯ ( t * , γ ) y 2 ¯ ( t * , γ ) = t * + t i y 2 ̲ ( t * , γ ) y 1 ̲ ( t * , γ )
Taking the differential transformation of the above system of equations, we get the recurrence relation for this system of equations between the two grid points t i and t i + 1 , represented as
Y 1 ̲ i ( k + 1 , γ ) = Y 2 ̲ i ( k , γ ) / ( k + 1 ) Y 1 ¯ i ( k + 1 , γ ) = Y 2 ¯ i ( k , γ ) / ( k + 1 ) Y 2 ̲ i ( k + 1 , γ ) = δ ( k 1 ) + t i Y 2 ¯ i ( ( k , γ ) ) Y 1 ¯ i ( k , γ ) / ( k + 1 ) Y 2 ¯ i ( k + 1 , γ ) = δ ( k 1 ) + t i Y 2 ̲ i ( k , γ ) Y 1 ̲ i ( k , γ ) / ( k + 1 )
with initial conditions, we get Y 1 ̲ 0 ( 0 , γ ) = 1 + γ ; Y 1 ¯ 0 ( 0 , γ ) = 1 γ and Y 2 ̲ 0 ( 0 , γ ) = γ Y 2 ¯ 0 ( 0 , γ ) = 2 γ . The whole interval is divided into the subintervals, and switching points are checked at each partition point. The differential transformation coefficients are determined from the above set of equations, and the lower and upper solutions are developed for each subdomain separately using the differential transformation method of order 20, namely, D T 20 . The final value of the preceding subdomain is used as the initial value of the next subdomain. For example, at t = 0.2 , y(t) in the first subinterval will be y ( t ) γ = [ 1.1144 γ 1.0204 , 1.2084 1.1144 γ ] such that γ [ 0 , 1 ] . The results are compared in Table 4, in which the computed solutions for the grid point 0 are shown for the changing fuzzification parameter γ.
The point t = 1.2 is the switching point, where y 2 ̲ = 0.7724 , y 2 ¯ = 0.7551 . For the next input, the values are interchanged, i.e., the upper solution becomes the lower solution and the lower solution changes to the upper solution by applying the min and max operators.
Example 2.
Consider the following FIVP:
y ( t ) 2 y ( t ) 8 y ( t ) = 8 t 2 + 20 t + 2 t > 0 y ( 0 ) = 0 y ( 0 ) = 8
over the interval [ 0 , 3 ] .
Solution:
The analytical solution of the problem is y ( t ) = 5 3 exp ( 2 t ) + 5 3 exp ( 4 t ) t 2 2 t .
In the first step, the domain of the analytic solution function is analyzed, and the Figure 2 of the solution function with its first and second derivatives is drawn. This graph shows that f ( t ) . f ( t ) > 0 for all t [ 0 , ] . As marked in Section 3, if f has a solution according to type-(i) derivatives, and if f ( t ) . f ( t ) > 0 for all t [ 0 , ] as a result, then f also has a derivative of type-(i). Thus, the problem has the D 1 , 1 2 solution in the D m , n system.
The solution of a second-order differential equation is determined by converting Differential Equation (18) into a system of first-order ODEs according to Equation (19).
y 1 ( t ) = y 2 ( t ) y 2 ( t ) = 2 y 2 ( t ) + 8 y 1 ( t ) + 8 t 2 + 20 t + 2
The initial conditions are also transferred according to the used transformation, y 1 ( 0 ) = 0 , y 2 ( 0 ) = 8 . Now, replacing t * = t t i ,
y 1 ( t * ) = y 2 ( t * ) y 2 ( t * ) = 2 y 2 ( t * ) + 8 y 1 ( t * ) + 8 ( t * ) 2 + 8 t i 2 + 20 t i + 2 + ( 16 t i + 20 ) ( t * )
The differential transformation of this ODE system at grid point t i is computed as Equation (21).
Y 1 ̲ i ( k + 1 , γ ) = Y 2 ̲ i ( k , γ ) / ( k + 1 ) Y 1 ¯ i ( k + 1 , γ ) = Y 2 ¯ i ( k , γ ) / ( k + 1 ) Y 2 ̲ i ( k + 1 , γ ) = 2 [ Y 2 ̲ i ( k , γ ) + 8 Y 1 ̲ i ( k , γ ) + 8 δ ( k 2 ) + ( 20 + 16 t i ) δ ( k 1 ) + ( 8 t i 2 + 20 t i + 2 ) δ ( k ) ] / ( k + 1 ) Y 2 ¯ i ( k + 1 , γ ) = [ 2 Y 2 ¯ i ( k , γ ) + 8 Y 1 ̲ i ( k , γ ) + 8 δ ( k 2 ) + ( 20 + 16 t i ) δ ( k 1 ) + ( 8 t i 2 + 20 t i + 2 ) δ ( k ) ] / ( k + 1 )
To solve the second-order initial value problem, the system of initial value problem is converted to the fuzzy initial value problem. For this, the initial conditions are fuzzified as y 1 ( 0 ) = ( 1 + γ , 1 γ ) and y 2 ( 0 ) = ( 7 + γ , 9 γ ) , and γ-level transformation of these initial conditions will be Y 1 ̲ ( 0 , γ ) = 1 + γ , Y 1 ¯ ( 0 , γ ) = 1 γ and Y 2 ̲ ( 0 , γ ) = 7 + γ Y 2 ¯ ( 0 , γ ) = 7 + γ . Similarly, the recurrence relation for the upper solution will be written as stated in Equation (21). The computed solution for multiple grid points is shown in Table 5 for γ = 1 .
The results show that there is no switching point, and the definition of the derivative is chosen according to the initial conditions.
Example 3.
Consider the following FIVP:
y ( t ) 6 y 7 y ( t ) = 9 e 2 t 0 t 2 y ( 0 ) = 2 , y ( 0 ) = 13
Solution:
y ( t ) = e t 2 e 7 t e 2 t is the analytical solution of Differential Equation (22).
The Figure 3 shows that the function and its first and second derivatives are always negative. This results in f ( t ) . f ( t ) > 0 for all t [ 0 , ] . As a result, f has a solution according to the type-(i) derivative, and f ( t ) . f ( t ) > 0 for all t [ 0 , ] ; as a result, f has a solution according to the type-(i) derivative, and the problem has the D 1 , 1 2 solution in the D m , n system.
For a numerical solution of a higher-order differential equation with fuzzy initial values, first convert Equation (22) in the system of first-order differential equations as Equation (23).
y 1 = y 2 y 2 = 6 y 2 + 7 y 1 9 e 2 t
The initial conditions become y 1 ( 0 ) = 2 , y 2 ( 0 ) = 13 . The system of initial value problem is converted to the fuzzy initial value problem, and these conditions become y 1 ( 0 , γ ) = ( 3 + γ , 1 γ ) and y 2 ( 0 , γ ) = ( 7 + γ , 9 γ ) . The differential equation for the system between the two grid points t i and t i + 1 can be represented as Equation (24), the recurrence relation for the differential transform method:
Y 1 ( K + 1 ) = Y 2 ( K ) Y 2 ( K + 1 ) = 6 Y 2 ( k ) + 7 Y 1 ( k ) 9 ( 2 ) k k ! e 2 t i
Now, for the fuzzy initial value solutions, the lower system of the recurrence relation can be written as Equation (25).
Y 1 ̲ i ( k + 1 , γ ) = Y 2 ̲ i ( k , γ ) / ( k + 1 ) Y 1 ̲ i ( k + 1 , γ ) = Y 2 ̲ i ( k , γ ) / ( k + 1 ) Y 2 ̲ i ( k + 1 , γ ) = [ 6 Y 2 ̲ i ( k , γ ) + 7 Y 1 ̲ i ( k , γ ) 9 ( 2 ) k k ! e 2 t i ] / ( k + 1 ) Y 2 ̲ i ( k + 1 , γ ) = [ 6 Y 2 ̲ i ( k , γ ) + 7 Y 1 ̲ i ( k , γ ) 9 ( 2 ) k k ! e 2 t i / ( k + 1 )
The transformation of initial conditions yields Y 1 ̲ 0 ( 0 , γ ) = 1 + γ , Y 1 ¯ 0 ( 0 , γ ) = 1 γ and Y 2 ̲ 0 ( 0 , γ ) = 7 + γ , Y 2 ¯ 0 ( 0 , γ ) = 9 γ . The series solution is computed for the upper and lower approximation in each subdomain separately. The initial values at the point t = 0.9 are the terminal values of the previous series, which is
y 1 ̲ = 2.1518 γ + 2.6452 , y 1 ¯ = 6.9488 2.1518 γ y 2 ̲ = 6.9108 γ + 19.1886 , y 2 ¯ = 33.0103 6.91086 γ
A similar method is followed for other points, and the results are represented in Table 6.
These results show that the solution has no switching point, and the system has D 1 , 1 2 differentiability in the domain.

6. Discussion

Second-order differential equations with fuzzy initial value problems are solved using a system of first-order differential equations. The second-order differential equation has four possible solutions in the D n , m 2 system. In this paper, we discuss the switching points, which are the points where the definition of differentiability switches from type-(i) to type-(ii) or vice versa. These sorts of solutions are developed for second-order differential equations. It is shown that only one definition of differentiability is valid for one point or subinterval. This definition switches to the second definition at the switching point. Three numerical examples are solved using the min, max operators for switching points, and promising results validate the method. These solutions are important for the development of automatic solutions to modeling and simulation of problems.

7. Conclusions

In this paper, we develop the solution of a second-order fuzzy initial value problem with the analysis of switching points using the differential transformation method. The definition of fuzzy differentiability switches from type-(i) to type-(ii) and vice versa when the product of the function and its derivative changes signs within the domain of the function. These switching points are handled by applying the min, max operators. This demonstrates the single solution of an FIVP at time t or for the interval. This method is used to solve second-order differential equations by converting second-order ODEs into an equivalent system of first-order ODEs. Application of the min, max operators will be applied for solving higher-order differential equations with fuzzy initial values in future works.

Author Contributions

The author, N.S., contributed mainly in the preparation and solution of the differential equations. The solution process was discussed a number of times with the second author, M.M. He also provided the simulation of results with M.M.S.M., M.M.S.M. provided the technical help regarding the code development and write-up of manuscript.

Funding

This article has no funding resources.

Acknowledgments

The authors are thankful to the reviewers for their valued time and constructive input which helped a lot to improve the quality of the manuscript.

Conflicts of Interest

The authors have no conflict of interest.

References

  1. Zadeh, L. Fuzzy sets. Inf. Control 1965, 8, 338–353. [Google Scholar] [CrossRef]
  2. Chang, S.S.L.; Zadeh, L. On fuzzy mapping and control. IEEE Trans. Syst. Man Cybern. 1972, 2, 30–34. [Google Scholar] [CrossRef]
  3. Dubois, D.; Prade, H. Towards fuzzy differential calculus: Part 3. differentiation. Fuzzy Sets Syst. 1982, 8, 225–233. [Google Scholar]
  4. Geotschel, R.; Voxman, W. Elementary fuzzy calculus. Fuzzy Sets Syst. 1986, 18, 31–43. [Google Scholar]
  5. Puri, M.L.; Ralescu, D.A. Differentials of fuzzy functions. J. Math. Anal. Appl. 1983, 91, 552–558. [Google Scholar] [CrossRef]
  6. Fard, O.S.; Ghal-Eh, N. Numerical solutions for linear system of first-order fuzzy differential equations with fuzzy constant coefficients. Inf. Sci. 2011, 181, 4765–4779. [Google Scholar] [CrossRef]
  7. Qiu, D.; Zhang, W. Symmetric fuzzy numbers and additive equivalence of fuzzy numbers. Soft Comput. 2013, 17, 1471–1477. [Google Scholar] [CrossRef]
  8. Qiu, D.; Lu, C.; Zhang, W.; Lan, Y. Algebraic properties and topological properties of the quotient space of fuzzy numbers based on Mares equivalence relation. Fuzzy Sets Syst. 2014, 245, 63–82. [Google Scholar] [CrossRef]
  9. Qiu, D.; Zhang, W.; Lu, C. On fuzzy differential equations in the quotient space of fuzzy numbers. Fuzzy Sets Syst. 2016, 295, 72–98. [Google Scholar] [CrossRef]
  10. Ahmadi, M.B.; Kiania, N.A. Differential transformation method for solving fuzzy differential inclusions by fuzzy partitions. Int. J. Ind. Math. 2013, 5, 237–249. [Google Scholar]
  11. Plotnikov, A.V.; Skripnik, N.V. Fuzzy Differential Equations with Generalized Derivative. J. Fuzzy Set Valued Anal. 2012, 2012, jfsva-00113. [Google Scholar] [CrossRef]
  12. Bede, B.; Gal, S.G. Solutions of fuzzy differential equations based on generalized differentiability. Commun. Math. Anal. 2010, 9, 22–41. [Google Scholar]
  13. Buckley, J.J.; Feuring, T. Fuzzy differential equations. Fuzzy Sets Syst. 2000, 110, 43–54. [Google Scholar]
  14. Georgiou, D.; Nieto, J.J.; Rodriguez-Lopez, R. Initial value problems for higher-order fuzzy differential equations. Nonlinear Anal. Theory Methods Appl. 2005, 63, 587–600. [Google Scholar]
  15. Bede, B.; Gal, S.G. Generalizations of the Differentiability of Fuzzy-number-valued Functions with Applications to Fuzzy Differential Equations. Fuzzy Sets Syst. 2005, 151, 581–599. [Google Scholar]
  16. Effati, S.; Pakdaman, M. Artificial Neural Network Approach for Solving Fuzzy Differential Equations. Inf. Sci. 2010, 180, 1434–1457. [Google Scholar]
  17. Jafari, H.; Saeidy, M.; Baleanu, D. The variational iteration method for solving n-th order fuzzy differential equations. Central Eur. J. Phys. 2012, 10, 76–85. [Google Scholar]
  18. Babolian, E.; Sadeghi, H.; Javadi, S. Numerically Solution of Fuzzy Differential Equations by Adomian Method. Appl. Math. Comput. 2004, 149, 547–557. [Google Scholar]
  19. Allahviranloo, T.; Ahmady, N.; Ahmady, E. Numerical solution of fuzzy differential equations by predictor-corrector method. Inf. Sci. 2007, 177, 1633–1647. [Google Scholar]
  20. Allahviranloo, T.; Ahmady, E.; Ahmady, N. A method for solving nth order fuzzy linear differential equations. Int. J. Comput. Math. 2009, 86, 730–742. [Google Scholar]
  21. Allahviranloo, T.; Shafiee, M.; Nejatbakhsh, Y. A Note on “Fuzzy Differential Equations and the Extension Principle”. Inf. Sci. 2009, 179, 2049–2051. [Google Scholar]
  22. Ma, M.; Friedman, M.; Kandel, A. Numerical solutions of fuzzy differential equations. Fuzzy Sets Syst. 1999, 105, 133–138. [Google Scholar]
  23. Allahviranloo, T.; Abbasbandy, S.; Salahshour, S.; Hakimzadeh, A. A New Method for Solving Fuzzy Linear Differential Equations. Computing 2011, 92, 181–197. [Google Scholar]
  24. Allahviranloo, T.; Salahshour, S. Euler Method for Solving Hybrid Fuzzy Differential Equation. Soft Comput. 2011, 15, 1247–1253. [Google Scholar]
  25. de Barros, L.C.; Pedro, F.S. Fuzzy differential equations with interactive derivative. Fuzzy Sets Syst. 2017, 309, 64–80. [Google Scholar] [CrossRef]
  26. Isa, S.; Abdul Majid, Z.; Ismail, F.; Rabiei, F. Diagonally Implicit Multistep Block Method of Order Four for Solving Fuzzy Differential Equations Using Seikkala Derivatives. Symmetry 2018, 10, 42. [Google Scholar] [CrossRef]
  27. Ahmadian, A.; Salahshour, S.; Chan, C.S.; Baleanu, D. Numerical solutions of fuzzy differential equations by an efficient Runge-Kutta method with generalized differentiability. Fuzzy Sets Syst. 2018, 331, 47–67. [Google Scholar] [CrossRef]
  28. Allahviranloo, T.; Chehlabi, M. Solving Fuzzy Differential Equations Based on the Length Function Properties. Soft Comput. 2015, 19, 307–320. [Google Scholar]
  29. Long, H.V.; Nieto, J.J.; Son, N.T.K. New approach for studying nonlocal problems related to differential systems and partial differential equations in generalized fuzzy metric spaces. Fuzzy Sets Syst. 2018, 331, 26–46. [Google Scholar] [CrossRef]
  30. Zhang, D.; Feng, W.; Zhao, Y.; Qiu, J. Global existence of solutions for fuzzy second-order differential equations under generalized H-differentiability. Comput. Math. Appl. 2010, 60, 1548–1556. [Google Scholar] [CrossRef]
  31. Rabiei, F.; Ismail, F.; Ahmadian, A.; Salahshour, S. Numerical Solution of Second-Order Fuzzy Differential Equation Using Improved Runge-Kutta Nystrom Method. Math. Probl. Eng. 2013, 2013, 803462. [Google Scholar]
  32. Najariyan, M.; Mazandarani, M. A note on “Numerical solutions for linear system of first-order fuzzy differential equations with fuzzy constant coefficients”. Inf. Sci. 2015, 305, 93–96. [Google Scholar] [CrossRef]
  33. Zhou, X. Differential Transformation and its Applications for Electric Circuits; Hauzhong University Press: Wuhan, China, 1986; Volume 24, pp. 1279–1289. [Google Scholar]
  34. Mirzaee, F. Differential Transformation Method for Solving the Linear and Nonlinear Systems of Ordinary Differential Equations. Appl. Math. Sci. 2011, 5, 3465–3472. [Google Scholar]
  35. Allahviranloo, T.; Kiani, N.; Motamedi, N. Solving fuzzy differential equations by differential transformation method. Inf. Sci. 2009, 179, 956–966. [Google Scholar]
  36. Akin, O.; Khaniyev, T.; Oruç, O.; Türksen, I. An algorithm for the solution of second order fuzzy initial value problems. Expert Syst. Appl. 2013, 40, 953–957. [Google Scholar]
  37. Guo, X.; Shang, D.; Lu, X. Fuzzy approximate solutions of second-order fuzzy linear boundary value problems. Bound. Value Probl. 2013, 2013, 212. [Google Scholar] [CrossRef]
  38. Kadkhoda, N.; Roushan, S.S.; Jafari, H. Differential Transform Method: A Tool for Solving Fuzzy Differential Equations. Int. J. Appl. Comput. Math. 2017, 4, 33. [Google Scholar] [CrossRef]
  39. Khastan, A.; Nieto, J.J. A boundary value problem for second order fuzzy differential equations. Nonlinear Anal. Theory Methods Appl. 2010, 72, 3583–3593. [Google Scholar] [CrossRef]
  40. Kadak, U.; Ozluk, M. A new approach for the sequence spaces of fuzzy level sets with the partial metric. J. Fuzzy Set Valued Anal. 2014, 2014, 1–13. [Google Scholar] [CrossRef]
  41. Khastan, A.; Bahrami, F.; Ivaz, K. New Results on Multiple Solutions for nth-Order Fuzzy Differential Equations under Generalized Differentiability. Bound. Value Probl. 2009, 2009, 395714. [Google Scholar] [CrossRef]
  42. Bede, B.; Stefanini, L. Generalized differentiability of fuzzy-valued functions. Fuzzy Sets Syst. 2013, 230, 119–141. [Google Scholar]
  43. Mansouri, S.S.; Ahmad, N. A Numerical Method For Solving Nth-Order Fuzzy Differential Equation by using Characterization Theorem. Commun. Numerical Anal. 2012, 2012, cna-00054. [Google Scholar] [CrossRef]
  44. Allahviranloo, T.; Kiani, N.; Barkhordari, M. Toward the existence and uniqueness of solutions of second-order fuzzy differential equations. Inf. Sci. 2009, 179, 1207–1215. [Google Scholar]
  45. Kadak, U.; Basar, F. Power series of fuzzy numbers with real or fuzzy coefficients. Filomat 2012, 26, 519–528. [Google Scholar]
Figure 1. Graph of the function and its first and second derivatives. The black line shows the function, the dashed line shows the first derivative, and the blue line shows the second derivative.
Figure 1. Graph of the function and its first and second derivatives. The black line shows the function, the dashed line shows the first derivative, and the blue line shows the second derivative.
Mathematics 07 00231 g001
Figure 2. Graph of the function and its first and second derivatives. The black line shows the function, the dashed line shows the first derivative, and the blue line shows the second derivative.
Figure 2. Graph of the function and its first and second derivatives. The black line shows the function, the dashed line shows the first derivative, and the blue line shows the second derivative.
Mathematics 07 00231 g002
Figure 3. Graph of the function and its first and second derivatives. The black line shows the function, the dashed line shows the first derivative, and the blue line shows the second derivative.
Figure 3. Graph of the function and its first and second derivatives. The black line shows the function, the dashed line shows the first derivative, and the blue line shows the second derivative.
Mathematics 07 00231 g003
Table 1. Nature of the function and type of differentiability; ↗ and ↘ stand for increasing and decreasing functions, respectively.
Table 1. Nature of the function and type of differentiability; ↗ and ↘ stand for increasing and decreasing functions, respectively.
S. N.f f Relationship with Derivatives of gDifferentiability
1 f < g type-(i)
2 f > g type-(ii)
Table 2. Nature of second-order derivatives and their relationship with derivatives of g.
Table 2. Nature of second-order derivatives and their relationship with derivatives of g.
S. N.f f f ( 2 ) Relationship with Derivatives of gType of Differentiability
1 f < g and f ( 2 ) < g ( 2 ) f and f have differentiability type-(I), and the system D 1 , 1 ( 2 ) is also type-(I)
2 f < g and f ( 2 ) > g ( 2 ) f follows type-(I), f has type-(II), and D 1 , 2 ( 2 ) is also type-(II)
3 f > g and f ( 2 ) > g ( 2 ) f follows type-(II), f has type-(I), and the system is represented by D 2 , 1 ( 2 )
4 f > g and f ( 2 ) > g ( 2 ) f and f have differentiability type-(II), and the system D 2 , 2 ( 2 ) is also type-(II)
Table 3. The functions and their differential transformation.
Table 3. The functions and their differential transformation.
Function FormDifferential Transform
f ( t ) = g ( t ) ± h ( x ) F ( j ) = G ( j ) ± G ( j )
f ( t ) = a h ( t ) F ( j ) = a H ( j )
f ( t ) = h ( t ) t F ( j ) = ( j + 1 ) H ( j + 1 )
f ( t ) = r h ( t ) t r F ( j ) = ( j + 1 ) ( j + 2 ) ( j + r ) H ( j + r )
F ( j ) = ( j + r ) ! j ! H ( j + r )
f ( t ) = h ( t ) . g ( t ) F ( j ) = r = 0 j G ( r ) H ( j r )
f ( t ) = t m F ( j ) = δ ( t m )
f ( t ) = e x p ( λ t ) F ( j ) = λ j j !
f ( t ) = ( 1 + t ) m F ( j ) = m ( m 1 ) ( m 2 ) ( m j + 1 ) j !
f ( t ) = sin ( ω t + γ ) F ( j ) = ω j j ! sin ( π j 2 ! + γ )
f ( t ) = cos ( ω t + γ ) F ( j ) = ω j j ! cos ( π j 2 ! + γ )
Table 4. Comparison of solutions for multiple t and γ = 1 .
Table 4. Comparison of solutions for multiple t and γ = 1 .
γ 0.20.40.60.811.21.41.61.82
y ̲ 0.09400.18070.26220.34020.41670.65820.56870.56940.64900.7331
y ¯ 0.09400.18070.26220.34030.41710.49230.65820.73460.81560.9028
Table 5. Computation of solution at multiple grid points for γ = 1 .
Table 5. Computation of solution at multiple grid points for γ = 1 .
γ 0.30.60.91.21.51.82.12.4
y ̲ 1.66774.797011.024223.486748.232896.9402192.1757334.9335
y ¯ 1.66774.797011.024323.487548.2340293.8502720.95831336
Table 6. Computation of lower and upper approximate solutions at multiple t for γ = 1 .
Table 6. Computation of lower and upper approximate solutions at multiple t for γ = 1 .
t00.40.81.21.62
y ̲ −2−27.0747−425.7578−6953.7e+03−1.1424e+05−1.9238e+06
y ¯ −2−27.0747−425.7578−6953.7e+03−1.1424e+05−1.8331e+06

Share and Cite

MDPI and ACS Style

Salamat, N.; Mustahsan, M.; Saad Missen, M.M. Switching Point Solution of Second-Order Fuzzy Differential Equations Using Differential Transformation Method. Mathematics 2019, 7, 231. https://doi.org/10.3390/math7030231

AMA Style

Salamat N, Mustahsan M, Saad Missen MM. Switching Point Solution of Second-Order Fuzzy Differential Equations Using Differential Transformation Method. Mathematics. 2019; 7(3):231. https://doi.org/10.3390/math7030231

Chicago/Turabian Style

Salamat, Nadeem, Muhammad Mustahsan, and Malik M. Saad Missen. 2019. "Switching Point Solution of Second-Order Fuzzy Differential Equations Using Differential Transformation Method" Mathematics 7, no. 3: 231. https://doi.org/10.3390/math7030231

APA Style

Salamat, N., Mustahsan, M., & Saad Missen, M. M. (2019). Switching Point Solution of Second-Order Fuzzy Differential Equations Using Differential Transformation Method. Mathematics, 7(3), 231. https://doi.org/10.3390/math7030231

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop