Numerical Solution for Hybrid Fuzzy Systems by Milne s Fourth Order Predictor-Corrector Method
|
|
- Bartholomew York
- 5 years ago
- Views:
Transcription
1 International Mathematical Forum, Vol. 9, 2014, no. 6, HIKARI Ltd, Numerical Solution for Hybrid Fuzzy Systems by Milne s Fourth Order Predictor-Corrector Method T. Jayakumar and K. Kanagarajan Department of Mathematics Sri Ramakrishna Mission Vidyalaya College of Arts and Science Coimbatore , Tamilnadu, India Copyright c 2014 T. Jayakumar and K. Kanagarajan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract In this paper we study numerical methods for hybrid fuzzy differential equations by an appllication of the Milne s fourth order predictor- corrector method for fuzzy differential equations. We prove a convergence result and give numerical examples to illustrate the theory. Keywords: Hybrid Systems; Fuzzy Differential Equations; Milne s fourth order predictor-corrector method 1. Introduction Fuzzy set theory is a powerful tool for modeling uncertainty and processing vague or subjective information in mathematical models. Its main directions of development and its applications to the very varied real-world problems have been diverse. Fuzzy differential equations are recently gaining more and more attention in the literature. The first and the most popular approach in dealing with FDEs is using the Hukuhara differentiability, or the Seikkala derivative for fuzzy-number-valued functions. Hybrid systems are devoted to modelling, design, and validation of interactive systems of computer programs and continuous systems. That is, control systems that are capable of controlling complex systems which have discrete event dynamics as well as continuous time dynamics can be modelled by hybrid systems. The differential systems containing fuzzzy valued functions and interaction with a discrete time controller are named as hybrid fuzzy differential systems.
2 274 T. Jayakumar and K. Kanagarajan In the last few years, many works have been performed by several authors in numerical solutions of fuzzy differential equations [1, 2, 3, 9, 10, 11, 12, 15, 18]. Pederson ans Sambandam [20, 21]. have investigated the numerical solution of hybrid fuzzy differential equations by using Runge-Kutta and Euler methods. Recently, the numerical solutions of fuzzy differential equations by predictorcorrector method has been studied in [4]. Allah viranloo [5] discussed numerical soultion of fuzzy differential equation by Adams-Bashforth Two-step method. In this paper we develop numerical method for hybrid fuzzy differential equation by an application of the Milne s fourth order predictor-corrector method. In Section 2 we list some basic definitions to fuzzy valued functions. Section 3 reviews hybrid fuzzy differential systems. Section 4 contains the Milne s fourth order predictor-corrector method for hybrid fuzzy differential equations and convergence theorem. Section 5, the MPC four-step algorithm is discussed. Section 6, contains numerical examples to illustrate this method. Finally, conclution is present in Section Preliminary Notes Let P k (R n ) denote the family of all nonempty compact, convex subsets of R n. If α, β R and A, B P k (R n ), then α(a + B) =αa + αb, α(βa) = (αβ)a, 1A = A and if α, β 0, then (α + β)a = αa + βa. Denote by E n the set of u : R n [0, 1] such that u satisesfies (i) (iv) mentionned below: (i) u is normal, that is, there exists an x 0 R n such that u(x 0 )=1; (ii) u is fuzzy convex, that is, for x, y R and 0 λ 1, u(λx +(1 λ)y) min[u(x),u(y)]; (iii) u is upper semicontinuous; (iv) [u] 0 = the closure of [x R : u(x) > 0] is compact. For 0 <r 1, we define [u] r =[x R n : u(x) r]. Then from (i) to (iv), it follows that r - level sets [u] r =[x p k (R n ) for 0 r 1. An example of a u E 1 is given by 4x 3, if x (0.75, 1], u(x) = 2x +3, if x (1, 1.5), (1) 0, if x / (0.75, 1.5). The r-level sets of u in (1) are given by [u] r =[ r, r]. (2) Let I be real interval. A mapping g : I E is called a fuzzy process and its r-level set is denoted by [y(t)] r =[y r (t), y r ], t I, r (0, 1].
3 Numerical solution for hybrid fuzzy systems 275 Definition 2.1. Consider the initial value problem y (t) =f(t, y(t)), y(0) = y 0, (3) where f :[a, b] R n R n. An m-step method for solving the initial-value problem is one whose difference equation for finding the approximation y(t i+1 ) at the mesh point t i+1 can be represented by the following equation : y i+1 = m 1 j=0 a m j 1yi j + h m 1 j=0 b mf(t i j+1,y i j+1 ) (4) for i = m 1,m,,N 1, such that a = t 0 t 1 t N = b, h = (b a) N = t i+1 t i and a 0,,a m 1,b 0,,b m are constants with the starting values y 0 = α 0, y 1 = α 1, y m 1 = α m 1. When b m = 0, the method is known as explicit, since Equation (4) gives y i+1 explicit in terms of previously determined values. When b m 0, the method is known as implicit, since y i+1 occurs on both sides of Equation (4) and is specfied only implicitly. A special case of multistep method is. Here, we set q y i+1 = y i 1 + h k m m f(t i,y i ), q =0, 1, 2,..., (5) m=0 where the constants ( ) n s k m =( 1) m ds 0 k are independent of f,t = t 0 + sh, f(t i,y i ) is the first backward difference of the f(t, y(t)) at the point of t = t i and higher backward difference are defined by k f(t, y i )= ( k 1 f(t i,y i )). The special case q = 3 of Milnes method are as follows: Milne explicit method: y 0 = α 0, y 1 = α 1, y 2 = α 2, y 3 = α 3, y i+1 = y i 3 + 4h 3 [2f(t i,y i ) f(t i 1,y i 1 )+2f(t i 2,y i 2 )], where i =3, 4, 5,...,N 1. Milne implicit method : y 1 = α 1, y 2 = α 2, y 3 = α 3, y i+1 = y i 1 + h 3 [f(t i 1,y i 1 )+4f(t i,y i )+f(t i+1,y i+1 )]. where i =3, 4, 5,,N 1.
4 276 T. Jayakumar and K. Kanagarajan Definition 2.2. Assosisted with the difference equation y i+1 = a m 1 y i + a m 2 y i a 0 y i+1 m + hf (t i,h,y i+1,y i,,y i+1 m ), y 0 = α 0, y 1 = α 1, y m 1 = α m 1, (6) the characteristic polynomial of the method is defined by p(λ) =λ m a m 1 λ m 1 a m 2 λ m 2 a 1 λ a 0. If λ i 1 for each i =1, 2,,m, and all roots with absolute value 1 are simple roots, then the difference method is said to satisfy the root condition. Theorem 2.1. A multistep method of the form (4) is stable if and only if it satisfies the roots condition. Notations used in this paper are as follows: A tilde is placed over a symbol to denote a fuzzy set so α 1, f(t), An arbitrary fuzzy number with an ordered pair of functions (u(r), u(r)). 0 r 1, which satisfy the following requirements is represented. 1. u(r) is a bounded left continuous non-decreasing function over [0, 1], 2. u(r) is a bounded right continuous non-increasing function over [0, 1], 3. u(r) u (r), 0 r 1. Definition 2.3. The supremum metric d on E is defined by d (U, V )=sup {d H ([U] r, [V ] r ):r I}, and (E,d ) is a complete metric space. Definition 2.4. A mapping F : T E is a Hukuhara differentiable at t 0 T R if for some h 0 > 0 the Hukuhara differences F (t 0 +Δt) h F (t 0 ), F(t 0 ) h F (t 0 Δt), exist in E for all 0 < Δt <h 0 and if there exist an F (t 0 ) E such that ( ) F lim d (t0 +Δt) h F (t 0 ) Δt 0 + =0 Δt, F (t 0 ) and ( ) F lim d (t0 ) h F (t 0 Δt) Δt 0 + =0 Δt, F (t 0 ) the fuzzy set F (t 0 ) is called the Hukuhara derivative of F at t 0. Recall that U h V = W E are defined on level sets, where [U] r h[v ] r = [W ] r for all r I. By consideration of definition of the metric d, all the level set mappings [F (.)] r are Hukuhara differentiable at t 0 with Hukuhara
5 Numerical solution for hybrid fuzzy systems 277 derivatives [F (t 0 )] r for each r I when F : T E is Hukuhara differentiable at t 0 with Hukuhara derivative F (t 0 ). Definition 2.5. The fuzzy integral b a y(t)dt, 0 a b 1, is defined by [ b r [ b ] b y(t)dt] = y r (t)dt, y r (t)dt, a a a provided the Lebesgue integrals on the right exist. Definition 2.6. If F : T E is Hukuhara differentialble and its Hukuhara derivative F is integrable over [0,1], then F (t) =F (t 0 )+ t for all values of t 0,t where 0 t 0 t 1. t 0 F (s)ds, Definition 2.7. A mapping y : I E is called a fuzzy process. We denote [y(t)] r =[y r (t), y r (t)], t I, r (0, 1]. The Seikkala derivative y (t) of a fuzzy process y is defined by [y(t)] r = [y r (t), y r (t)], t I, r (0, 1] provided the equation defines a fuzzy number y (t) E. Definition 2.8. If y : I E is Seikkala differentiable and its Seikkala derivative y is integrable over [0, 1], then y(t) =y(t 0 )+ for all values of t 0,t where t 0,t I. t t 0 y (s)ds, 3. The Hybrid Fuzzy Differential System Consider the hybrid fuzzy differential system x (t) =f(t, x(t),λ k (x k )), t [t k,t k+1 ], x(t k )=x k, (7)
6 278 T. Jayakumar and K. Kanagarajan where denotes Seikkala differentiation, 0 t 0 <t 1 <... < t k <..., t k, f C[R + E 1 E 1,E 1 ],λ k C[E 1,E 1 ]. To be specific the system would look like x (t) = x 0 (t) =f(t, x 0 (t),λ 0 (x 0 )), x 0 (t 0 )=x 0, t 0 t t 1, x 1 (t) =f(t, x 1 (t),λ 1 (x 1 )), x 1 (t 1 )=x 1, t 1 t t 2,... x k (t) =f(t, x k (t),λ k (x k )), x k (t k )=x k, t k t t 2,... Assuming that the existence and uniquness of solution of (7) hold for each [t k,t k+1 ], by the solution of (7) we mean the following function: x 0 (t), t 0 t t 1, x 1 (t), t 1 t t 2,... x (t) = x k (t), t k t t k+1..., We note that the solution of (7) are piecewise differentiable in each interval for t [t k,t k+1 ] for a fixed x k E 1 and k =0, 1, 2,... Using a representation of fuzzy numbers studied by Goestschel and Woxman [8] and Wu and Ma [23], we may represent x E 1 by a pair of fuctions (x(r), x(r)), 0 r 1, such that (i) x(r) is bounded, left continuous, and nondecreasing, (ii) x(r) is bounded, left continuons, and nonincreasing, and (iii) x(r) x(r), 0 r 1. For example, u E 1 given in (1) is represented by (u(r), u(r))=( r, r), 0 r 1, which is similar to [u] α given by (2). Therefore we may replace (7) by an equivalent system x (t) =f(t, x, λ k (x k )) F k (t, x, x), x(t k )=x k, (8) x (t) =f(t, x, λ k (x k )) G k (t, x, x), x(t k )=x k, which possesses a unique solution (x, x) which is a fuzzy function. That is for each t, the pair [x(t; r), x(t; r)] is a fuzzy number, where x(t; r), x(t; r) are respectively the solutions of the parametric form given by x (t; r) =F k [t, x(t; r), x(t; r)], x(t k ; r) =x k (r), x (t; r) =G k [t, x(t; r), x(t; r)], x(t k ; r) =x k (r), (9)
7 Numerical solution for hybrid fuzzy systems Milne s fourth order predictor-corrector method In this section for a hybrid fuzzy differential equation (7). We develop the Milne s fourth order predictor-corrector method in [3]. When f and λ k in (7) can be obtained via, the Zadeh extension principal from f [R + R R, R] and λ k [R, R]. We assume that the existance and uniqueness of solution (7) hold for each [t k,t k+1 ]. For a fixed r, to integrate the system in (7) in [t 0,t 1 ], [t 1,t 2 ],..., [t k,t k+1 ],...,we replace each interval by a set of N k + 1 discreate equally spaced grid points (including the end points) at which the exact solution (x(t; r), x(t; r)) is approximated by some (y k (t; r), y k (t; r)). For the chosen grid points on [t k,t k+1 ] at t k,n = t k + nh k,h k = t k+1 t k N k, 0 n N k. Let (Y k (t; r), Y k (t; r)) (y(t, r), y(t, r)), (Y k (t; r), Y k (t; r))and (y k (t; r), y k (t; r)) may we denoted respctively by (Y k (t; r), Y k (t; r))and (y k (t; r), y k (t; r)). The Milne s fourth order predictor-corrector method is approximation of y k (t; r), and y k (t; r) which can be written as y 0 (r) =[α 0 (r), α 0 (r)], y 2 (r) =[α 2 (r), α 2 (r)], y 1 (r) =[α 1 (r), α 1 (r)] y 3 (r) =[α 3 (r), α 3 (r)] Y p n+1 = Y n 3 (r)+ 4h 3 [2f(t n,y n (r),λ k (u k )) f(t n 1,y n 1 (r),λ k (u k )) +2f(t n 2,y n 2 (r),λ k (u k ))], Y p n+1 = Y n 3(r)+ 4h 3 [2f(t n,y n (r),λ k (u k )) f(t n 1,y n 1 (r),λ k (u k )) +2f(t n 2,y n 2 (r),λ k (u k ))], (10) Y c n+1 = Y n 1(r)+ h[f(t 3 n 1,y n 1 (r),λ k (u k ))+4f(t n,y n (r),λ k (u k )) +f(t n+1,y n+1 (r),λ k (u k ))], Y c n+1 = Y n 1 (r)+ h [f(t 3 n 1,y n 1 (r),λ k (u k ))+4f(t n,y n (r),λ k (u k )) +f(t n+1,y n+1 (r),λ k (u k ))], where Y n (r) =[y n (r), y n (r)] f(t n,y n (r),λ k (u k )) = min{f(t, u, λ k (u k )) u [y(r), y n (r)]}, f(t n,y n (r),λ k (u k )) = max{f(t, u, λ k (u k )) u [y n (r), y n (r)]}.
8 280 T. Jayakumar and K. Kanagarajan Now we define F p [t, y(r)] = [2f(t n,y n (r),λ k (u k )) f(t n 1,y n 1 (r),λ k (u k )) +2f(t n 2,y n 2 (r),λ k (u k ))], G p [t, y(r)] = [2f(t n,y n (r),λ k (u k )) f(t n 1,y n 1 (r),λ k (u k )) +2f(t n 2,y n 2 (r),λ k (u k ))], F c [t, y(r)] = [f(t n 1,y n 1 (r),λ k (u k ))+4f(t n,y n (r),λ k (u k )) (11) +f(t n+1,y n+1 (r),λ k (u k ))], G c [t, y(r)] = [f(t n 1,y n 1 (r),λ k (u k ))+4f(t n,y n (r),λ k (u k )) +f(t n+1,y n+1 (r),λ k (u k ))]. The exact and approximate solutions at t n,0 n N are denoted by [Y (t n )] r = [Y (t; r), Y (t; r)] and [y(t n )] r = [y(t; r), y(t; r)] respectively. By (10),(11) we have Y c (t n+1 ; r) Y c (t n ; r)+ h 3 F [t n; Y c (t n ; r), Y c (t n ; r)] Y c (t n+1 ; r) Y c (t n ; r)+ h 3 G[t n; Y c (t n ; r), Y c (t n ; r)] (12) we define y c (t n+1 ; r) =y c (t n ; r)+ h 3 F [t n; y c (t n ; r), y c (t n ; r)] y c (t n+1 ; r) =y c (t n ; r)+ h 3 G[t n; y c (t n ; r), y c (t n ; r)] (13) However, (12) will use y 0,0 = x 0 (r), y 0,0 = x 0 (r), and y k.0 (r) =y k 1 (r)n k 1, y k.0 (r) =y k 1 (r)n k 1. If k 1, then (12) reperesents an approximation of y k (t; r), y k (t; r) for each of intervals t 0 t t 1,t 1 t t 2,...,t k t t k+1,... Lemma 4.1. Suppose i Z +,ɛ i > 0,r [0, 1], and h i < 1 are fixed. Let {Z i,n (r)} N i n=0 be the Milne s fourth order predictor-corrector approximation with N = N i to the fuzzy IVP: { x (t) =f(t, x(t),λ i (x i )), t [t i,t i+1 ], (14) x(t i )=x i.
9 Numerical solution for hybrid fuzzy systems 281 If {y i,n (r)} N i n=0 denotes the result of (10) from some y i,0 (r), then there exist a δ i > 0 such that z i,0 (r) y i,0 (r) <δ i, z i,0 (r) y i,0 (r) <δ i implies z i,0 (r) y i,0 (r) <ɛ i, z i,0 (r) y i,0 (r) <ɛ i proof: Proof is similar to lemma 3.1 [19] Theorem 4.1. r [0, 1], Proof. See [19]. Consider the systems (9) and (11). For a fixed k Z + and lim (r) =x(t h 0,...,h k 0 k,n k+1 ; r), k (15) lim k,n h 0,...,h k 0 k (r) =x(t k+1 ; r). (16) 5. Algorithm The following algorithm is based on Milne s fourth order predictor-corrector method. To approximate the solution of following fuzzy initial value problem y (t) =f(t, y(t),λ k (u k )), t 0 t T, y α = α 0, y α = α 1 y α = α 2, y α = α 3, y α = α 4, y α = α 5 y α = α 6, y α = α 7,, positive integer N is chosen Step 1. Let h = (T t 0), N y (t) =f(t, y(t),λ k (u k )), t 0 t T, w α = α 0, w α = α 1 w α = α 2, w α = α 3, w α = α 4, w α = α 5 w α = α 6, w α = α 7, Step 2. For n =0, 1, 2, do Steps 3-4. Compute starting values using fourth order Runge - Kutta method Step 3. Let k 1 (t k,n ; y k,n (r)) = min{h k.f(t k,n,u,λ k (u k )) u [y k,n (r), y k,n (r)]}, k 1 (t k,n,y k,n (r)) = max{h.f(t k,n,u,λ k (u k )) u, [y k,n (r), y k,n (r)]}, k 2 (t k,n,y k,n (r)) = min{h.f(t k,n + h/2,u,λ k (u k ) u [z 1 (t k,n,y k,n (r)), z 1 (t k,n,y k,n (r))]}, k 2 (t k,n,y k,n (r)) = max{h.f(t k,n + h/2,u,λ k (u k )) u [z 1 (t k,n,y k,n (r)), z 1 (t k,n,y k,n (r))]},
10 282 T. Jayakumar and K. Kanagarajan k 3 (t k,n,y k,n (r)) = min{h.f(t k,n + h/2,u,λ k (u k )) u [z 2 (t k,n,y k,n (r)), z 2 (t k,n,y k,n (r))]}, k 3 ( k,n,y k,n (r)) = max{h.f(t k,n + h/2,u,λ k (u k )) u [z 2 (t k,n,y k,n (r)), z 2 (t k,n,y k,n (r))]} k 4 (t k,n,y k,n (r)) = min{h.f(t k,n + h/2,u,λ k (u k )) u [z 3 (t k,n,y(r)), z 3 (t k,n,y k,n (r))]}, k 4 (t k,n,y k,n (r)) = max{h.f(t k,n + h, u, λ k (u k )) u [z 3 (t k,n,y k,n (r)), z 3 (t k,n,y k,n (r))]}, where z 1 (t k,n,y k,n (r)) = y k,n (r)+ 1 2 k 1(t k,n,y k,n (r)), z 1 (t k,n,y k,n (r)) = y k,n (r)+ 1 2 k 1(t, y k,n (r)), z 2 (t k,n,y k,n (r)) = y k,n (r)+ 1 2 k 2(t k,n,y k,n (r)), z 2 (t k,n,y k,n (r)) = y k,n (r)+ 1 2 k 2(t k,n,y k,n (r)), z 3 (t k,n,y k,n (r)) = y k,n (r)+k 3 (t k,n,y k,n (r)), z 3 (t k,n,y k,n (r)) = y k,n (t)+k 3 (t k,n,y k,n (r)). step 4. Let y(t k,n+1,r)=y(t k,n,y k,n (r)) + 1 (k 6 1(t k,n,y k,n (r))+2k 2 (t k,n,y k,n (r)) +2k 3 (t k,n,y k,n (r)) + k 4 (t k,n,y k,n (r))), y(t k,n+1,r)=y(t n,y k,n (r)) + 1 (k 6 1(t k,n,y k,n (r))+2k 2 (t k,n,y k,n (r)) +2k 3 (t k,n,y k,n (r)) + k 4 (t k,n,y k,n (r))). Step 5. For i=3,4,...,n-1. Let y p = y + 4h[2f(t i+1 i 3 3 n,y i,λ k (u k )) f(t i 1,y i 1,λ k (u k ))+2f(t i 2,y i 2,λ k (u k )) 9f(t i 3,y i 3 )λ k (u k )], y p i+1 = y i 3 + 4h 3 [2f(t i,y i,λ k (u k )) f(t i 1,y i 1,λ k (u k ))+2f(t i 2,y i 2,λ k (u k )) 9f(t i 3,y i 3 ),λ k (u k )].
11 Numerical solution for hybrid fuzzy systems 283 step 6. Let t i+2 = t 0 +(i +2)h Step 7. Let y c i+1 = y i 1 + h 3 [f(t i 1,y i 1,λ k (u k ))+4f(t i,y i,λ k (u k )) + f(t i+1,y i+1,λ k (u k )], y c i+1 = y i 1 + h 3 [f(t i 1,y i 1,λ k (u k ))+4f(t i,y i,λ k (u k )) + f(t i+1,y i+1 ),λ k (u k )] Step 8. For i = i +1 Step 9. If i N 2 go to step 3. Step 10. Algorithm will be completed and (w α (T ), w α (T )) approximates real value of (Y α (T ), Y α (T )). 6. Numerical Examples Before illustrating the numerical solution of hybrid fuzzy IVP, first we recall the Example 3 of [18] studying the fuzzy IVP; x (t) =x(t), x(0; r) =[ r, r], 0 r 1. (17) By Milne s predictor-corrector method with N=10. Since the exact solution of (17) is x(t; r) = [( r)e t, ( r)e t ], 0 r 1, we see that x(1; r) = [( r)e, ( r)e], 0 r 1. Example 6.1. Next consider the folllowing hybrid fuzzy IVP, { x (t) =x(t)+m(t)λ k x(t k ), t [t k,t k+1 ], t k = k, k =0, 1, 2, 3,..., x(0; r) = [( r), ( r)], 0 r 1, (18) where m(t) = { 2(t(mod1)) if t(mod1) 0.5, 2(1 t(mod1) if t(mod1) > 0.5, λ k (μ) = { ˆ0, if k =0 μ, if k {1, 2,...}. The hybrid fuzzy IVP (18) is equivalent to the following systems of fuzzy IVPs: x 0(t) =x 0 (t), t [0, 1], x(0; r) = [( r), ( r)], 0 r 1, x i (t) =x i (t)+m(t)x i 1 (t),t [t i,t i+1 ],x i (t) =x i 1 (t i ),i=1, 2,...,
12 284 T. Jayakumar and K. Kanagarajan In (14) x(t)+m(t)λ k (x(t k ) is continous function of t, x and λ k (x(t k ). Therefore by Example 6.1 of Kaleva [14], for each k =0, 1, 2,..., the fuzzy IVP { x (t) =x(t)+m(t)λ k (x(t k )), t [t k,t k+1 ],t k = k, x(t k )=x tk, (19) has a unique solution on [t k,t k+1 ]. Therefore x(1; r) = [( r)e, ( r)e],x(1; r)(3 e 3). Then x(1.5; 1) is approximately and y 1 (1.5; r) is approximitly For t [1.5; 2] the exact solution satisfies Therefore x(t; r) =x(1; r)(2t 2+e t 1.5 (3 e 4). x(2.0; r) =x(1; r)(2 + 3e 4 e). Then x(2.0; r) is approximately and y(2.0; 1) is approximately The exact and approximate solution by using Milne s predictor-corrector method. With N = 10 the following results are obtained. The exact and approximate solution by fourth order Runge-Kutta method and Milne s Predictor Corrector method are compared at t = 2 see Table 1 and Figure 1, Error in fourth order Runge-Kutta method and Milne s predictor-corrector method see Table 2. Table 1. r RK-Fourth Order Predictor Corrector Exact y 1 (t i ; r) y 2 (t i ; r) y 1 (t i ; r) y 2 (t i ; r) Y 1 (t i ; r) Y 2 (t i ; r)
13 Numerical solution for hybrid fuzzy systems 285 Table 2. r Rk 4 th -Order Predictor Corrector y 1 (t i ; r) y 2 (t i ; r) y 1 (t i ; r) y 2 (t i ; r) Predictor corrector method Exact.RK order4 0 PC Method r y Figure 1: h=0.1 Example 6.2. Next consider the following hybrid fuzzy IVP, x 0(t) =x 0 (t)+m(t)λ k x(t k ), t [0, 1], x(0; r) = [( r), ( r)], 0 r 1, x i (t) =x i (t)+m(t)x i 1 (t),t [t i,t i+1 ],x i (t) =x i 1 (t i ),i=1, 2,..., where m(t) = sin(πt), k =0, 1, 2,..., λ k (μ) = { ˆ0, if k =0 μ, if k {1, 2,...}.
14 286 T. Jayakumar and K. Kanagarajan Then x(t)+m(t)λ k (x(t k )) is continous function of t, x, and λ k (x(t k )). Therefore by Example of 6.1 of Kaleva [14], for each k=0,1,2,..., the fuzzy IVP { x (t) =x(t)+m(t)λ k (x(t k )), t [t k,t k+1 ],t k = k, x(t k )=x tk, (20) has a unique solution on [t k,t k+1 ]. To numerically solve the hybrid fuzzy IVP (18) we will apply the Milne s predictor-corrector method for hybrid fuzzy differential equation with N = 10 to obtain y 1,2 (r) approximating x(2.0; r). The exact solution of (20) satisfies x(t; r) =x(1; r) πcos(πt)+sin(πt) π 2 ( + et e x(1; r) 1+ π ), π 2 +1 x(t; r) =x(1; r) πcos(πt)+sin(πt) ( + et π 2 e x(1; r) 1+ π ). π 2 +1 Therefore x(1; r) = [( r)e, ( r)e], ( ( π x(2; r) = π e 1+ π )) x(1; r). π 2 +1 Then x(2; 1) is approximately where as y 1 (2.0; 1) is approximately see in Fig.2. Table 3. r RK-Fourth Order Predictor Corrector Exact y 1 (t i ; r) y 2 (t i ; r) y 1 (t i ; r) y 2 (t i ; r) y 1 (t i ; r) y 2 (t i ; r)
15 Numerical solution for hybrid fuzzy systems 287 Table 4. r Rk 4 th -Order Predictor Corrector y 1 (t i ; r) y 2 (t i ; r) y 1 (t i ; r) y 2 (t i ; r) Predictor corrector method Exact.RK order4 0 PC Method r y Figure 2: h= Conclution In this paper, we have applied itrative solution of Milne s predictor-corrector fourth order method for finding the numerical solution of hybrid fuzzy differential equations. Comparision of solution of Example (6.1) and (6.2) shows that our proposed method gives better solution then fourth order Runge-Kutta method.
16 288 T. Jayakumar and K. Kanagarajan References [1] S.Abbasbandy, T.Allah Viranloo, Numerical solution of fuzzy differential equation by Talor method, Journal of Computational Methods in Applied mathematics, 2 (2002) [2] S.Abbasbandy, T.Allah Viranloo, Numerical solution of fuzzy differential equation, Mathematical and Computational Applications, 7 (2002) [3] S.Abbasbandy, T.Allah Viranloo, Numerical solution of fuzzy differential equation by Runge-Kutta method, Nonlinear Studies, 11 (2004) [4] T.Allah Viranloo, N.Ahmady, E.Ahmady, Numerical solution of fuzzy differential equation by predictor-corrector method, Information Sciences 177 (2007) [5] T.Allah Viranloo, Numerical solution of fuzzy differential equation by Adams-Bashforth Two-Step Method, Journal of Applied Mathematics Islamic Azad University Lahijan (2004), [6] B.Bede, Note on Numerical solutions of fuzzy differential equations by predictor corrector method, Information science 178 (2008) [7] J.J.Bukley, T.Feuring, Fuzzy Differential Equations, Fuzzy Sets and Systems, 110 (2000) [8] R.Goetschel, Woxman, Elementary fuzzy calculus, Fuzzy Sets and Systems, 18 (1986) [9] T. Jayakumar, K. Kanagarajan, Numerical Solution for Hybrid Fuzzy System by Improved Euler Method, International Journal of Applied Mathematical Science, 38 (2012) [10] T. Jayakumar, K. Kanagarajan, S. Indrakumar, Numerical Solution of N th -Order Fuzzy Differential Equation by Runge-Kutta Nystrom Method, International Journal of Mathematical Engineering and Science, 1(5) [11] 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) [12] 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)
17 Numerical solution for hybrid fuzzy systems 289 [13] A.Khastan, K.Ivaz, Numerical Solution for Fuzzy Differential Equations by Nystrom method, Chaos, Solitons and Fractals 41 (2009) [14] O.Kaleva, Fuzzy Differential Equations, Fuzzy Sets and Systems, 24 (1987) [15] K. Kanagarajan, M.Sambath Numerical Solution for Fuzzy Differential Equations by Third Order Runge-Kutta Method, International Journal of Applied Mathematics and Computation, 2(4) (2010) [16] V.Lakshmikantham, X.Z.Liu, Impulsive Hybrid systems and stability theory,international Journalof Nonlinear Differential Equations, 5 (1999) [17] V.Lakshmikantham, R.N.Mohapatra, Theory of Fuzzy Differential Equations and Inclusions, Taylor and Francis, United Kingdom, [18] M.Ma, M.Friedman, A.Kandel, Numerical solutions of fuzzy differential equations, Fuzzy Sets and Systems, 105 (1999) [19] S.Pederson, M.Sambandham, Numerical solution to Hybrid fuzzy systems, Matheatical and Computer Modelling, 45(2007) [20] S.Pederson, M.Sambandham, The Runge-Kutta method for Hybrid Fuzzy Differential Equations, Nonlinear Analysis Hybrid Systems, 2 (2008) [21] M.Sambandham, Perturbed Lyapunov-like functions and Hybrid Fuzzy Differential Equations, International Journal of Hybrid Systems, 2 (2002) [22] S.Seikkala, On the fuzzy initial value problem, Fuzzy Sets and Systems, 24 (1987) [23] C.X.Wu, M.Ma, Embedding problem of fuzzy number space Part I, Fuzzy Sets and Systems, 44 (1991) Received: December 21, 2013
Numerical Solution of Hybrid Fuzzy Dierential Equation (IVP) by Improved Predictor-Corrector Method
Available online at http://ijim.srbiau.ac.ir Int. J. Industrial Mathematics Vol. 1, No. 2 (2009)147-161 Numerical Solution of Hybrid Fuzzy Dierential Equation (IVP) by Improved Predictor-Corrector Method
More informationTwo Step Method for Fuzzy Differential Equations
International Mathematical Forum, 1, 2006, no. 17, 823-832 Two Step Method for Fuzzy Differential Equations T. Allahviranloo 1, N. Ahmady, E. Ahmady Department of Mathematics Science and Research Branch
More informationNumerical Solution of Fuzzy Differential Equations
Applied Mathematical Sciences, Vol. 1, 2007, no. 45, 2231-2246 Numerical Solution of Fuzzy Differential Equations Javad Shokri Department of Mathematics Urmia University P.O. Box 165, Urmia, Iran j.shokri@mail.urmia.ac.ir
More informationNumerical Solution of Fuzzy Differential Equations of 2nd-Order by Runge-Kutta Method
Journal of Mathematical Extension Vol. 7, No. 3, (2013), 47-62 Numerical Solution of Fuzzy Differential Equations of 2nd-Order by Runge-Kutta Method N. Parandin Islamic Azad University, Kermanshah Branch
More informationM.S.N. Murty* and G. Suresh Kumar**
JOURNAL OF THE CHUNGCHEONG MATHEMATICAL SOCIETY Volume 19, No.1, March 6 THREE POINT BOUNDARY VALUE PROBLEMS FOR THIRD ORDER FUZZY DIFFERENTIAL EQUATIONS M.S.N. Murty* and G. Suresh Kumar** Abstract. In
More informationNUMERICAL SOLUTIONS OF FUZZY DIFFERENTIAL EQUATIONS BY TAYLOR METHOD
COMPUTATIONAL METHODS IN APPLIED MATHEMATICS, Vol.2(2002), No.2, pp.113 124 c Institute of Mathematics of the National Academy of Sciences of Belarus NUMERICAL SOLUTIONS OF FUZZY DIFFERENTIAL EQUATIONS
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationNUMERICAL SOLUTION OF A BOUNDARY VALUE PROBLEM FOR A SECOND ORDER FUZZY DIFFERENTIAL EQUATION*
TWMS J. Pure Appl. Math. V.4, N.2, 2013, pp.169-176 NUMERICAL SOLUTION OF A BOUNDARY VALUE PROBLEM FOR A SECOND ORDER FUZZY DIFFERENTIAL EQUATION* AFET GOLAYOĞLU FATULLAYEV1, EMINE CAN 2, CANAN KÖROĞLU3
More informationA Method for Solving Fuzzy Differential Equations Using Runge-Kutta Method with Harmonic Mean of Three Quantities
A Method for Solving Fuzzy Differential Equations Using Runge-Kutta Method with Harmonic Mean of Three Quantities D.Paul Dhayabaran 1 Associate Professor & Principal, PG and Research Department of Mathematics,
More informationNumerical Solving of a Boundary Value Problem for Fuzzy Differential Equations
Copyright 2012 Tech Science Press CMES, vol.86, no.1, pp.39-52, 2012 Numerical Solving of a Boundary Value Problem for Fuzzy Differential Equations Afet Golayoğlu Fatullayev 1 and Canan Köroğlu 2 Abstract:
More informationNumerical Method for Solving Second-Order. Fuzzy Boundary Value Problems. by Using the RPSM
International Mathematical Forum, Vol., 26, no. 4, 643-658 HIKARI Ltd, www.m-hikari.com http://d.doi.org/.2988/imf.26.6338 Numerical Method for Solving Second-Order Fuzzy Boundary Value Problems by Using
More informationResearch Article Adams Predictor-Corrector Systems for Solving Fuzzy Differential Equations
Mathematical Problems in Engineering Volume 2013, Article ID 312328, 12 pages http://dx.doi.org/10.1155/2013/312328 Research Article Adams Predictor-Corrector Systems for Solving Fuzzy Differential Equations
More informationSolution of the Fuzzy Boundary Value Differential Equations Under Generalized Differentiability By Shooting Method
Available online at www.ispacs.com/jfsva Volume 212, Year 212 Article ID jfsva-136, 19 pages doi:1.5899/212/jfsva-136 Research Article Solution of the Fuzzy Boundary Value Differential Equations Under
More informationMultistep Methods for IVPs. t 0 < t < T
Multistep Methods for IVPs We are still considering the IVP dy dt = f(t,y) t 0 < t < T y(t 0 ) = y 0 So far we have looked at Euler s method, which was a first order method and Runge Kutta (RK) methods
More informationA Novel Numerical Method for Fuzzy Boundary Value Problems
Journal of Physics: Conference Series PAPER OPEN ACCESS A Novel Numerical Method for Fuzzy Boundary Value Problems To cite this article: E Can et al 26 J. Phys.: Conf. Ser. 77 253 Related content - A novel
More informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationInitial-Value Problems for ODEs. Introduction to Linear Multistep Methods
Initial-Value Problems for ODEs Introduction to Linear Multistep Methods Numerical Analysis (9th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University
More informationMathematics for chemical engineers. Numerical solution of ordinary differential equations
Mathematics for chemical engineers Drahoslava Janovská Numerical solution of ordinary differential equations Initial value problem Winter Semester 2015-2016 Outline 1 Introduction 2 One step methods Euler
More informationNumerical Differential Equations: IVP
Chapter 11 Numerical Differential Equations: IVP **** 4/16/13 EC (Incomplete) 11.1 Initial Value Problem for Ordinary Differential Equations We consider the problem of numerically solving a differential
More informationInitial value problems for ordinary differential equations
Initial value problems for ordinary differential equations Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 IVP
More informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More informationRemark on the Sensitivity of Simulated Solutions of the Nonlinear Dynamical System to the Used Numerical Method
International Journal of Mathematical Analysis Vol. 9, 2015, no. 55, 2749-2754 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2015.59236 Remark on the Sensitivity of Simulated Solutions of
More informationOne-Step Hybrid Block Method with One Generalized Off-Step Points for Direct Solution of Second Order Ordinary Differential Equations
Applied Mathematical Sciences, Vol. 10, 2016, no. 29, 142-142 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.6127 One-Step Hybrid Block Method with One Generalized Off-Step Points for
More informationFUZZY SOLUTIONS FOR BOUNDARY VALUE PROBLEMS WITH INTEGRAL BOUNDARY CONDITIONS. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXXV 1(26 pp. 119 126 119 FUZZY SOLUTIONS FOR BOUNDARY VALUE PROBLEMS WITH INTEGRAL BOUNDARY CONDITIONS A. ARARA and M. BENCHOHRA Abstract. The Banach fixed point theorem
More informationA Study on Linear and Nonlinear Stiff Problems. Using Single-Term Haar Wavelet Series Technique
Int. Journal of Math. Analysis, Vol. 7, 3, no. 53, 65-636 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.988/ijma.3.3894 A Study on Linear and Nonlinear Stiff Problems Using Single-Term Haar Wavelet Series
More informationOn Nonlinear Methods for Stiff and Singular First Order Initial Value Problems
Nonlinear Analysis and Differential Equations, Vol. 6, 08, no., 5-64 HIKARI Ltd, www.m-hikari.com https://doi.org/0.988/nade.08.8 On Nonlinear Methods for Stiff and Singular First Order Initial Value Problems
More information8.1 Introduction. Consider the initial value problem (IVP):
8.1 Introduction Consider the initial value problem (IVP): y dy dt = f(t, y), y(t 0)=y 0, t 0 t T. Geometrically: solutions are a one parameter family of curves y = y(t) in(t, y)-plane. Assume solution
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationCHAPTER 5: Linear Multistep Methods
CHAPTER 5: Linear Multistep Methods Multistep: use information from many steps Higher order possible with fewer function evaluations than with RK. Convenient error estimates. Changing stepsize or order
More informationmultistep methods Last modified: November 28, 2017 Recall that we are interested in the numerical solution of the initial value problem (IVP):
MATH 351 Fall 217 multistep methods http://www.phys.uconn.edu/ rozman/courses/m351_17f/ Last modified: November 28, 217 Recall that we are interested in the numerical solution of the initial value problem
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationJim Lambers MAT 772 Fall Semester Lecture 21 Notes
Jim Lambers MAT 772 Fall Semester 21-11 Lecture 21 Notes These notes correspond to Sections 12.6, 12.7 and 12.8 in the text. Multistep Methods All of the numerical methods that we have developed for solving
More informationMath 128A Spring 2003 Week 12 Solutions
Math 128A Spring 2003 Week 12 Solutions Burden & Faires 5.9: 1b, 2b, 3, 5, 6, 7 Burden & Faires 5.10: 4, 5, 8 Burden & Faires 5.11: 1c, 2, 5, 6, 8 Burden & Faires 5.9. Higher-Order Equations and Systems
More informationResearch Article Diagonally Implicit Block Backward Differentiation Formulas for Solving Ordinary Differential Equations
International Mathematics and Mathematical Sciences Volume 212, Article ID 767328, 8 pages doi:1.1155/212/767328 Research Article Diagonally Implicit Block Backward Differentiation Formulas for Solving
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationOn interval predictor-corrector methods
DOI 10.1007/s11075-016-0220-x ORIGINAL PAPER On interval predictor-corrector methods Andrzej Marcinia 1,2 Malgorzata A. Janowsa 3 Tomasz Hoffmann 4 Received: 26 March 2016 / Accepted: 3 October 2016 The
More informationECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations
ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent
More informationTHIRD ORDER RUNGE-KUTTA METHOD FOR SOLVING DIFFERENTIAL EQUATION IN FUZZY ENVIRONMENT. S. Narayanamoorthy 1, T.L. Yookesh 2
International Journal of Pure Applied Mathematics Volume 101 No. 5 2015, 795-802 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu PAijpam.eu THIRD ORDER RUNGE-KUTTA
More informationSolving Ordinary Differential Equations
Solving Ordinary Differential Equations Sanzheng Qiao Department of Computing and Software McMaster University March, 2014 Outline 1 Initial Value Problem Euler s Method Runge-Kutta Methods Multistep Methods
More information1. Introduction: 1.1 Fuzzy differential equation
Numerical Solution of First Order Linear Differential Equations in Fuzzy Environment by Modified Runge-Kutta- Method and Runga- Kutta-Merson-Method under generalized H-differentiability and its Application
More informationZ. Omar. Department of Mathematics School of Quantitative Sciences College of Art and Sciences Univeristi Utara Malaysia, Malaysia. Ra ft.
International Journal of Mathematical Analysis Vol. 9, 015, no. 46, 57-7 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ijma.015.57181 Developing a Single Step Hybrid Block Method with Generalized
More informationNumerical solution of first order linear fuzzy differential equations using Leapfrog method
IOSR Journal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue 5 Ver. I (Sep-Oct. 4), PP 7- Numerical solution of first order linear fuzz differential equations using Leapfrog method
More informationLinear Multistep Methods I: Adams and BDF Methods
Linear Multistep Methods I: Adams and BDF Methods Varun Shankar January 1, 016 1 Introduction In our review of 5610 material, we have discussed polynomial interpolation and its application to generating
More informationRunge-Kutta Method for Solving Uncertain Differential Equations
Yang and Shen Journal of Uncertainty Analysis and Applications 215) 3:17 DOI 1.1186/s4467-15-38-4 RESEARCH Runge-Kutta Method for Solving Uncertain Differential Equations Xiangfeng Yang * and Yuanyuan
More informationNew Nonlinear Conditions for Approximate Sequences and New Best Proximity Point Theorems
Applied Mathematical Sciences, Vol., 207, no. 49, 2447-2457 HIKARI Ltd, www.m-hikari.com https://doi.org/0.2988/ams.207.7928 New Nonlinear Conditions for Approximate Sequences and New Best Proximity Point
More information2 Numerical Methods for Initial Value Problems
Numerical Analysis of Differential Equations 44 2 Numerical Methods for Initial Value Problems Contents 2.1 Some Simple Methods 2.2 One-Step Methods Definition and Properties 2.3 Runge-Kutta-Methods 2.4
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Linear Multistep methods Gustaf Söderlind and Carmen Arévalo Numerical Analysis, Lund University Textbooks: A First Course in the
More informationChapter 2. Fuzzy function space Introduction. A fuzzy function can be considered as the generalization of a classical function.
Chapter 2 Fuzzy function space 2.1. Introduction A fuzzy function can be considered as the generalization of a classical function. A classical function f is a mapping from the domain D to a space S, where
More informationResearch Article On λ-statistically Convergent Double Sequences of Fuzzy Numbers
Hindawi Publishing Corporation Journal of Inequalities and Applications Volume 2008, Article ID 47827, 6 pages doi:0.55/2008/47827 Research Article On λ-statistically Convergent Double Sequences of Fuzzy
More informationResearch Article Solution of Fuzzy Matrix Equation System
International Mathematics and Mathematical Sciences Volume 2012 Article ID 713617 8 pages doi:101155/2012/713617 Research Article Solution of Fuzzy Matrix Equation System Mahmood Otadi and Maryam Mosleh
More informationPartial Averaging of Fuzzy Differential Equations with Maxima
Advances in Dynamical Systems and Applications ISSN 973-5321, Volume 6, Number 2, pp. 199 27 211 http://campus.mst.edu/adsa Partial Averaging o Fuzzy Dierential Equations with Maxima Olga Kichmarenko and
More informationDisconvergent and Divergent Fuzzy Sequences
International Mathematical Forum, Vol. 9, 2014, no. 33, 1625-1630 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2014.49167 Disconvergent and Divergent Fuzzy Sequences M. Muthukumari Research
More information1 Error Analysis for Solving IVP
cs412: introduction to numerical analysis 12/9/10 Lecture 25: Numerical Solution of Differential Equations Error Analysis Instructor: Professor Amos Ron Scribes: Yunpeng Li, Mark Cowlishaw, Nathanael Fillmore
More informationSolving Two-Dimensional Fuzzy Partial Dierential Equation by the Alternating Direction Implicit Method
Available online at http://ijim.srbiau.ac.ir Int. J. Instrial Mathematics Vol. 1, No. 2 (2009) 105-120 Solving Two-Dimensional Fuzzy Partial Dierential Equation by the Alternating Direction Implicit Method
More informationSOLVING FUZZY DIFFERENTIAL EQUATIONS BY USING PICARD METHOD
Iranian Journal of Fuzzy Systems Vol. 13, No. 3, (2016) pp. 71-81 71 SOLVING FUZZY DIFFERENTIAL EQUATIONS BY USING PICARD METHOD S. S. BEHZADI AND T. ALLAHVIRANLOO Abstract. In this paper, The Picard method
More informationReview Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester
HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE
More informationMTH 452/552 Homework 3
MTH 452/552 Homework 3 Do either 1 or 2. 1. (40 points) [Use of ode113 and ode45] This problem can be solved by a modifying the m-files odesample.m and odesampletest.m available from the author s webpage.
More information9.6 Predictor-Corrector Methods
SEC. 9.6 PREDICTOR-CORRECTOR METHODS 505 Adams-Bashforth-Moulton Method 9.6 Predictor-Corrector Methods The methods of Euler, Heun, Taylor, and Runge-Kutta are called single-step methods because they use
More informationCOSC 3361 Numerical Analysis I Ordinary Differential Equations (II) - Multistep methods
COSC 336 Numerical Analysis I Ordinary Differential Equations (II) - Multistep methods Fall 2005 Repetition from the last lecture (I) Initial value problems: dy = f ( t, y) dt y ( a) = y 0 a t b Goal:
More informationThe Generalized Viscosity Implicit Rules of Asymptotically Nonexpansive Mappings in Hilbert Spaces
Applied Mathematical Sciences, Vol. 11, 2017, no. 12, 549-560 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.718 The Generalized Viscosity Implicit Rules of Asymptotically Nonexpansive
More informationRemarks on Fuzzy Differential Systems
International Journal of Difference Equations ISSN 097-6069, Volume 11, Number 1, pp. 19 6 2016) http://campus.mst.edu/ijde Remarks on Fuzzy Differential Systems El Hassan Eljaoui Said Melliani and Lalla
More informationKKM-Type Theorems for Best Proximal Points in Normed Linear Space
International Journal of Mathematical Analysis Vol. 12, 2018, no. 12, 603-609 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2018.81069 KKM-Type Theorems for Best Proximal Points in Normed
More informationProx-Diagonal Method: Caracterization of the Limit
International Journal of Mathematical Analysis Vol. 12, 2018, no. 9, 403-412 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2018.8639 Prox-Diagonal Method: Caracterization of the Limit M. Amin
More informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationSequences of (ψ, φ)-weakly Contractive Mappings and Stability of Fixed Points
Int. Journal of Math. Analysis, Vol. 7, 2013, no. 22, 1085-1096 HIKARI Ltd, www.m-hikari.com Sequences of (ψ, φ)-weakly Contractive Mappings and Stability of Fixed Points S. N. Mishra 1, Rajendra Pant
More informationCommon Fixed Point Theorems for Ćirić-Berinde Type Hybrid Contractions
International Journal of Mathematical Analysis Vol. 9, 2015, no. 31, 1545-1561 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2015.54125 Common Fixed Point Theorems for Ćirić-Berinde Type
More informationMorera s Theorem for Functions of a Hyperbolic Variable
Int. Journal of Math. Analysis, Vol. 7, 2013, no. 32, 1595-1600 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2013.212354 Morera s Theorem for Functions of a Hyperbolic Variable Kristin
More informationStrong Convergence of the Mann Iteration for Demicontractive Mappings
Applied Mathematical Sciences, Vol. 9, 015, no. 4, 061-068 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ams.015.5166 Strong Convergence of the Mann Iteration for Demicontractive Mappings Ştefan
More informationA method for solving first order fuzzy differential equation
Available online at ttp://ijim.srbiau.ac.ir/ Int. J. Industrial Matematics (ISSN 2008-5621) Vol. 5, No. 3, 2013 Article ID IJIM-00250, 7 pages Researc Article A metod for solving first order fuzzy differential
More informationDiophantine Equations. Elementary Methods
International Mathematical Forum, Vol. 12, 2017, no. 9, 429-438 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/imf.2017.7223 Diophantine Equations. Elementary Methods Rafael Jakimczuk División Matemática,
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationSymmetric Properties for the (h, q)-tangent Polynomials
Adv. Studies Theor. Phys., Vol. 8, 04, no. 6, 59-65 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/astp.04.43 Symmetric Properties for the h, q-tangent Polynomials C. S. Ryoo Department of Mathematics
More informationA Generalization of Generalized Triangular Fuzzy Sets
International Journal of Mathematical Analysis Vol, 207, no 9, 433-443 HIKARI Ltd, wwwm-hikaricom https://doiorg/02988/ijma2077350 A Generalization of Generalized Triangular Fuzzy Sets Chang Il Kim Department
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationAN OVERVIEW. Numerical Methods for ODE Initial Value Problems. 1. One-step methods (Taylor series, Runge-Kutta)
AN OVERVIEW Numerical Methods for ODE Initial Value Problems 1. One-step methods (Taylor series, Runge-Kutta) 2. Multistep methods (Predictor-Corrector, Adams methods) Both of these types of methods are
More informationNew Multi-Step Runge Kutta Method For Solving Fuzzy Differential Equations
Abstract New Multi-Step Runge Kutta Method For Solving Fuzzy Differential Equations Nirmala. V 1* Chenthur Pandian.S 2 1. Department of Mathematics, University College of Engineering Tindivanam (Anna University
More informationLecture 5: Single Step Methods
Lecture 5: Single Step Methods J.K. Ryan@tudelft.nl WI3097TU Delft Institute of Applied Mathematics Delft University of Technology 1 October 2012 () Single Step Methods 1 October 2012 1 / 44 Outline 1
More information5.6 Multistep Methods
5.6 Multistep Methods 1 Motivation: Consider IVP: yy = ff(tt, yy), aa tt bb, yy(aa) = αα. To compute solution at tt ii+1, approximate solutions at mesh points tt 0, tt 1, tt 2, tt ii are already obtained.
More informationFuzzy Sequences in Metric Spaces
Int. Journal of Math. Analysis, Vol. 8, 2014, no. 15, 699-706 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.4262 Fuzzy Sequences in Metric Spaces M. Muthukumari Research scholar, V.O.C.
More informationOrdinary Differential Equations
CHAPTER 8 Ordinary Differential Equations 8.1. Introduction My section 8.1 will cover the material in sections 8.1 and 8.2 in the book. Read the book sections on your own. I don t like the order of things
More informationRemark on a Couple Coincidence Point in Cone Normed Spaces
International Journal of Mathematical Analysis Vol. 8, 2014, no. 50, 2461-2468 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.49293 Remark on a Couple Coincidence Point in Cone Normed
More informationGeneralized Functions for the Fractional Calculus. and Dirichlet Averages
International Mathematical Forum, Vol. 8, 2013, no. 25, 1199-1204 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2013.3483 Generalized Functions for the Fractional Calculus and Dirichlet Averages
More informationA STUDY OF GENERALIZED ADAMS-MOULTON METHOD FOR THE SATELLITE ORBIT DETERMINATION PROBLEM
Korean J Math 2 (23), No 3, pp 27 283 http://dxdoiorg/568/kjm232327 A STUDY OF GENERALIZED ADAMS-MOULTON METHOD FOR THE SATELLITE ORBIT DETERMINATION PROBLEM Bum Il Hong and Nahmwoo Hahm Abstract In this
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationOrdinary differential equation II
Ordinary Differential Equations ISC-5315 1 Ordinary differential equation II 1 Some Basic Methods 1.1 Backward Euler method (implicit method) The algorithm looks like this: y n = y n 1 + hf n (1) In contrast
More informationLecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More informationNotes for Numerical Analysis Math 5466 by S. Adjerid Virginia Polytechnic Institute and State University (A Rough Draft) Contents Numerical Methods for ODEs 5. Introduction............................
More informationSolving Systems of Fuzzy Differential Equation
International Mathematical Forum, Vol. 6, 2011, no. 42, 2087-2100 Solving Systems of Fuzzy Differential Equation Amir Sadeghi 1, Ahmad Izani Md. Ismail and Ali F. Jameel School of Mathematical Sciences,
More informationThe collocation method for ODEs: an introduction
058065 - Collocation Methods for Volterra Integral Related Functional Differential The collocation method for ODEs: an introduction A collocation solution u h to a functional equation for example an ordinary
More information2.1 NUMERICAL SOLUTION OF SIMULTANEOUS FIRST ORDER ORDINARY DIFFERENTIAL EQUATIONS. differential equations with the initial values y(x 0. ; l.
Numerical Methods II UNIT.1 NUMERICAL SOLUTION OF SIMULTANEOUS FIRST ORDER ORDINARY DIFFERENTIAL EQUATIONS.1.1 Runge-Kutta Method of Fourth Order 1. Let = f x,y,z, = gx,y,z be the simultaneous first order
More informationInitial value problems for ordinary differential equations
AMSC/CMSC 660 Scientific Computing I Fall 2008 UNIT 5: Numerical Solution of Ordinary Differential Equations Part 1 Dianne P. O Leary c 2008 The Plan Initial value problems (ivps) for ordinary differential
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
More informationApproximating solutions of nonlinear second order ordinary differential equations via Dhage iteration principle
Malaya J. Mat. 4(1)(216) 8-18 Approximating solutions of nonlinear second order ordinary differential equations via Dhage iteration principle B. C. Dhage a,, S. B. Dhage a and S. K. Ntouyas b,c, a Kasubai,
More informationThe Greatest Common Divisor of k Positive Integers
International Mathematical Forum, Vol. 3, 208, no. 5, 25-223 HIKARI Ltd, www.m-hiari.com https://doi.org/0.2988/imf.208.822 The Greatest Common Divisor of Positive Integers Rafael Jaimczu División Matemática,
More informationInternational Mathematical Forum, Vol. 9, 2014, no. 36, HIKARI Ltd,
International Mathematical Forum, Vol. 9, 2014, no. 36, 1751-1756 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2014.411187 Generalized Filters S. Palaniammal Department of Mathematics Thiruvalluvar
More informationOn Pointwise λ Statistical Convergence of Order α of Sequences of Fuzzy Mappings
Filomat 28:6 (204), 27 279 DOI 0.2298/FIL40627E Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat On Pointwise Statistical Convergence
More informationMath 128A Spring 2003 Week 11 Solutions Burden & Faires 5.6: 1b, 3b, 7, 9, 12 Burden & Faires 5.7: 1b, 3b, 5 Burden & Faires 5.
Math 128A Spring 2003 Week 11 Solutions Burden & Faires 5.6: 1b, 3b, 7, 9, 12 Burden & Faires 5.7: 1b, 3b, 5 Burden & Faires 5.8: 1b, 3b, 4 Burden & Faires 5.6. Multistep Methods 1. Use all the Adams-Bashforth
More informationSolving Ordinary Differential equations
Solving Ordinary Differential equations Taylor methods can be used to build explicit methods with higher order of convergence than Euler s method. The main difficult of these methods is the computation
More informationCS520: numerical ODEs (Ch.2)
.. CS520: numerical ODEs (Ch.2) Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca people.cs.ubc.ca/ ascher/520.html Uri Ascher (UBC) CPSC 520: ODEs (Ch. 2) Fall
More information