Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 803462, 10 pages

http://dx.doi.org/10.1155/2013/803462

## Numerical Solution of Second-Order Fuzzy Differential Equation Using Improved Runge-Kutta Nystrom Method

^{1}Mathematics Department, Science Faculty, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia^{2}Institute for Mathematical Research, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia^{3}Young Researchers and Elite Club, Mobarakeh Branch, Islamic Azad University, Mobarakeh, Iran

Received 5 February 2013; Accepted 25 March 2013

Academic Editor: Tofigh Allahviranloo

Copyright © 2013 Faranak Rabiei et al. 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

We develop the Fuzzy Improved Runge-Kutta Nystrom (FIRKN) method for solving second-order fuzzy differential equations (FDEs) based on the generalized concept of higher-order fuzzy differentiability. The scheme is two-step in nature and requires less number of stages which leads to less number of function evaluations in comparison with the existing Fuzzy Runge-Kutta Nystrom method. Therefore, the new method has a lower computational cost which effects the time consumption. We assume that the fuzzy function and its derivative are Hukuhara differentiable. FIRKN methods of orders three, four, and five are derived with two, three, and four stages, respectively. The numerical examples are given to illustrate the efficiency of the methods.

#### 1. Introduction

Fuzzy differential equations serve as mathematical models for many exciting *real-world* problems, not only in science and technology but also in such diverse fields as population models [1], civil engineering [2], and modelling hydraulic [3].

Initially, the derivative of fuzzy-valued functions was first introduced by Chang and Zadeh [4]. It was followed by Dubois and Prade [5], who used the extension principle in their approach. Other methods have been discussed by Puri and Ralescu [6]; they generalized and extended the concept of Hukuhara differentiability (-derivative) from set-valued mappings to the class of fuzzy mappings.

Subsequently, using -derivative, Kaleva [7, 8] and Seikkala [9] developed the theory of fuzzy differential equations.

In the last few years, many researchers have worked on theoretical and numerical solution of FDEs [10–23], specially some authors considered the second-order fuzzy differential equations [24–26].

To our best knowledge up to now, a few investigations have been devoted to the numerical solution of second order fuzzy differential equations. In this paper, a novel Runge-kutta Nystrom method is proposed for solving fuzzy differential equations of second order which is constructed based on the used of the previous values of acquired in the last stage. This Algorithm initially has been proposed by Rabiei et al. [27] for solving second-order ordinary differential equations which was the extension of the crisp concept of this method in solving first-order ODEs given in [28, 29]. The most important advantage of this method is that it has a lower computational cost in comparison with the previous findings for methods of the same order which improved the efficiency.

The aim of this paper is to introduce the fuzzy extension of Improved Runge-Kutta Nystrom method in solving second-order ODEs given in [27]. The second-order FODEs are assumed under Hukuhara differentiability. Actually, we suppose that fuzzy function and its derivative are -differentiable. We, therefore, motivated our interest in the examples under this assumption. The accuracy of the proposed algorithms is demonstrated by test problems. Additionally, many formulas corresponding to mentioned references previously are applied for solving other kinds of FODEs and are traceless in the literature for second-order ODEs in the fuzzy sense which is another motivation for developing Improved Runge-Kutta Nystrom method for solving second order FODEs.

The paper is organized as follows. In Section 2, we give some basic definitions and theorem on FDEs. In Section 3, Fuzzy Improved Runge-Kutta Nystrom method of orders 3, 4, and 5 are proposed. In Section 4, the numerical examples are provided to illustrate the validity and applicability of the new method. Finally, some conclusions are given.

#### 2. Preliminaries

We give some definitions and introduce the necessary notation which will be used throughout the paper; see [30, 31].

We consider , the set of all real numbers. A fuzzy number is mapping with the following properties:(a)is upper semicontinuous;(b) is fuzzy convex, that is, for all ;(c) is normal, that is, for which ;(d) is the support of the , and its closure is compact.

