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

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.


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.
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.

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.

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) > 0for 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: and . The subtraction of fuzzy numbers h − g is defined as addition, and the parameterized subtraction is defined as The scalar multiplication of a scalar a with a fuzzy number is computed as follows: ah(γ) = ah(γ), ah(γ) = ah(γ) when a ≥ 0 and ah(γ) = ah(γ), ah(γ) = ah(γ) when a < 0, and the scalar multiplication will become ah γ = [ah(γ), ah(γ)] Multiplication: Let h, g be two fuzzy numbers; then, multiplication of fuzzy numbers is computed as follows: hg = [hg, hg], where (hg Division: Division between two fuzzy numbers is defined conditionally. g(γ)] and 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 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 Definition 5 (Hausdorff Distance).Let F : R f × R f → R + ∪ {0} be a function.Then, the Hausdorff distance is defined as 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: 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 for all h < 0 If h is sufficiently small and h < 0, then 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 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.
If f g (x 0 ) exists, then f is called differentiable in the generalized sense.
where ( f (x) γ and g(x) γ are differentiable real-valued functions with respect to x.Then, F(x) is g_ differentiable, and mathematically, 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].) If f , g are differentiable, then we can write F (2) in a similar fashion: 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.
] be an interval-valued function, where f is an increasing function and g is a decreasing function, and let the system D (1) 2 f : I → R f be two fuzzy functions.Then, 1.
Let f and g be two differentiable functions.If F and D (1)

2.
Let f and g be two differentiable functions.If F is (I)-differentiable and D (1)

3.
Let f and g be two differentiable functions.If F is (II)-differentiable and D (1) Let f and g be two differentiable functions.If F and D (1) Proof.For details, see [41].
Theorem 4 (Theorem 3.1 in [44]).Let t 0 ∈ [a, b] and assume that 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: , where x and x are (i)-differentials, s))ds)ds, where x and x are (ii)-differentials, 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 gH differentiability, so it is also gH differentiable.

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.

S. N. f f
Relationship with Derivatives of g Differentiability 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.
Table 2. Nature of second-order derivatives and their relationship with derivatives of g.

S. N.
f f f (2) Relationship with Derivatives of g Type of Differentiability 1 f < g and f (2) < g (2)  f and f have differentiability type-(I), and the system D 1,1 is also type-(I) 2 f < g and f (2) > g (2)  f follows type-(I), f has type-(II), and D 1,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 f > g and f (2) > g (2)  f and f have differentiability type-(II), and the system 2,2 is also type-(II) 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.
This second-order differential equation is decomposed into the system of first-order differential equation as where transformation y1(t) = y(t) and y2(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: 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: Definition 12.The differential transform of the n th -order derivative of a function f (x) is defined as Then, F(n), the inverse transform of the above relation, is defined as Equation ( 9) represents the Taylor series expansion of function f (x) at x = x 0 .

Differential Transformation Method
First, the function and its differential transformation is provided in Table 3.
Table 3.The functions and their differential transformation.

Function Form Differential Transform
These basic functions are used for developing the differential transformation solution of the problem.

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 f is (i)-differentiable: 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 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 Then the function x(t, r) can be rewritten as 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 The differential transformation method is a numerical technique based on the Taylor series.
Definition 14.The function is transformed by differential transform as, where Y(n) is a fuzzy function, and inverse of differential transform for the function [Y(n)] γ is defined as where Equation ( 11) is known as the fuzzy Taylor series expansion of fuzzy function y(x) at fuzzy point x = x 0 .
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 Y1(0, r) = y1 0 (r), Y1(0, r) = y1 0 (r) and Y2(0, r) = y2 0 (r), Y2(0, r) = y2 0 (r) for the first sub-domain, the Taylor series can be written as 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 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 y1(t 1 , r) = min(y1(t 1 , r), y1(t 1 , r)) and y1(t 1 , r) = max(y1(t 1 , r), y1(t 1 , r)).If y1(t 1 , r) = y1(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 y2.As a result, the initial conditions for the series at t 1 are recomputed from the series of the previous subdomain.
Continuing the same process and series at the point t i+1 will be

Numerical Examples
In this section, we consider the numerical examples of first and second order.

Solution:
The analytical solution of this problem is y(t) = exp( −t 2 )(cos( 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 2 1,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 2  1,1 within the interval.This shows that the solution switches from type D 2  1,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).
The initial conditions also change to y1(0) = 0, y2(0) = 1.Replacing t * = t − t i , then above equation will become 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 Similarly, the second value is fuzzified as 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 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 with initial conditions, we get Y1 0 (0, γ) = −1 + γ; Y1 0 (0, γ) = 1 − γ and Y2 0 (0, γ) = γ Y2 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, DT20.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 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 y2 = 0.7724, y2 = 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.

Solution:
The analytical solution of the problem is y(t) = −5 3 exp(−2t) + 5 3 exp(4t) − t 2 − 2t.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 2  1,1 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).
The initial conditions are also transferred according to the used transformation, y The differential transformation of this ODE system at grid point t i is computed as Equation ( 21).
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 y1(0) = (−1 + γ, 1 − γ) and y2(0) = (7 + γ, 9 − γ), and γ-level transformation of these initial conditions will be Y1(0, γ) = −1 + γ, Y1(0, γ) = 1 − γ and Y2(0, γ) = 7 + γ Y2(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: Solution: 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 2  1,1 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).
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: Now, for the fuzzy initial value solutions, the lower system of the recurrence relation can be written as Equation (26).

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 2 n,m 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.

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.

Figure 1 .
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 2 .
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.

5 - 8 × 10 9 - 6 × 10 9 - 4 × 10 9 - 2 × 10 9 Figure 3 .
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.

Table 1 .
Nature of the function and type of differentiability;

Table 4 .
Comparison of solutions for multiple t and γ = 1.

Table 5 .
Computation of solution at multiple grid points for γ = 1.

Table 6 .
Computation of lower and upper approximate solutions at multiple t for γ = 1.