International Journal of Mathematical Engineering and Science
ISSN : 2277-6982 http://www.ijmes.com/
Volume 1 Issue 5(May 2012) https://sites.google.com/site/ijmesjournal/
Numerical Solution of N th -Order Fuzzy Differential Equation by Runge-Kutta Nystrom Method T. Jayakumar, K. Kanagarajan, S. Indrakumar Department of Mathematics Sri Ramakrishna Mission Vidyalaya College of Arts and science Coimbatore-641 020, Tamilnadu, India.
[email protected],
[email protected] [email protected] Abstract : In this paper we study a numerical solution for N th -order fuzzy differential equations based on Seikkala derivative with initial value problem. The Runge-Kutta Nystrom method is used for the numerical solution of this problem and the convergence and stability of the method is proved. By this method, we can obtain strong fuzzy solution. This method is illustrated by solving examples. Keywords : Fuzzy Numbers, N th -order FDEs, RK-Nystrom method.
1 Introduction The topic of fuzzy differential equations (FDEs) have been rapidly growing in recent years. The concept of fuzzy derivative was first introduced by Chang and Zadeh [10], it was followed up by Dubois and Prade [12] by using the extension principle in their approach. Other methods have been discussed by Puri and Ralescu [26] and Goetschel and Voxman [15]. Kandel and Byatt [25] applied the concept of fuzzy differential equation (FDE) to the analysis of fuzzy dynamical problems. The FDE and the initial value problem (Cauhy problem) were rigorously treated by Kaleva [21, 22], Seikkala [28], He and Yi [16], Kloeden [23]and by other researchers (see [6, 20]). The numerical methods for solving fuzzy differential equations are introduced in [1, 2, 3]. Buckley and Feuring [9] introduced two analytical methods for solving N th -order linear differential equations with fuzzy initial value conditions. Their first method of solution was to fuzzy if the crisp solution and then check to see if it satiesfies the differential equation with fuzzy initial value conditions; and the second method was the reverse of the first method, they first solved the fuzzy initial value condition and the checked to see if it defined a fuzzy function. In this paper, a numerical method to solve N th -order linear differential equations with fuzzy initial conditions is presented. The structure of the paper is organized as follows: In Section 2, we give some basic results on fuzzy numbers and define a fuzzy derivative and a fuzzy integral. then the fuzzy initial values is treated in Section 3 using the extension principle of Zadeh and the concept of fuzzy derivative. It is shown that the fuzzy initial value problem has a unique fuzzy solution when f satiesfies Lipschitz condition which guarantees a unique solution to the deterministic initial value problem. In Section 4, the Runge-Kutta 1
International Journal of Mathematical Engineering and Science
ISSN : 2277-6982 http://www.ijmes.com/
Volume 1 Issue 5(May 2012) https://sites.google.com/site/ijmesjournal/
Nystrom method for solving N th -order fuzzy differential equations is introduced. In Section 5 Convergence and Stability are illustrated. In Section 6 the proposed method is illustrated by solving several examples, and the conclusion is drawn in section 7.
2 Preliminaries An arbitrary fuzzy number is represented by an ordered pair of functions (u(r), u(r)) for all r ∈ [0, 1], which satisfy the following requirements [12]: (i) u(r) is a bounded left continuous non-decreasing function over [0, 1], (ii) u(r) is a bounded right continuous non-increasing function over [0, 1], (iii) u(r) ≤ u(r) ∀ r ∈ [0, 1], let E be the set of all upper semi-continuous normal convex fuzzy numbers with bounded α-level intervals.
Lemma 2.1 Let [v(α), v(α)], α ∈ (0, 1] be a given family of non-empty intervals. If (i) [v(α), v(α)] ⊃ [v(β), v(β)] for 0 < α ≤ β, and (ii) [ lim v(αk ), lim v(αk )] = [v(α), v(α)], k→∞
k→∞
whenever (αk ) is a non-decreasing sequence converging to α ∈ (0, 1], then the family [v(α), v(α)], α ∈ (0, 1], represent the α-level sets of a fuzzy number v ∈ E. Conversely if [v(α), v(α)], α ∈ (0, 1], are α-level sets of a fuzzy number v ∈ E, then the conditions (i) and (ii) hold true.
Definition 2.1 Let I be a real interval. A mapping v : I → E is called a fuzzy process and we denoted the α-level set by [v(t)]α = [v(t, α), v(t, α)]. The Seikkala derivative v ′ (t) of v is defined by ′
′
[v(t)]α = [v (t, α), v (t, α)], provided that is a equation defines a fuzzy number v ′ (t) ∈ E.
Definition 2.2 Suppose u and v are fuzzy sets in E. Then their Hausdroff D : E × E → R+ ∪ {0}, n o D(u, v) = supα∈[0,1] max |u(α) − v(α)|, |u(α) − v(α)| ,
i.e., D(u, v) is maximal distance between α-level sets of u and v.
2
International Journal of Mathematical Engineering and Science
ISSN : 2277-6982 http://www.ijmes.com/
Volume 1 Issue 5(May 2012) https://sites.google.com/site/ijmesjournal/
3 Fuzzy Initial Value Problem Now we consider the initial value problem (n) x (t) = ψ(t, x, x′ , . . . , x(n−1) ), x(0) = a1 , . . . , x(n−1) (0) = an ,
(1)
where ψ is a continuous mapping from R+ ×Rn into R and ai (0 ≤ i ≤ n) are fuzzy numbers in E. The mentioned nth order fuzzy differential equation by changining variables y1 (t) = x(t), y2 (t) = x′ (t), . . . , yn (t) = x(n−1) (t), converts to the following fuzzy system ′ y1 (t) = f1 (t, y1 , . . . , yn ), .. . ′ y (t) = fn (t, y1 , . . . , yn ), n [0] [0] y1 (0) = y1 = a1 , . . . , yn (0) = yn = an ,
(2)
[0]
are fuzzy
y ′i (t, α) = min{fi (t, u1 , . . . , un ); uj ∈ [y j (t, α), y j (t, α)]} = f i (t, y(t, α)),
(3)
y ′i (t, α) = max{fi (t, u1 , . . . , un ); uj ∈ [y j (t, α), y j (t, α)]} = f i (t, y(t, α)),
(4)
where fi (1 ≤ i ≤ n) are continuous mapping from R+ × Rn into R and yi numbers in E with α-level intervals [0]
[0]
(α), y i (α)] f or i = 1, . . . , n, and 0 < α ≤ 1. [yi ]α = [y [0] i We call y = (y1 , . . . , yn )t is a fuzzy solution of (2) on an interval I, if
and y i (0, α) = y [0] (α), i
[0]
y i (0, α) = y i (α).
Thus for fixed α we have a system of initial value problem in R2n . If we can solve it (uniquely), we have only to verify that the intervals, [y j (t, α), y j (t, α)] define a fuzzy number [0]
[0]
yi (t) ∈ E. Now let y [0] (α) = (y [0] (α), . . . , y [0] (α))t and y [0] (α) = (y 1 (α), . . . , y n (α))t , with 1 n respect to the above mentioned indicators, system (2) can be written as with assumption ′ y (t) = F (t, y(t)), (5) y(0) = y [0] ∈ E n . With assumption y(t, α) = [y(t, α), y(t, α)] and y ′ (t, α) = [y ′ (t, α), y ′ (t, α)] where y(t, α) = (y(t, α), . . . , y(t, α))t , y(t, α) = (y(t, α), . . . , y(t, α))t ,
(6) (7)
y ′ (t, α) = (y ′ (t, α), . . . , y ′ (t, α))t , y ′ (t, α) = (y ′ (t, α), . . . , y ′ (t, α))t ,
(8) (9)
3
International Journal of Mathematical Engineering and Science
ISSN : 2277-6982 http://www.ijmes.com/
Volume 1 Issue 5(May 2012) https://sites.google.com/site/ijmesjournal/
and with assumption F (t, y(t, α)) = [F (t, y(t, α)), F (t, y(t, α))], where F (t, y(t, α)) = (f 1 (t, y(t, α)), . . . , f n (t, y(t, α))t , t
F (t, y(t, α)) = (f 1 (t, y(t, α)), . . . , f n (t, y(t, α)) ,
(10) (11)
y(t) is a fuzzy solution of (5) on an interval I for all α ∈ (0, 1], if ′ y (t, α) = F (t, y(t, α)); y ′ (t, α) = F (t, y(t, α)); y(0, α) = y [0] (α), y(0, α) = y [0] (α).
or
′ y (t, α) = F (t, y(t, α)),
(12)
(13)
y(0, α) = y [0] (α).
Now we show that under the assumption for functions fi , for i = 1, . . . , n how we can obtain a unique fuzzy solution for system (2).
Theorem 3.1 If fi (t, u1 , . . . , un ) for i = 1, . . . , n are continuous function of t and satiesfies the Lipschitz condition in u = (u1 , . . . , un )t in the region D = {(t, u)| t ∈ [0, 1], −∞ < ui < ∞ for i = 1, . . . , n} with constant Li then the initial value problem (2) has a unique fuzzy solution in each case. Proof. Denote G = (F , F )t = (f 1 , . . . , f n , f 1 , . . . , f n )t where f i (t, u) = min{fi (t, u1 , . . . , un ); uj ∈ [y j , y j ], f or j = 1, . . . , n},
(14)
f i (t, u) = max{fi (t, u1 , . . . , un ); uj ∈ [y j , y j ], f or j = 1, . . . , n},
(15)
and y = (y, y)t = (y 1 , . . . , y n , y 1 , . . . , y n )t ∈ R2n . It can be shown that Lipschitz condition of functions fi imply ||F (t, z) − F (t, z ∗ )|| ≤ L||z − z ∗ ||. This guarantees the existence and uniqueness solution of y ′ (t) = F (t, y(t)), y(0) = y [0] = (y [0] , y [0] )2t ∈ R2n .
(16)
Also for any continuous function y [1] : R+ → R2n the successive approximations y
[m+1]
(t) = y
[0]
+
Z
t
F (s, y [m] (s))ds, t ≥ 0, m = 1, 2, . . . , 0
4
(17)
International Journal of Mathematical Engineering and Science
ISSN : 2277-6982 http://www.ijmes.com/
Volume 1 Issue 5(May 2012) https://sites.google.com/site/ijmesjournal/
converge uniformly on closed subintervals of R+ to the solution of (16). In other word we have the following successive approximations Z t [0] f i (s, y [m] (s))ds, f or i = 1, 2, . . . , n. (t) = y + y [m+1] (18) i i 0
[m+1]
yi
[0]
(t) = y i +
Z
t 0
f i (s, y [m] (s))ds, f or i = 1, 2, . . . , n.
(19)
By choosing y [0] = (y [0] (α), y [0] (α)) in (16) we get a unique solution y α (t) = (y(t, α), y(t, α)) to (3) and (4) for each α ∈ (0, 1]. Next we will show that the y(t, α) = [y(t, α), y(t, α)], defines a fuzzy number in E n for each 0 ≤ t ≤ T , i.e. that y = (y1 , . . . , yn )t is a fuzzy solution to (14) and (15). Thus we will show that the intervals [y i (t, α), y i (t, α)], for i = 1, . . . , n satisfy the conditions of Lemma (2.1). The successive approximations y [m] = y [0] ∈ E n , Z t [m+1] [0] F (s, y [m] (s))ds, t ≥ 0, m = 1, 2, . . . , (20) y (t) = y + 0
where the integrals are the fuzzy integrals, define a sequence of fuzzy numbers [m] [m] [m] [m] y [m] (t) = (y1 (t), . . . , yn (t))t for each 0 ≤ t ≤ T . Hence [yi (t)]β ⊂ [yi (t)]α , if 0 < α ≤ β ≤ 1, which implies that [m]
[m]
(t, β), y i (t, β)] ⊂ [y [m] (t, α), y i (t, α)], 0 < α ≤ β ≤ 1, [y [m] i i [m]
since, by convergence of sequences (16) and (19), the end points of [yi (t)]α converge to y i (t, α) and y i (t, α) that means [m]
y [m] (t, α) → y i (t, α) and i
y i (t, α) → y i (t, α).
(21)
Thus the inclusion property (i) of Lemma (2.1) holds for the intervals [y i (t, α), y i (t, α)], for 0 < α ≤ 1 . For the proof of the property (ii) of Lemma (2.1), let (αp ) be a non-decreasing sequence in (0,1] converging to α. Then y [0] (αp ) → y [0] (α) and y [0] (αp ) → y [0] (α), because of y [0] ∈ E n . But by the continuous dependence on the initial value of the solution (16), y(t, αp ) → y(t, α) and y(αp ) → y(α), this means (ii) holds for the interval [y(t, α), y(t, α)], for 0 < α ≤ 1. Hence by Lemma (2.1), y(t) ∈ E n and so y is a fuzzy solution of (1). the uniqueness follows from the uniqueness of the solution of (16).
4 The Runge-Kutta Nystrom Method With before assumptions, the initial values problem (2) has a unique solution, such as y = (y1 , . . . , yn )t ∈ E n . for found an approximate solution for (2) with the Runge-Kutta Nystrom method, we first define y(tn+1 ; α) − y(tn ; α) =
3 P
wi ki,1 (tn , y(tn ; α)),
i=1
y(tn+1 ; α) − y 2 (tn ; α) =
3 P
i=1
5
(22) wi ki,2 (tn , y(tn ; α)),
International Journal of Mathematical Engineering and Science
ISSN : 2277-6982 http://www.ijmes.com/
Volume 1 Issue 5(May 2012) https://sites.google.com/site/ijmesjournal/
[ki (t, y(t; α))]α = [ki,1 (t, y(t; α)), ki,2 (t, y(t; α))], i = 1, 2, 3 ki,1 (tn , y(tn ; α)) = h.f (tn + ai h, y1 (tn ) +
i−1 P
bij kj,1 (tn , y(tn ; α))),
j=1
ki,2 (tn , y(tn ; α)) = h.f (tn + ai h, y2 (tn ) +
i−1 P
(23)
bij kj,2 (tn , y(tn ; α))),
j=1
and k1,1 (t, y(t; α)) = min{h.f (t, s1 , . . . , sn )|si ∈ [y1 (t; α), y2 (t; α)]}, (i = 1, 2, . . . , n) k1,2 (t, y(t; α)) = max{h.f (t, s1 , . . . , sn )|si ∈ [y1 (t; α), y2 (t; α)]}, k2,1 (t, y(t; α)) = min{h.f (t + 23 h, s1 , . . . , sn )|si ∈ [z1,1 (t, y(t; α), h), z1,2 (t, y(t; α), h)]},
(24)
k2,2 (t, y(t; α)) = max{h.f (t + 32 h, s1 , . . . , sn )|si ∈ [z1,1 (t, y(t; α), h), z1,2 (t, y(t; α), h)]}, k3,1 (t, y(t; α)) = min{h.f (t + 23 h, s1 , . . . , sn )|si ∈ [z2,1 (t, y(t; α), h), z2,2 (t, y(t; α), h)]}, k3,2 (t, y(t; α)) = max{h.f (t + 32 h, s1 , . . . , sn )|si ∈ [z2,1 (t, y(t; α), h), z2,2 (t, y(t; α), h)]}, where in the Runge-Kutta Nystrom method, z1,1 (t, y(t; α), h) = y1 (t; α) + 23 k1,1 (t, y(t; α)), z1,2 (t, y(t; α), h) = y2 (t; α) + 32 k1,2 (t, y(t; α)),
(25)
z2,1 (t, y(t; α), h) = y1 (t; α) + 32 k2,1 (t, y(t; α), h), z2,2 (t, y(t; α), h) = y2 (t; α) + 23 k2,2 (t, y(t; α), h). Define, F (t, y(t; α), h) = 2k1,1 (t, y(t; α), h) + 3k2,1 (t, y(t; α), h) + 3k3,1 (t, y(t; α), h), (26) G(t, y(t; α), h) = 2k1,2 (t, y(t; α), h) + 3k2,2 (t, y(t; α), h) + 3k3,2 (t, y(t; α), h). and suppose that the discrete equally spaced grid points {t0 = 0, t1 , . . . , tN = T } is a partition for interval [0, T ]. If the exact and approximate solution in the i-th α cut at [m]
[m]
[m]
(α), y i (α)] and [Y i (α), Y i (α)] respectively, then tm , 0 ≤ m ≤ N are denoted by [y [m] i the numerical method for solution approximation in the i-th coordinate α cut, with the Runge-Kutta Nystrom method is [m+1]
Yi Y
[m]
(α) = Y i (α) + h8 Fi (tm , Y m (α), h),
[m+1] (α) i
=Y
[m] i (α)
+
h m 8 Gi (tm , Y (α), h),
6
[0]
(α), Y i (α) = y [0] i Y
[0] i (α)
=
[0] y i (α),
(27)
International Journal of Mathematical Engineering and Science
ISSN : 2277-6982 http://www.ijmes.com/
Volume 1 Issue 5(May 2012) https://sites.google.com/site/ijmesjournal/
where [Yi (t)]α = [Y i (t, α), Y i (t, α)], Y [m] (α) = [Y [m] (α), Y [m]
t Y [m] (α) = (Y 1 (α), . . . , Y [m] n (α)) , and Y
[m]
[m]
(α)]
[m]
[m]
(α) = (Y 1 (α), . . . , Y n (α))t .
(28)
Now we input F ∗ (t, Y [m] (α), h) = 18 (F1 (t, Y [m] (α), h), . . . , Fn (t, Y [m] (α), h))t ,
(29)
G∗ (t, Y [m] (α), h) = 18 (G1 (t, Y [m] (α), h), . . . , Gn (t, Y [m] (α), h))t . The Runge-Kutta Nystrom method for solutions approximation α-cut of differential equation (13) is as follows Y [m+1] (α) = Y [m] (α) + hH(tm , Y [m] (α), h),
Y [0] (α) = y [0] (α)
(30)
where H(tm , Y [m] (α), h) = [F ∗ (tm , Y [m] (α), h), G∗ (tm , Y [m] (α), h)],
(31)
and F ∗ (tm , Y [m] (α), h) =
1 [2k (tm , Y [m] (α), h) + 3k 2 (tm , Y [m] (α), h) 8 1 +3k 3 (tm , Y [m] (α), h)],
1 G∗ (tm , Y [m] (α), h) = [2k 1 (tm , Y [m] (α), h) + 3k 2 (tm , Y [m] (α), h) 8
(32)
+3k 3 (tm , Y [m] (α), h)] and also k j (t, Y [m] (α), h) = (k 1j (t, Y [m] (α), h), . . . , k nj (t, Y [m] (α), h))t , (33) k j (t, Y [m] (α), h) = (k 1j (t, Y [m] (α), h), . . . , k nj (t, Y [m] (α), h))t .
5 Convergence and Stability Definition 5.1 A one-step method for approximating the solution of a differential equation ′ y (t) = F (t, y(t)), y(0) = y [0] ∈ Rn ,
(34)
which F is a n-th ordered as follows f = (f1 , . . . , fn )t and fi : R+ × Rn → R(1 ≤ i ≤ n), is a method which can be written in the form Y [n+1] = Y [n] + hψ(tn , Y [n] , h), 7
(35)
International Journal of Mathematical Engineering and Science
ISSN : 2277-6982 http://www.ijmes.com/
Volume 1 Issue 5(May 2012) https://sites.google.com/site/ijmesjournal/
where the increment function ψ is determined by F and is a function tn , Y [n] and h only.
Theorem 5.1 If ψ(t, y, h) satisfies a Lipschitz condition in y, then the method given by (35) is stable.
Theorem 5.2 In relation (5), if F (t, y) satisfies a Lipschitz condition in y, then the method given by (30) is stable.
Theorem 5.3 If Y [m+1] (α) = Y [m] (α) + hψ(tm , Y [m] (α), h),
Y [0] = y [m] (α)
(36)
where ψ(tm , Y [m] (α), h) = [ψ1 (tm , Y [m] (α), h), ψ2 (tm , Y [m] (α), h)] is a numerical method for approximation of differential equation (13), and ψ1 and ψ2 are continuous in t, y, h for 0 ≤ t ≤ T ,0 ≤ h ≤ h0 and all y, and if they satisfy a Lipschitz condition in the region D = {(t, u, v, h)| 0 ≤ t ≤ T, −∞ < ui ≤ vi , −∞ < vi ≤ +∞, 0 ≤ h ≤ h0 i = 1, . . . , n}, necessary and sufficient conditions for convergence above mentioned method is ψ(t, y(t, α), h) = F (t, y(t, α)).
(37)
Proof. Suppose that ψ(t, y(t, α), 0) = F (t, y(t, α)), since, F (t, y(t, α)) satisfying the conditions of theory (1), then the following equation ′ y (t) = F (t, y(t)), (38) y(0) = y [0] (α),
has a unique solution such as y(t, α) = [y(t, α), y(t, α)], where y(t, α) = (y(t, α), . . . , y(t, α))t and y(t, α) = (y(t, α), . . . , y(t, α))t . We will show that the numerical solutions given by (36) convergent to the y(t). By the mean value theorem, = y [m] + hf i (tm + θi h, y(tm + θi h)), y [m+1] i i [m+1]
yi
[m]
= yi
+ hf i (tm + θi h, y(tm + θi h)), t
f or 0 < θi < 1,
(39)
f or 0 < θi < 1,
(40)
t
with assumption ψ = (ψ 1 , . . . , ψ n ) and ψ = (ψ 1 , . . . , ψ n ) . From equation (36) obtain the following relations [m+1]
Yi
[m+1]
Yi
[m]
(α) = Y i (α) + hψ i (tm , Y [m] (α), h), [m]
(α) = Y i (α) + hψ i (tm , Y [m] (α), h),
and subtracting (39), (40) from (41), (42) respectively, and setting e[m] (α) = [e[m] (α), e[m] (α)], 8
(41) (42)
International Journal of Mathematical Engineering and Science
ISSN : 2277-6982 http://www.ijmes.com/
Volume 1 Issue 5(May 2012) https://sites.google.com/site/ijmesjournal/
where e[m] (α) = e(tm , α) = Y [m] (α) − y [m] (α) and e[m] (α) = e(tm , α) = Y
[m]
(α) − y [m] (α),
we get [m+1]
ei
(α)
[m]
ei (α) + h{ψ i (tm , Y [m] (α), h) − ψ i (tm , y [m] (α), h) + ψ i (tm , y [m] (α), h)
=
−ψ i (tm , y [m] (α), 0) + ψ i (tm , y [m] (α), 0) − f i (tm + θi h, y(tm + θi h)} on the other way, with respect to the relation of ψ i (tm , y [m] (α), 0) − f i (tm , y [m] (α)) we can write |ψ i (tm , y [m] (α), 0) − f i (tm + θi h, y(tm + θi h))| n n P P ≤ hL1 θi + L1 |y i (tm + θi h) − y i (tm )| |y i (tm + θi h) − y i (tm )| + L1 = hL1 θi + L1
then [m+1] |ei (α)|
≤
[m] |ei (α)
≤
[m] |ei (α)
i=1 n P i=1
|y ′i (tm + ξ i θi h)θi h| + L1
+ hL1
X n
[m] |ej (α)|
+
1≤j≤n
n
i=1
[m] ej (α)|
j=1
+ nhL1 max
n P
i=1
|y ′i (tm + ξ i θi h)θi h| = hL2 ,
+ h 2 L1 + h 2 L2
n o o [m] [m] |ej (α)| + nhL1 max |ej (α)| + h2 (L1 + L2 ).
On the other hand [m] [m] max |ej (α)| = ki |ei (α)|, 1≤j≤n
1≤j≤n
max
1≤j≤n
[m] [m] |ej (α)| = ki′ |ei (α)|
with assumption k1 = max ki , ki′ and M = L1 + L2 , we can write 1≤i≤n
[m+1]
|ei
(α)|
[m] [m] [m] ≤ |ei (α)| + nhk1 L1 |ei (α)| + |ei (α)| + M1 h2 [m] [m] [m] ≤ |ei (α)| + 2nhk1 L1 max |ei (α)|, |ei (α)| + M1 h2 ,
(43)
similarly, we can obtain the following relation [m+1]
|ei
[m]
[m]
[m]
(α)| ≤ |ei (α)| + 2nhk2 L′1 max{|ei (α)|, |ei (α)|} + M2 h2 .
(44)
Now, we input L = max{L1 , L′1 } and M = max{M1 , M2 } so the relations (43) and (44) can be written as follow [m] [m+1] [m] [m] |ei (α)| ≤ |ei (α)| + 2nhkL max |ei (α)|, |ei (α)| + M h2 , [m] [m+1] [m] [m] (α)| ≤ |ei (α)| + 2nhkL max |ei (α)|, |ei (α)| + M h2 |ei [m]
Denote ei
[m]
[m]
= |ei (α)| + |ei (α)|. Then by virture of lemma (5.7) [m]
[0]
ei (α) ≤ (1 + 4nhkL)m ei (α) + 2M h2 9
(1 + 4nhkL)m − 1 , 4nhkL
International Journal of Mathematical Engineering and Science
ISSN : 2277-6982 http://www.ijmes.com/
[0]
[0]
Volume 1 Issue 5(May 2012) https://sites.google.com/site/ijmesjournal/
[0]
where ei = |ei (α)| + |ei (α)|. Then [m]
|ei (α)| [m]
|ei (α)|
e4mnklh h 2nhkL e4mnklh [0] ≤ e4mnkh × ei + M h. 2nhkL [0]
≤ e4mnkh × ei + M
In particular [N ]
[0]
|ei (α)| ≤ e4N nkh × ei + M [0]
e4N nklh e4N nklh [N ] [0] h and |ei (α)| ≤ e4N nkh × ei + M h. 2nhkL 2nhkL
[0]
Since ei (α) = ei (α) = 0 and h =
T N
[N ]
we obtain ||e1 (α)|| ≤ M
e4N nklh h and 2nhkL
e4N nklh e4N nklh h. Then ||e[N ] (α)|| ≤ 2M h if h → 0 we get ||e[N ] (α)|| → 0, 2nhkL 2nhkL so the numerical solution (36) converge to the solutions (38). Conversely, suppose that the numerical method (36) convergent to the solution of the system (38). With absurd hypothesis we suppose that (37) is not correct. Then ψ(t, y(t, α), 0) = g(t, y(t, α)) 6= F (t, y(t, α)). Similarly, we can proof that the numerical method of (36) is convergent to the solution of following system ′ u (t) = g(t, y(t)), (45) u(0) = y [0] (α), [N ]
||e2 (α)|| ≤ M
then y(t, α) = u(t). Since g(t, y(t, α)) 6= F (t, y(t, α)), suppose that F and g differ at some point (tα , y(tα , α)). If we consider the initial value problem (38) and (45) starting from (tα , y(tα , α)) we have y ′ (tα , α) = F (tα , y(tα , α)) 6= g(tα , y(tα , α)) = g(tα , u(tα )) = u′ (tα ) which is a contradiction.
Corollary 5.1 The Runge-Kutta proposed method by (30) and is convergent to the solution of the system (13) respectively.
6 Numerical Example Example 6.1 Consider the following fuzzy differential equation with fuzzy initial value ′′′ y (t) = 2y ′′ (t) + 3y ′ (t) y(0) = (3 + α, 5 − α), y ′ (0) = (−1 − α, −3 + α), ′′ y (0) = (8 + α, 10 − α).
(0 ≤ t ≤ 1),
the eigen value -eigen vector solution is as follows: 1 7 11 1 7 19 y(t; r) = − + e3t + ( + α)e−t , − + e3t + ( − α)e−t 3 12 4 3 12 4 10
International Journal of Mathematical Engineering and Science
ISSN : 2277-6982 http://www.ijmes.com/
Volume 1 Issue 5(May 2012) https://sites.google.com/site/ijmesjournal/
The solution of Runge-Kutta Nystrom method is as follows and Figure 1 and Table 1 show the obtain the results: Y
[m+1] 1
Y
[m+1] 1
Y
[m+1] 2
Y
[m+1] 2
Y
[m+1] 3
Y
[m+1] 3
=Y
[m] 1
=Y
[m] 1
=Y
[m] 2
=Y
[m] 2
=Y
[m] 3
=Y
[m] 3
+
! h3 [m] h+ Y2 + 2
! h2 h3 [m] Y3 , + 2 3
+
! h3 [m] h+ Y2 + 2
! h2 h3 [m] Y3 , + 2 3
+
! 3h2 [m] 3 +h Y2 + 2
! 7h3 [m] h+h + Y3 , 6
+
! 3h2 [m] + h3 Y 2 + 2
! 3 7h [m] h + h2 + Y3 , 6
+
! 7h3 [m] 3h + 3h + Y2 + 2
! 7h2 10h3 [m] 2h + Y3 , + 2 3
+
! 3 7h [m] 3h + 3h2 + Y2 + 2
! 7h2 10h3 [m] 2h + Y3 . + 2 3
2
2
Table 1. r 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Approximation y1 (ti ; r) y2 (ti ; r) 12.3925 13.0544 12.4293 13.0175 12.4659 12.9807 12.5025 12.9439 12.5394 12.9073 12.5763 12.8692 12.6130 12.8338 12.6499 12.7971 12.6865 12.7600 12.7233 12.7233
Exact y1 (ti ; r) y2 (ti ; r) 12.4317 13.0939 12.4685 13.0571 12.5053 13.0203 12.5421 12.9835 12.5788 12.9467 12.6156 12.9099 12.6524 12.8731 12.6892 12.8363 12.7260 12.7996 12.7628 12.7628
7 Conclusion In this paper a numerical method for solving nth -order fuzzy differential equation with fuzzy initial condition is presented. In this method nth -order fuzzy differential equation is converted to a fuzzy system which will be solved with the Runge-Kutta Nystrom method. 11
International Journal of Mathematical Engineering and Science
ISSN : 2277-6982 http://www.ijmes.com/
Volume 1 Issue 5(May 2012) https://sites.google.com/site/ijmesjournal/
RK Nystrom Method 1 0.9 * RK Nystrom − Exact
0.8 0.7
r
0.6 0.5 0.4 0.3 0.2 0.1 12.3
12.4
12.5
12.6
12.7
12.8
12.9
13
13.1
13.2
y
Figure 1: h=0.1
References [1] S. Abbasbandy, T. Allah Viranloo, Numerical solution of Fuzzy Differential Equations by Taylor method, J. Comput. Methods Appl. Math. 2(2002)113-124. [2] S. Abbasbandy, T. Allah Viranloo, Oscar Lopez-Pouso, J.J. Neito, Numerical Methods for Fuzzy Differential Inclusions, Comput. Math. Appl. 48(2004) 1633-1641. [3] T. Allah Viranloo, N. Ahmady, E. Ahmady, Numerical solution of Fuzzy Differential Equations by predictor-corrector method, Infor. Sci. 177(2007) 1633-1647. [4] T. Allah Viranloo, N. Ahmady, E. Ahmady, A method for solving nth -order Fuzzy Linear Differential Equations, Comput. Math. Appl 86(2009) 730-742. [5] S. Abbasbandy, T. Allah Viranloo, Numerical solution of N th -Order Fuzzy Differential Equations by Runge-Kutta Method Comput. Methods Appl. Math. 4 (2011) 935-946. [6] B. Bede, I.J. Rudas, A. L. Bencsik, First order Linear Fuzzy Differential Equations generalized differentiability, Inform. Sci. 177 (2007) 1648-1662. [7] J. J. Buckley, Thomas Feuring, Fuzzy Differential Equations, Fuzzy Sets and Systems 105 (1999), 241-248. [8] J. J. Buckley, Thomas Feuring, Introduction to Fuzzy Partial Differential Equations, Fuzzy Sets and Systems 110 (2000), 43-54. [9] J. J. Buckley, T. Feuring, Fuzzy initial value problem for N th -order fuzzy linear differential equations, Fuzzy Sets and Systems 121 (2001) 247-255. [10] S. L. Chang and A. Zadeh,On Fuzzy Mapping and Control, IEEEE Trans. Systems Man Cybernet., 2 (1972) 30-34. 12
International Journal of Mathematical Engineering and Science
ISSN : 2277-6982 http://www.ijmes.com/
Volume 1 Issue 5(May 2012) https://sites.google.com/site/ijmesjournal/
[11] W. Congxin, S. Shiji, Existence theorem to the Cauchy problem of Fuzzy Differential Equations under compactness-type conditions, Infor. Sci. 108(1998) 123-134. [12] D. Dubois and H. Prade,Towards Fuzy Differential Calculus: Part 3, Differentiation, Fuzzy Sets and Systems, 8 (1982) 225-233. [13] M. Friedman, M. Ming, A. Kandel, Fuzzy linear systems, Fuzzy Set and Systems 96 (1998) 201-209. [14] C.W. Gear, Numerical Initial Value Problem In Ordinary Differential Equations, Prentice Hall(1971). [15] R. Goetschel and W. Voxman, Elementary Fuzzy Calculs, Fuzzy Sets and Systems, 18 (1986) 31-43. [16] O. He, W. Yi, On fuzzy differential equations, Fuzzy Sets Systems, 24 (1989) 321-325. [17] T. Jayakumar, K. Kanagarajan, Numerical Solution for Hybrid Fuzzy System by Improved Euler Method, International Journal of Applied Mathematical Science 38 (2012) 1847-1862. [18] T. Jayakumar, D. Maheshkumar, K. Kanagarajan, Numerical Solution of Fuzzy Differential Equations by Runge-Kutta Method of Order Five, International Journal of Applied Mathematical Science 6 (2012) 2989-3002. [19] T. Jayakumar, K. Kanagarajan, Numerical Solution for Hybrid Fuzzy System by Runge-Kutta Method of Order Five, International Journal of Applied Mathematical Science 6 (2012) 3591-3606. [20] L.J. Jowers, J.J. Buckley. K.D. Reilly, Simulating continuous fuzzy systems, Infor. Sci. 177 (2007) 436-448. [21] O. Kaleva, Fuzzy Differential Equations, Fuzzy Sets Systems, 24 (1987) 301-317. [22] O. Kaleva, The Cauchy Problem for Fuzzy Differential Equations, Fuzzy Sets and Systems, 35 (1990) 389-396. [23] P. Kloeden, remarks on Peano-like theorems for fuzzy differential Equations, Fuzzy Sets Systems, 44 (1991) 161-164. [24] K. Kanagarajan, M.Sambath, Runge-Kutta Nystrom, method for solving fuzzy differential equations, Journal of computational methods in Applied mathemtics 10 (2010) 195-203. [25] A. Kandel, W.J. Byatt, Fuzzy differential equations, in Proceedings of the International Conference on Cybernetics and Society, Tokyo, (1978) 1213-1216. [26] M. L. Puri, D. A. Ralescu, Differentials of fuzzy functions, J. Math. Anal. Application 91 (1983),321-325. [27] D. Ralescu, A survey of the representation of fuzzy concepts and its applications, in: M.M Gupta, R.K. Ragade, R.R. Yager, Eds., Advances in fuzzy Set Theory and Applications, North-Holland, Amsterdam(1979) 77-91. [28] S. Seikkala, On the Fuzzy Initial Value Problem, Fuzzy Sets Systems, 24 (1987) 319-330.
13