Let be the set of all fuzzy numbers on . The -level set of a fuzzy number , denoted by , is defined as It is clear that -level set of a fuzzy number is a closed and bounded interval , where denotes the left-hand endpoint of and denotes the right-hand endpoint of . Since each can be regarded as a fuzzy number is defined by

For and , the sum and the product are defined by

, for all , where means that usual addition of two intervals (subsets) of and means the usual product between a scalar and a subset of .

The Hausdorff distance fuzzy numbers are given by : It is easy to see that is a metric in and has the following properties:(i), for all ,(ii), for all ,(iii), for all ,(iv) is a complete metric space.

*Definition 1. *Let be a fuzzy valued function. If for arbitrary fixed and such that
is said to be continuous.

Initially the -derivative (Hukuhara differentiability) for fuzzy mappings was introduced by Puri and Ralescu [6] which is based on the -difference sets, as follows.

*Definition 2. *Let . If there exists such that , then is called the -difference of and , and it is denoted by .

In this paper, the sign “” stands for -difference, and also note that .

*Definition 3. *Let be a fuzzy function. We say is differentiable at , if there exists an element such that limits
exist and are equal to . Here the limits are taken in the metric space , since we have defined and .

*Definition 4 (see [32]). *Let and . We say that is strongly generalized differential at . If there exists an element , such that(i) for all sufficiently small, and the limits (in the metric )
(ii) for all sufficiently small, and the limits (in the metric )
(iii) for all sufficiently small, and the limits (in the metric )
(iv) for all sufficiently small, and the limits (in the metric )

*Remark 5. *In [32], the authors consider four cases for derivatives. Here we only consider the two first cases of Definition 4. In the other cases, the derivative is trivial because it is reduced to a crisp element. We say is (1)-differentiable on if is differentiable with the meaning (i) of Definition 4 and also (2)-differentiable that satisfies in the Definition 7 case (ii).

Theorem 6 (see [33]). *Let be a function and denote , for each . Then*(1)* if is (1)-differentiable, then and are differentiable functions and
*(2)* if is (2)-differentiable, then and are differentiable functions and
*

##### 2.1. Fuzzy Initial Value Problem

Consider the second-order fuzzy initial value problem: where is a fuzzy function with -level sets of initial values We write and where By using the extension principle, when is a fuzzy number we have the membership function It follows thatwhere

*Definition 7 (see [26]). *Let and . We say that is strongly generalized differentiable of the second order at , if there exists an element , such that

(i) for all sufficiently small, and the limits (in the metric )

or (ii) for all sufficiently small, and the limits (in the metric )

or (iii) for all sufficiently small, and the limits (in the metric )

or (iv) for all sufficiently small, and the limits (in the metric )

Proposition 8. *For a supposed fuzzy function , one has two possibilities, according to Definitions 4 and 7, to obtain the derivative of over : and . Then for each of these two derivatives, one has again two possibilities: and , respectively.*

*Remark 9. *In the rest of paper, we assume that and are (1)-differentiable.

#### 3. Fuzzy Improved Runge-Kutta Nystrom Method (FIRKN)

Improved Runge-Kutta (IRK) method for solving first-order ordinary differential equations was constructed by Rabiei and Ismail [28, 29]. Later, they developed the Improved Runge-Kutta Nystrom (IRKN) method for solving special second order ODEs. This scheme is two-step in nature and requires less number of stages which leads to less number of function evaluations compared with existing Runge-Kutta Nystrom methods. The general form of IRKN method with -stages is given by (see [27]) where and depends on both and while and depend on the values of and for . Here is the number of function evaluations performed at each step and increases with the order of local accuracy of the IRK method. In each step we only need to evaluate the values of , while are calculated from the previous step.

Based on IRKN method in (22), we proposed the fuzzy version of IRKN method which is denoted by FIRKN. Here, FIRKN methods of orders three, four, and five with two, three, and four stages, respectively, are derived. The coefficients of FIRKN methods are the same as the coefficients of the IRKN methods.

Let the exact solution be approximated by . We define Note that the values of and in each step are replaced by and from the previous step; therefore, there is no need to evaluate them again.

