Research Article  Open Access
Sankar Prasad Mondal, Susmita Roy, Biswajit Das, "Numerical Solution of FirstOrder Linear Differential Equations in Fuzzy Environment by RungeKuttaFehlberg Method and Its Application", International Journal of Differential Equations, vol. 2016, Article ID 8150497, 14 pages, 2016. https://doi.org/10.1155/2016/8150497
Numerical Solution of FirstOrder Linear Differential Equations in Fuzzy Environment by RungeKuttaFehlberg Method and Its Application
Abstract
The numerical algorithm for solving “firstorder linear differential equation in fuzzy environment” is discussed. A scheme, namely, “RungeKuttaFehlberg method,” is described in detail for solving the said differential equation. The numerical solutions are compared with (i)gH and (ii)gH differential (exact solutions concepts) system. The method is also followed by complete error analysis. The method is illustrated by solving an example and an application.
1. Introduction
Fuzzy Differential Equation. In modeling of real natural phenomena, differential equations play an important role in many areas of discipline, exemplary in economics, biomathematics, science, and engineering. Many experts in such areas widely use differential equations in order to make some problems under study more comprehensible. In many cases, information about the physical phenomena related is always immanent with uncertainty.
Today, the study of differential equations with uncertainty is instantaneously growing as a new area in fuzzy analysis. The terms such as “fuzzy differential equation” and “fuzzy differential inclusion” are used interchangeably in mention to differential equations with fuzzy initial values or fuzzy boundary values or even differential equations dealing with functions on the space of fuzzy numbers. In the year 1987, the term “fuzzy differential equation (FDE)” was introduced by Kandel and Byatt [1]. There are different approaches to discuss the FDEs: (i) the Hukuhara derivative of a fuzzy number valued function is used, (ii) Hüllermeier [2] and Diamond and Watson [3–5] suggested a different formulation for the fuzzy initial value problems (FIVP) based on a family of differential inclusions, (iii) in [6, 7], Bede et al. defined generalized differentiability of the fuzzy number valued functions and studied FDE, and (iv) applying a parametric representation of fuzzy numbers, Chen et al. [8] established a new definition for the differentiation of a fuzzy valued function and used it in FDE.
Solution of Fuzzy Differential Equation by Numerical Techniques. Numerical methods are the methods by which we can find the solution of differential equation where the exact solution is critical to find. There exist various numerical methods for solving differential equation such as Setia et al. [9], Liu [10], and Setia et al. [11]. Our aim is to find the numerical techniques by which the solution of a linear or nonlinear firstorder fuzzy differential equation comes easily and the solution is very close to the exact solution. There exist many techniques of numerical methods for finding the solution of fuzzy differential equation. Authors applied the method in certain types of fuzzy differential equation which shows that their techniques are best fit for that particular problem. The first paper on fuzzy differential equation and numerical analysis was published in 1999 by Ma et al. [12]. Allahviranloo et al. [13] apply the twostep method on fuzzy differential equations. Allahviranloo et al. [14] find the numerical solution by using predictorcorrector method. Allahviranloo et al. [15] find an algorithm for finding the solution thorder fuzzy linear differential equations using numerical techniques. Pederson and Sambandham in [16] use characterization theorem on hybrid fuzzy initial value problem. A soft computing technique, namely, artificial neural network, is implicated for solving FDE by Effati and Pakdaman [17]. Duraisamy and Usha [18] used modified Euler’s method. The extension principle method was compared by Euler’s method in Saberi Najafi et al. [19] article. Rostami et al. [20] find a numerical algorithm for solving nonlinear fuzzy differential equations. Moghadam and Dahaghin [21] apply twostep methods for numerical solution of FDE. Batiha in [22] finds an iterative solution of multispecies predatorprey model by variational iteration method. Ahmad and Hasan [23] proposed a new fuzzy version of Euler’s method for solving differential equations with fuzzy initial value. Nirmala and Chenthur Pandian [24] give an idea for improving the numerical result on FDE. Shafiee et al. [25] use predictorcorrector method for nonlinear fuzzy Volterra integral equations. Comparison results on some numerical techniques on firstorder fuzzy differential equation are illustrated by Ghanbari [26]. The use of variational iteration method for solving thorder fuzzy differential equations is shown by Jafari et al. [27]. Tapaswini and Chakraverty [28] discuss a new approach to fuzzy initial value problem by Improved Euler method. The solution of FIVP is compared by Least Square method and Adomian Decomposition method by Ahmed and Fadhel [29]. Solution of differential equation by Euler’s method using fuzzy concept is developed by Saikia [30]. Ezzati et al. [31] find the numerical solution of VolterraFredholm integral equations with the help of inverse and direct discrete fuzzy transforms and collocation technique. The Adomian method is applied on secondorder FDE by Wang and Guo [32]. Fard [33] uses iterative scheme to find the solution of generalized system of linear FDE, whereas Block method is used by Mehrkanoon et al. [34]. Asady and Alavi [35] apply a numerical method for solving thorder linear fuzzy differential equation.
Solution of Fuzzy Differential Equation by RungeKutta Method. RungeKutta method is well known for finding the approximate or numerical solution. In the last decade RungeKutta method is applied in fuzzy differential equation for finding the numerical solution. The researchers are giving various types of view to apply these methods. Someone changes the order and someone applies different types on FDE, a comparison of another method to RungeKutta method. The details of published work done in RungeKutta method are summarized below.
Numerical Solution of Fuzzy Differential Equations by RungeKutta method of order three is developed by Duraisamy and Usha [36]. Solution techniques for fourthorder RungeKutta method with higher order derivative approximations are developed by Nirmala and Chenthur Pandian [37]. RungeKutta method of order five is developed by Jayakumar et al. [38]. The techniques extended RungeKuttalike formulae of order four are developed by Ghazanfari and Shakerami [39]. Thirdorder RungeKutta method is developed by Kanagarajan and Sambath [40]. RungeKuttaFehlberg method for hybrid fuzzy differential equation is solved by Jayakumar and Kanagarajan [41]. A different approach followed by RungeKutta method is applied by Akbarzadeh Ghanaie and Mohseni Moghadam [42]. “Numerical Solution of Fuzzy IVP with Trapezoidal and Triangular Fuzzy Numbers by Using FifthOrder RungeKutta Method” is solved by Ghanbari [43]. “New MultiStep RungeKutta Method for Solving Fuzzy Differential Equations” is solved by Nirmala and Chenthur Pandian [44]. “Numerical Solution of Fuzzy Hybrid Differential Equation by Third Order RungeKutta Nystrom Method” is solved by Saveetha and Chenthur Pandian [45]. A new approach to solve fuzzy differential equation by using thirdorder RungeKutta method is developed by Deshmukh [46]. RungeKutta method of order four is developed by Duraisamy and Usha [47] and order five is developed by Jayakumar and Kanagarajan [48].
Application of Fuzzy Differential Equation. Fuzzy differential equations play a significant role in the fields of biology, engineering, and physics as well as among other fields of science, for example, in population models [49], civil engineering [50], bioinformatics and computational biology [51], quantum optics and gravity [52], modeling hydraulic [53], HIV model [54], decay model [55], predatorprey model [56], population dynamics model [57], friction model [58], growth model [59], bacteria culture model [60], bank account and drug concentration problem [61], barometric pressure problem [62], concentration problem [63], weight loss and oil production model [64], arm race model [65], vibration of mass [66], and fractional predatorprey equation [67].
Novelties. Although some developments are done, some new interest and new work have been done by ourselves which are mentioned below:(i)Differential equation is solved in fuzzy environment by numerical techniques; that is, coefficients and initial condition both are taken as fuzzy number of a differential equation and solved by numerical techniques.(ii)The numerical solution is compared with the exact solution ((i)gH and (ii)gH both cases).(iii)RungeKuttaFehlberg method for solving fuzzy differential equation is used.(iv)For application purpose a mixture problem is considered.(v)The solutions are found using different step length for better accuracy of the result.(vi)The necessary algorithm for numerical solution is given.
Structure of the Paper. The paper is organized as follows: in Preliminary Concepts, the preliminary concepts and basic concepts on fuzzy number and fuzzy derivative are given. The method for finding the exact solution is discussed in Exact Solution of Fuzzy Differential Equation. In Numerical Solution of Fuzzy Differential Equation we proposed RungeKuttaFehlberg method in fuzzy environment. The convergence of the said method and algorithm for finding the numerical results are also discussed in this section. Numerical Example shows a numerical example. In Application an important application, namely, mixture problem, is illustrated in fuzzy environment. Finally conclusions and future research scope of this paper are drawn in last section, Conclusion.
2. Preliminary Concepts
Definition 1 (fuzzy set). A fuzzy set is defined by . In the pair the first element belongs to the classical set and the second element belongs to the interval , called membership function.
Definition 2 (cut of a fuzzy set). The level set (or interval of confidence at level or cut) of the fuzzy set of is a crisp set that contains all the elements of that have membership values in greater than or equal to ; that is, .
Definition 3 (fuzzy number). The basic definition of fuzzy number is as follows [30]: if we denote the set of all real numbers by and the set of all fuzzy numbers on is indicated by then a fuzzy number is mapping such that , which satisfies the following four properties:(i) is upper semicontinuous.(ii) is a fuzzy convex; that is, for all .(iii) is normal; that is, for which .(iv) is support of and the closure of () is compact.
Definition 4 (parametric form of fuzzy number [31]). A fuzzy number is represented by an ordered pair of functions , , that satisfy the following condition:(1) is a bounded left continuous nondecreasing function for any .(2) is a bounded left continuous nonincreasing function for any .(3) for any .
Note. If , then is a crisp number.
Definition 5 (generalized Hukuhara difference [20]). The generalized Hukuhara difference of two fuzzy numbers is defined as follows:Consider ; then and .
Here the parametric representation of a fuzzy valued function is expressed by , .
Definition 6 (generalized Hukuhara derivative for first order [20]). The generalized Hukuhara derivative of a fuzzy valued function at is defined as If satisfying (2) exists, we say that is generalized Hukuhara differentiable at .
Also we say that is (i)gH differentiable at if and is (ii)gH differentiable at if
Definition 7 (see [6]). For arbitrary and , the quantityis the distance between fuzzy numbers and .
Definition 8 (triangular fuzzy number). A triangular fuzzy number (TFN) denoted by is defined as where the membership function
Definition 9 (cut of a fuzzy set ). The cut of is given by
Definition 10 (fuzzy ordinary differential equation (FODE)). Consider a simple 1storder linear ordinary differential equation as follows:The above ordinary differential equation is called fuzzy ordinary differential equation if any one of the following three cases holds:(i)Only is a fuzzy number (TypeI).(ii)Only is a fuzzy number (TypeII).(iii)Both and are fuzzy numbers (TypeIII).
3. Exact Solution of Fuzzy Differential Equation
Consider the fuzzy initial value problemwhere is a continuous mapping from into and with level setsWe write and , .
Because of we have the following.
When is (i)gH differentiableWhen is (ii)gH differentiablewhere, by using extension principle, we have the membership function So fuzzy number is . From this it follows that where
Note. (1) Both cases ((i)gH and (ii)gH) can be applied to a FDE for finding exact solution.
(2) After taking cut of the given FDE, it transforms to system of ordinary differential equation.
4. Numerical Solution of Fuzzy Differential Equation
4.1. RungeKuttaFehlberg Method for Ordinary (Crisp) Differential Equation
Consider the initial value problem .
The RungeKuttaFehlberg method (denoted as RKF45) is one way to try to resolve this problem.
The problem is to solve the initial value problem in above equation by means of RungeKutta methods of order 4 and order 5.
First we need some definitions: Then an approximation to the solution of initial value problem is made using RungeKutta method of order 4:A better value for the solution is determined using a RungeKutta method of order 5:The optimal step size can be determined by multiplying the scalar times the step size . The scalar iswhere is the specified error control tolerance.
Note that RK4 requires 4 function evaluations and RK5 requires 6 evaluations, that is, 10 for RK4 and RK5. Fehlberg devised a method to get RK4 and RK5 results using only 6 function evaluations by using some of values in both methods.
4.2. RungeKuttaFehlberg Method for Solving Fuzzy Differential Equations
Let be the exact solution and let be the approximated solution of the fuzzy initial value problem.
Let .
Throughout this argument, the value of is fixed. Then the exact and approximated solution at are, respectively, denoted byThe grid points at which the solution is calculated are .
Then we obtainedwherewhere
4.3. Convergence of Fuzzy RungeKuttaFehlberg Method
The solution is calculated by grid points at and .
Therefore, we haveClearly, and converge to and , respectively, whenever ; that is,
Proof. Before we go to the main proof we need to know some results.
Lemma 11. Let the sequence of numbers satisfyfor some given positive constants and . Then
Lemma 12. Let the sequence of numbers and satisfyfor some given positive constants and , and denoteThen where and .
Let and be obtained by substituting in (21) and (23); that is,The domain where and are defined is as
Theorem 13. Let and belong to and let the partial derivative of and be bounded over . Then for arbitrary fixed , the numerical solution of (9), converges to the exact solution .
Proof (see [46]). By using Taylor’s theorem we getwhere .
Now if we denotethen the above two expressions converted toHence we can writewhere for and is a bound for the partial derivative of and .
Therefore we can writewhere .
In particular,Since , we haveThus, if , we get and , which completes the proof.
4.4. Algorithm for Finding the Numerical Solution
Step 1. “Function to be supplied” “Function to be supplied”
Step 2. Read , limit.
Step 3. For limit , + − − − −
Step 4.
Step 5.
Step 6. . Write , .
Step 7. Next
Step 8. End.
5. Numerical Example
Example. Solve with initial condition . Then find the solution at .
Solution. For (i)gH differentiable case the exact solution isand for (ii)gH differentiable case the exact solution is
Remark 14. From Figure 1 and Table 1 we conclude that the lower exact solution ((i)gH case) is approximately equal to the numerical solution when we take the step length (for is nearly equal), whereas the lower exact solution ((ii)gH case) is approximately equal to the numerical solution when we take the step length .

Remark 15. From Figure 2 and Table 1 we conclude that the upper exact solution ((i)gH case) is approximately equal to the numerical solution when we take the step length (for is nearly equal), whereas the upper exact solution ((ii)gH case) is approximately equal to the numerical solution when we take the step length .
6. Application
Problem. A tank initially contains gals of brine which has dissolved in lbs of salt. Coming into the tank at gals/min is brine with concentration lbs salt/gals and the well stirred mixture leaves at the rate gals/min. Let lbs be the salt in the tank at any time . Then , with , if the initial condition is being modeled as fuzzy numbers and . Find solution at .
Solution. For (i)gH differentiable case the exact solution isFor (ii)gH differentiable case the exact solution is
Remark 16. From Figure 3 and Table 2 we conclude that the lower exact solution ((i)gH case) is approximately equal to the numerical solution when we take the step length (for is nearly equal), whereas the lower exact solution ((ii)gH case) is not equal to any numerical solution.

Remark 17. From Figure 4 and Table 2 we conclude that the upper exact solution ((i)gH case) is approximately equal to the numerical solution when we take the step length of and . For it is nearly equal, whereas the upper exact solution ((ii)gH case) is not equal to any numerical solution.
7. Conclusion
In this paper we applied RungeKuttaFehlberg method for finding the numerical solution of firstorder linear differential equation in fuzzy environment. The numerical solution is compared with the exact solution ((i)gH and (ii)gH both cases). The results presented in the contribution show that RungeKuttaFehlberg method is a powerful mathematical tool for solving firstorder linear differential equation in fuzzy environment. The convergence of RungeKuttaFehlberg method has been discussed. The process method is applied to a mechanical problem in fuzzy environment which shows that it is a promising method to solve the said types of differential equation. In the future we can apply these methods for solving higher order linear and nonlinear differential equation in fuzzy environment.
Competing Interests
The authors declare that there are no competing interests.
References
 A. Kandel and W. J. Byatt, “Fuzzy differential equations,” in Proceedings of the International Conference Cybernetics and Society, pp. 1213–1216, Tokyo, Japan, 1978. View at: Google Scholar  MathSciNet
 E. Hüllermeier, “An approach to modelling and simulation of uncertain dynamical systems,” International Journal of Uncertainty, Fuzziness and KnowledgeBased Systems, vol. 5, pp. 117–137, 1997. View at: Publisher Site  Google Scholar
 P. Diamond, “Timedependent differential inclusions, cocycle attractors and fuzzy differential equations,” IEEE Transactions on Fuzzy Systems, vol. 7, no. 6, pp. 734–740, 1999. View at: Publisher Site  Google Scholar
 P. Diamond, “Stability and periodicity in fuzzy differential equations,” IEEE Transactions on Fuzzy Systems, vol. 8, no. 5, pp. 583–590, 2000. View at: Publisher Site  Google Scholar
 P. Diamond and P. Watson, “Regularity of solution sets for differential inclusions quasiconcave in a parameter,” Applied Mathematics Letters, vol. 13, no. 1, pp. 31–35, 2000. View at: Publisher Site  Google Scholar  MathSciNet
 B. Bede and S. G. Gal, “Generalizations of the differentiability of fuzzynumbervalued functions with applications to fuzzy differential equations,” Fuzzy Sets and Systems, vol. 151, no. 3, pp. 581–599, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 B. Bede, I. J. Rudas, and A. L. Bencsik, “First order linear fuzzy differential equations under generalized differentiability,” Information Sciences, vol. 177, no. 7, pp. 1648–1662, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 M. Chen, C. Wu, X. Xue, and G. Liu, “On fuzzy boundary value problems,” Information Sciences, vol. 178, no. 7, pp. 1877–1892, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 A. Setia, Y. Liu, and A. S. Vatsala, “Numerical solution of fractional integrodifferential equations with nonlocal boundary conditions,” Journal of Fractional Calculus and Applications, vol. 5, no. 2, pp. 155–165, 2014. View at: Google Scholar
 Y. Liu, “An efficient algorithm for solving differential equations to facilitate modeling and simulation of aerospace systems,” SAE Technical Paper 2015012402, 2015. View at: Publisher Site  Google Scholar
 A. Setia, Y. Liu, and A. S. Vatsala, “Solution of linear fractional Fredholm integrodifferential equation by using second kind Chebyshev wavelet,” in Proceedings of the 11th International Conference on Information Technology (ITNG '14), pp. 465–469, Las Vegas, Nev, USA, April 2014. View at: Publisher Site  Google Scholar
 M. Ma, M. Friedman, and A. Kandel, “Numerical solutions of fuzzy differential equations,” Fuzzy Sets and Systems, vol. 105, no. 1, pp. 133–138, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 T. Allahviranloo, N. Ahmady, and E. Ahmady, “Two step method for fuzzy differential equations,” International Mathematical Forum, vol. 1, no. 17–20, pp. 823–832, 2006. View at: Google Scholar  MathSciNet
 T. Allahviranloo, N. Ahmady, and E. Ahmady, “Numerical solution of fuzzy differential equations by predictorcorrector method,” Information Sciences, vol. 177, no. 7, pp. 1633–1647, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 T. Allahviranloo, E. Ahmady, and N. Ahmady, “Nthorder fuzzy linear differential equations,” Information Sciences, vol. 178, pp. 1309–1324, 2008. View at: Google Scholar
 S. Pederson and M. Sambandham, “Numerical solution of hybrid fuzzy differential equation IVPs by a characterization theorem,” Information Sciences, vol. 179, no. 3, pp. 319–328, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 S. Effati and M. Pakdaman, “Artificial neural network approach for solving fuzzy differential equations,” Information Sciences, vol. 180, no. 8, pp. 1434–1457, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 C. Duraisamy and B. Usha, “Another approach to solution of fuzzy differential equations by modified Euler's method,” in Proceedings of the International Conference on Communication and Computational Intelligence (INCOCCI '10), pp. 52–55, Kongu Engineering College, Erode, India, December 2010. View at: Google Scholar
 H. Saberi Najafi, F. Ramezani Sasemasi, S. Sabouri Roudkoli, and S. Fazeli Nodehi, “Comparison of two methods for solving fuzzy differential equations based on Euler method and Zadehs extension,” The Journal of Mathematics and Computer Science, vol. 2, no. 2, pp. 295–306, 2011. View at: Google Scholar
 M. Rostami, M. Kianpour, and E. bashardoust, “A numerical algorithm for solving nonlinear fuzzy differential equations,” The Journal of Mathematics and Computer Science, vol. 2, no. 4, pp. 667–671, 2011. View at: Google Scholar
 M. Sh. Dahaghin and M. M. Moghadam, “Analysis of a twostep method for numerical solution of fuzzy ordinary differential equations,” Italian Journal of Pure and Applied Mathematics, vol. 27, pp. 333–340, 2010. View at: Google Scholar
 K. Batiha, “Numerical solutions of the multispecies predatorprey model by variational iteration method,” Journal of Computer Science, vol. 3, no. 7, pp. 523–527, 2007. View at: Publisher Site  Google Scholar
 M. Z. Ahmad and M. K. Hasan, “A new fuzzy version of Euler's method for solving differential equations with fuzzy initial values,” Sains Malaysiana, vol. 40, no. 6, pp. 651–657, 2011. View at: Google Scholar  Zentralblatt MATH
 V. Nirmala and S. Chenthur Pandian, “A method to improve the numerical solution of fuzzy initial value problems,” International Journal of Mathematical Archive, vol. 3, no. 3, pp. 1231–1240, 2012. View at: Google Scholar
 M. Shafiee, S. Abbasbandy, and T. Allahviranloo, “Predictorcorrector method for nonlinear fuzzy volterra integral equations,” Australian Journal of Basic and Applied Sciences, vol. 5, no. 12, pp. 2865–2874, 2011. View at: Google Scholar
 M. Ghanbari, “Solution of the first order linear fuzzy differential equations by some reliable methods,” Journal of Fuzzy Set Valued Analysis, vol. 2012, Article ID jfsva00126, 20 pages, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 H. Jafari, M. Saeidy, and D. Baleanu, “The variational iteration method for solving nth order fuzzy differential equations,” Central European Journal of Physics, vol. 10, no. 1, pp. 76–85, 2012. View at: Publisher Site  Google Scholar
 S. Tapaswini and S. Chakraverty, “A new approach to fuzzy initial value problem by improved Euler method,” Fuzzy Information and Engineering, vol. 4, no. 3, pp. 293–312, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. A. Ahmed and F. S. Fadhel, “Solution of fuzzy initial value problems using Least square and Adomian decomposition method,” Journal of Kufa for Mathematics and Computer, vol. 1, no. 4, pp. 49–60, 2011. View at: Google Scholar
 R. K. Saikia, “Solution of differential equation by Euler's method using fuzzy concept,” International Journal of Computer Technology & Applications, vol. 3, no. 1, pp. 226–230, 2012. View at: Google Scholar
 R. Ezzati, F. Mokhtari, and M. Maghasedi, “Numerical solution of volterrafredholm integral equations with the help of inverse and direct discrete fuzzy transforms and collocation technique,” International Journal of Industrial Mathematics, vol. 4, no. 3, pp. 221–229, 2012. View at: Google Scholar
 L. Wang and S. Guo, “Adomian method for secondorder fuzzy differential equation,” World Academy of Science, Engineering and Technology, vol. 76, pp. 979–982, 2011. View at: Google Scholar
 O. S. Fard, “An iterative scheme for the solution of generalized system of linear fuzzy differential equations,” World Applied Sciences Journal, vol. 7, no. 12, pp. 1597–1604, 2009. View at: Google Scholar
 S. Mehrkanoon, M. Suleiman, and Z. A. Majid, “Block method for numerical solution of fuzzy differential equations,” International Mathematical Forum, vol. 4, no. 46, pp. 2269–2280, 2009. View at: Google Scholar
 B. Asady and M. Alavi, “Numerical method for solving nthorder linear fuzzy differential equations,” Advances in Fuzzy Sets and Systems, vol. 6, no. 2, pp. 95–106, 2010. View at: Google Scholar  MathSciNet
 C. Duraisamy and B. Usha, “Numerical solution of fuzzy differential equations by RungeKutta method of order three,” International Journal of Advanced Scientific and Technical Research, vol. 1, no. 1, pp. 1–13, 2011. View at: Google Scholar
 V. Nirmala and S. Chenthur Pandian, “Numerical solution of fuzzy differential equation by fourth order RungeKutta method with higher order derivative approximations,” European Journal of Scientific Research, vol. 62, no. 2, pp. 198–206, 2011. View at: Google Scholar
 T. Jayakumar, D. Maheskumar, and K. Kanagarajan, “Numerical solution to fuzzy differential equations by Runge Kutta method of order five,” Applied Mathematical Sciences, vol. 6, no. 57–60, pp. 2989–3002, 2012. View at: Google Scholar  MathSciNet
 B. Ghazanfari and A. Shakerami, “Numerical solutions of fuzzy differential equations by extended RungeKuttalike formulae of order 4,” Fuzzy Sets and Systems, vol. 189, pp. 74–91, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 K. Kanagarajan and M. Sambath, “Numerical solution of fuzzy differential equations by third order RungeKutta method,” International Journal of Applied Mathematics and Computation, vol. 2, no. 4, pp. 1–8, 2010. View at: Google Scholar
 T. Jayakumar and K. Kanagarajan, “Numerical solution for hybrid fuzzy differential equations by RugeKutta Fehlberg method,” International Journal of Mathematical Analysis, vol. 6, no. 5356, pp. 2619–2632, 2012. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 Z. Akbarzadeh Ghanaie and M. Mohseni Moghadam, “Solving fuzzy differential equations by RungeKutta method,” The Journal of Mathematics and Computer Science, vol. 2, no. 2, pp. 208–221, 2011. View at: Google Scholar
 M. Ghanbari, “Solution of the first order linear fuzzy differential equations by some reliable methods,” Journal of Fuzzy Set Valued Analysis, vol. 2012, Article ID jfsva00126, 20 pages, 2012. View at: Publisher Site  Google Scholar
 V. Nirmala and S. Chenthur Pandian, “New multistep RungeKutta method for solving fuzzy differential equations,” Mathematical Theory and Modeling, vol. 1, no. 3, pp. 16–22, 2011. View at: Google Scholar
 N. Saveetha and S. Chenthur Pandian, “Numerical solution of fuzzy hybrid differential equation by third order Runge Kutta Nystrom method,” Mathematical Theory and Modeling, vol. 2, no. 4, pp. 8–17, 2012. View at: Google Scholar
 N. M. Deshmukh, “A new approach to solve Fuzzy differential equation by using third order RungeKutta method,” ISSN No20315063, Vol. 1, Issue 3, pp. 1–4, 2011. View at: Google Scholar
 C. Duraisamy and B. Usha, “Numerical solution of fuzzy differential equations by RungeKutta method of order four,” European Journal of Scientific Research, vol. 67, no. 2, pp. 324–337, 2012. View at: Google Scholar
 T. Jayakumar and K. Kanagarajan, “Numerical solution for hybrid fuzzy system by Runge Kutta method of order five,” Applied Mathematical Sciences, vol. 6, no. 72, pp. 3591–3606, 2012. View at: Google Scholar
 J. J. Buckley, T. Feuring, and Y. Hayashi, “Linear systems of first order ordinary differential equations: fuzzy initial conditions,” Soft Computing, vol. 6, no. 6, pp. 415–421, 2002. View at: Publisher Site  Google Scholar
 M. Oberguggenberger and S. Pittschmann, “Differential equations with fuzzy parameters,” Mathematical and Computer Modelling of Dynamical Systems, vol. 5, no. 3, pp. 181–202, 1999. View at: Publisher Site  Google Scholar
 J. Casasnovas and F. Rossell, “Averaging fuzzy biopolymers,” Fuzzy Sets and Systems, vol. 152, no. 1, pp. 139–158, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 M. S. El Naschie, “From experimental quantum optics to quantum gravity via a fuzzy Kähler manifold,” Chaos, Solitons & Fractals, vol. 25, no. 5, pp. 969–977, 2005. View at: Publisher Site  Google Scholar
 A. Bencsik, B. Bede, J. Tar, and J. Fodor, “Fuzzy differential equations in modeling hydraulic differential servo cylinders,” in Proceedings of the 3rd RomanianHungarian Joint Symposium on Applied Computational Intelligence (SACI '06), Timisoara, Romania, 2006. View at: Google Scholar
 H. Zarei, A. V. Kamyad, and A. A. Heydari, “Fuzzy modeling and control of HIV infection,” Computational and Mathematical Methods in Medicine, vol. 2012, Article ID 893474, 17 pages, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 G. L. Diniz, J. F. R. Fernandes, J. F. C. A. Meyer, and L. C. Barros, “A fuzzy Cauchy problem modelling the decay of the biochemical oxygen demand in water,” in Proceedings of the Joint 9th IFSA World Congress and 20th NAFIPS International Conference, IEEE, Vancouver, Canada, July 2001. View at: Google Scholar
 M. Z. Ahmad and B. De Baets, A PredatorPrey Model with Fuzzy Initial Populations, IFSAEUSFLAT, 2009.
 L. C. Barros, R. C. Bassanezi, and P. A. Tonelli, “Fuzzy modelling in population dynamics,” Ecological Modelling, vol. 128, no. 1, pp. 27–33, 2000. View at: Publisher Site  Google Scholar
 B. Bede, I. J. Rudas, and J. Fodor, “Friction model by using fuzzy differential equations,” in Proceedings of the 12th International Fuzzy Systems Association World Congress (IFSA '07), Cancun, Mexico, June 2007, vol. 4529 of Lecture Notes in Computer Science, pp. 23–32, Springer, Berlin, Germany, 2007. View at: Google Scholar
 S. P. Mondal, S. Banerjee, and T. K. Roy, “First order linear homogeneous ordinary differential equation in fuzzy environment,” International Journal of Pure and Applied Sciences and Technology, vol. 14, no. 1, pp. 16–26, 2013. View at: Google Scholar
 S. P. Mondal and T. K. Roy, “First order linear non homogeneous ordinary differential equation in fuzzy environment,” Mathematical theory and Modeling, vol. 3, no. 1, pp. 85–95, 2013. View at: Google Scholar
 S. P. Mondal and T. K. Roy, “First order linear homogeneous ordinary differential equation in fuzzy environment based on laplace transform,” Journal of Fuzzy Set Valued Analysis, vol. 2013, Article ID jfsva00174, 18 pages, 2013. View at: Publisher Site  Google Scholar
 S. P. Mondal and T. K. Roy, “First order linear homogeneous fuzzy ordinary differential equation based on lagrange multiplier method,” Journal of Soft Computing and Applications, vol. 2013, Article ID jsca00032, 17 pages, 2013. View at: Publisher Site  Google Scholar
 S. P. Mondal and T. K. Roy, “First order linear non homogeneous ordinary differential equation in fuzzy environment based on Laplace transform,” Journal of Mathematics and Computer Science, vol. 3, no. 6, pp. 1533–1564, 2013. View at: Google Scholar
 S. Prasad Mondal and T. Kumar Roy, “First order homogeneous ordinary differential equation with initial value as triangular intuitionistic fuzzy number,” Journal of Uncertainty in Mathematics Science, vol. 2014, Article ID jums00003, 17 pages, 2014. View at: Publisher Site  Google Scholar
 S. P. Mondal and T. K. Roy, “System of differential equation with initial value as triangular intuitionistic fuzzy number and its application,” International Journal of Applied and Computational Mathematics, vol. 1, no. 3, pp. 449–474, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 S. Chakraverty and S. Tapaswini, “Standerd orthogonal polynomial based solution of fuzzy differential equation,” The International Journal of Fuzzy Computation and Modelling, vol. 1, no. 1, 2014. View at: Google Scholar
 S. Tapaswini and S. Chakraverty, “Numerical solution of fuzzy arbitrary order predatorprey equations,” Applications & Applied Mathematics, vol. 8, no. 2, pp. 647–672, 2013. View at: Google Scholar
Copyright
Copyright © 2016 Sankar Prasad Mondal 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.