##### 3.1. Fuzzy Improved Runge-Kutta Nystrom Method of Order Three

Based on formulas (22), we define the FIRKN3 with two stages as follows:

whereHere, the coefficients of FIRKN3 are the same as the coefficients of IRKN3 which are given as follows:

##### 3.2. Fuzzy Improved Runge-Kutta Nystrom Method of Order Four

From formulas (22), consider the FIRKN4 with three stages as follows: where

The coefficients of FIRKN4 are the same as the coefficients of IRKN4 which are given as follows:

##### 3.3. Fuzzy Improved Runge-Kutta Nystrom Method of Order Five

From formulas (22), consider Fuzzy Improved Runge-Kutta Nystrom method of order five (FIRKN5) with four stages as follows: where

Again, the coefficients of FIRKN5 are the same as the coefficients of IRKN5 which are given as the following:

#### 4. Numerical Examples

In this section, we solved the fuzzy initial value problems to show the efficiency and accuracy of the proposed methods. Let the exact solution be . The absolute error formula, considered in Tables 1–4, is as follows:

*Problem 1 (see [13]). *Consider

The exact solution using (1)-differentiability is:

*Problem 2 (see [34]). *Consider

The exact solution under (1)-differentiability: where
To illustrate the efficiency of FIRKN methods, we compared the numerical results of the new methods with the existing method and numerical results are given in Tables 1–4 and Figures 1–4. In addition, the following abbreviations are used in Tables 1–4.(i)FIRKN3: third-order Fuzzy Improved Runge-Kutta Nystrom method with two stages derived in this paper. (ii)FIRKN4: fourth-order Fuzzy Improved Runge-Kutta Nystrom method with three stages derived in this paper. (iii)FIRKN5: fifth-order Fuzzy Improved Runge-Kutta Nystrom method with four stages derived in this paper. (iv)FRKN4: fourth-order Fuzzy Runge-Kutta Nystrom method with four stages based on the fourth-order Runge-Kutta Nystrom method with order ten dispersion and dissipation order five, derived by Houwen and Sommeijer (see [35]).(v): number of stages.

For Problem 1, the absolute error of FIRKN3, FIRKN4, and FIRKN5 compared with FRKN4 for and are given in Tables 1 and 2, respectively. FIRKN3 used only two stages, but the numerical results for both and are comparable with the results produced by FRKN4 with 4 stages. We can say that FIRKN3 is as accurate as FRKN4 but with less number of stages; hence, it is computationally more efficient.

FIRKN4 with three stages given accuracy up to which is two orders higher than FRKN4. Hence, FIRKN4 is more accurate with less number of functions of evaluations compared to FRKN4.

FIRKN5 with four stapes gives accuracy up to while FRKN4 with the same number of stages is accurate up to , that indicates the efficiency of FIRKN methods in solving Problem 1.

Figure 1, shows the approximate solutions using FIRKN methods tend the exact solutions. In Figure 2, the graph of fuzzy function is plotted and we can see that the approximated solutions of are between the fuzzy solutions and ; thus, the approximated solutions are valid for fuzzy level set.

From the numerical results of Problem 2 in Tables 3 and 4, it is observed that FIRKN3 which used less number of stages than FRKN4 gives a suitable accuracy. FIRKN4 and FIRKN5 give produced accuracy up to while the accuracy of FRKN4 is only up to . Also, Figure 3, shows that the approximated solutions by FIRKN methods are close to the exact solutions and the validity of the approximated solutions of for fuzzy level set in solving Problem 2 is given by Figure 4.

#### 5. Conclusion

In this paper we developed the Fuzzy Improved Runge-Kutta Nystrom methods for solving second-order fuzzy differential equations. The scheme is based on the Improved Runge-Kutta Nystrom method for solving second-order ordinary differential equations. The methods of orders three, four, and five with two, three and four stages, respectively, are extended.

FIRKN methods used less number of stages which leads to the less number of function evaluations have a lower computational cost. Therefore, we can conclude that the Fuzzy Improved Runge-Kutta Nystrom methods with high accuracy and less number of function evaluations compared with the existing fuzzy Runge-Kutta Nystrom methods are more efficient for solving second-order fuzzy differential equations.

The presented method in this research may be useful if the coefficients, initial values are fuzzy and it can be applied by choosing the different types of fuzzy derivatives except -differentiability.

#### Acknowledgments

The authors thank the reviewers and Leader Guest Editor, Professor Allahviranloo, for their useful comments which helped improve the quality of the paper.

#### References

- M. Guo, X. Xue, and R. Li, “The oscillation of delay differential inclusions and fuzzy biodynamics models,”
*Mathematical and Computer Modelling*, vol. 37, no. 7-8, pp. 651–658, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Guo, X. Xue, and R. Li, “Impulsive functional differential inclusions and fuzzy population models,”
*Fuzzy Sets and Systems*, vol. 138, no. 3, pp. 601–615, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Bencsik, B. Bede, J. Tar, and J. Fodor, “Fuzzy differential equations in modeling hydraulic differential servo cylinders,” in
*Proceedings of the 3rd Romanian-Hungarian Joint Symposium on Applied Computational Intellidence (SACI '06)*, Timisoara, Romania, 2006. - S. S. L. Chang and L. A. Zadeh, “On fuzzy mapping and control,”
*IEEE Transactions on Systems, Man, and Cybernetics*, vol. 2, pp. 30–34, 1972. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Dubois and H. Prade, “Towards fuzzy differential calculus. III. Differentiation,”
*Fuzzy Sets and Systems*, vol. 8, no. 3, pp. 225–233, 1982. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. L. Puri and D. A. Ralescu, “Differentials of fuzzy functions,”
*Journal of Mathematical Analysis and Applications*, vol. 91, no. 2, pp. 552–558, 1983. View at Google Scholar · View at Scopus - O. Kaleva, “Fuzzy differential equations,”
*Fuzzy Sets and Systems*, vol. 24, no. 3, pp. 301–317, 1987. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - O. Kaleva, “The Cauchy problem for fuzzy differential equations,”
*Fuzzy Sets and Systems*, vol. 35, no. 3, pp. 389–396, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Seikkala, “On the fuzzy initial value problem,”
*Fuzzy Sets and Systems*, vol. 24, no. 3, pp. 319–330, 1987. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Abbasbandy and T. Allah Viranloo, “Numerical solution of fuzzy differential equation,”
*Mathematical & Computational Applications*, vol. 7, no. 1, pp. 41–52, 2002. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Abbasbandy, T. A. Viranloo, Ó. López-Pouso, and J. J. Nieto, “Numerical methods for fuzzy differential inclusions,”
*Computers & Mathematics with Applications*, vol. 48, no. 10-11, pp. 1633–1641, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - O. Solaymani Fard and T. AliAbdoli Bidgoli, “The Nystrom method for hybrid fuzzy differential equation IVPs,”
*Journal of King Saud University Science*, vol. 23, pp. 371–379, 2011. View at Google Scholar - T. Allahviranloo, N. Ahmady, and E. Ahmady, “Numerical solution of fuzzy differential equations by predictor-corrector method,”
*Information Sciences*, vol. 177, no. 7, pp. 1633–1647, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. Allahviranloo, E. Ahmady, and N. Ahmady, “$n$th-order fuzzy linear differential equations,”
*Information Sciences*, vol. 178, no. 5, pp. 1309–1324, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Barkhordari Ahmadi and M. Khezerloo, “Fuzzy bivariate Chebyshev method for solving fuzzy Volterra-Fredholm integral equations,”
*International Journal of Industrial Mathematics*, vol. 3, pp. 67–77, 2011. View at Google Scholar - A. M. Bertone, R. M. Jafelice, L. C. de Barros, and R. C. Bassanezi, “On fuzzy solutions for partial differential equations,”
*Fuzzy Sets and Systems*, vol. 219, pp. 68–80, 2013. View at Google Scholar - H. Kim and R. Sakthivel, “Numerical solution of hybrid fuzzy differential equations using improved predictor-corrector method,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 17, no. 10, pp. 3788–3794, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. T. Malinowski, “Existence theorems for solutions to random fuzzy differential equations,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 73, no. 6, pp. 1515–1532, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Mazandarani and A. V. Kamyad, “Modified fractional Euler method for solving fuzzy fractional initial value problem,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 18, no. 1, pp. 12–21, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Ghazanfari and A. Shakerami, “Numerical solutions of fuzzy differential equations by extended Runge-Kutta-like formulae of order 4,”
*Fuzzy Sets and Systems*, vol. 189, pp. 74–91, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Salahshour, T. Allahviranloo, S. Abbasbandy, and D. Baleanu, “Existence and uniqueness results for fractional differential equations with uncertainty,”
*Advances in Difference Equations*, vol. 2012, 112 pages, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - S. Salahshour, T. Allahviranloo, and S. Abbasbandy, “Solving fuzzy fractional differential equations by fuzzy Laplace transforms,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 17, no. 3, pp. 1372–1381, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Ziari, R. Ezzati, and S. Abbasbandy, “Numerical solution of linear fuzzy fredholm integral equations of the second kind using fuzzy Haar wavelet,” in
*Advances in Computational Intelligence*, vol. 299 of*Communications in Computer and Information Science*, pp. 79–89, Springer, New York, NY, USA, 2012. View at Google Scholar - Ö. Akın, T. Khaniyev, O. Oruç, and I. B. Türksen, “An algorithm for the solution of second order fuzzy initial value problems,”
*Expert Systems with Applications*, vol. 40, pp. 953–957, 2013. View at Google Scholar - D. N. Georgiou, J. J. Nieto, and R. Rodríguez-López, “Initial value problems for higher-order fuzzy differential equations,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 63, no. 4, pp. 587–600, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Khastan, F. Bahrami, and K. Ivaz, “New results on multiple solutions for $n$th-order fuzzy differential equations under generalized differentiability,”
*Boundary Value Problems*, vol. 2009, Article ID 395714, 13 pages, 2009. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. Rabiei, F. Ismail, S. Norazak, and N. Abasi, “Construction of improved Runge-Kutta Nystrom method for solving second-order ordinary differential equations,”
*World Applied Sciences Journal*. In press. - F. Rabiei and F. Ismail, “Third-order improved Runge-Kutta method for solving ordinary differential equation,”
*International Journal of Applied Physics and Mathematics*, vol. 1, pp. 191–194, 2011. View at Google Scholar - F. Rabiei and F. Ismail, “Fifth-order improved Runge-Kutta method for solving ordinary differential equation,”
*Australian Journal of Basic and Applied Sciencs*, vol. 6, pp. 97–105, 2012. View at Google Scholar - R. Goetschel, Jr. and W. Voxman, “Elementary fuzzy calculus,”
*Fuzzy Sets and Systems*, vol. 18, no. 1, pp. 31–43, 1986. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H.-J. Zimmermann,
*Fuzzy Set Theory and Its Applications*, Kluwer Academic Publishers, Boston, Mass, USA, 2nd edition, 1992. View at MathSciNet - B. Bede and S. G. Gal, “Generalizations of the differentiability of fuzzy-number-valued functions with applications to fuzzy differential equations,”
*Fuzzy Sets and Systems*, vol. 151, no. 3, pp. 581–599, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Chalco-Cano and H. Román-Flores, “On new solutions of fuzzy differential equations,”
*Chaos, Solitons & Fractals*, vol. 38, no. 1, pp. 112–119, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. M. Franco, “A class of explicit two-step hybrid methods for second-order IVPs,”
*Journal of Computational and Applied Mathematics*, vol. 187, no. 1, pp. 41–57, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. J. van der Houwen and B. P. Sommeijer, “Explicit Runge-Kutta (-Nyström) methods with reduced phase errors for computing oscillating solutions,”
*SIAM Journal on Numerical Analysis*, vol. 24, no. 3, pp. 595–617, 1987. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet