New Approximation Methods Based on Fuzzy Transform for Solving SODEs: II

: In this research, three approximation methods are used in the new generalized uniform fuzzy partition to solve the system of differential equations (SODEs) based on fuzzy transform (FzT). New representations of basic functions are proposed based on the new types of a uniform fuzzy partition and a subnormal generating function. The main properties of a new uniform fuzzy partition are examined. Further, the simpler form of the fuzzy transform is given alongside some of its fundamental results. New theorems and lemmas are proved. In accordance with the three conventional numerical methods: Trapezoidal rule (one step) and Adams Moulton method (two and three step modiﬁcations), new iterative methods (NIM) based on the fuzzy transform are proposed. These new fuzzy approximation methods yield more accurate results in comparison with the above-mentioned conventional methods


Introduction
Differential equation is particularly useful for different areas of applied sciences and engineering.Many differential equations have no closed form solutions.Thus, many researchers are developing approximation methods for solving differential equations, for example [1][2][3].In this paper, we continue the study of approximation methods based on FzT to solutions of differential equations.
The core idea of FzT is a fuzzy partition of a universe into fuzzy subsets.The first fuzzy partition of FzT with the Ruspini condition was introduced by [4] and was extensively investigated by [5].This condition implies normality of the fuzzy partition.In addition, the fuzzy partition with the generalized Ruspini condition (fuzzy r-partition) was introduced by [6].This fuzzy partition was achieved by replacing the partition of unity by fuzzy r-partition.This type of partition was used by [6,7] for smoothing or filtering data based on the inverse FzT.Further, a generalized fuzzy partition appeared in connection with the notion of the FzT, where FzT components are polynomials of degree m [8].By [9], different types of fuzzy partitions are taken into consideration such as B-splines, Shepard kernels, Bernstein basis polynomials and Favard-Szasz-Mirakjan operators.Later, the higher degree FzT based on B-splines was proposed [10] to improve the quality of the function approximation of two variables.
A generalized fuzzy partition was implicitly introduced by [11] with the purpose of meeting the requirements of image compression.In addition, a generalized fuzzy partition can also be considered in connection with radial membership functions [12].Further, necessary and sufficient conditions for modeling the generalized fuzzy partition was provided by [13].Recently, a new representation formula for basic functions of FzT and a new fuzzy numerical method based on block pulse functions for numerical solution of integral equations were presented by [14].The approximation method based on the FzT with Shepard-type basic functions for linear Fredholm integral equations was discussed by [15].New representations of the generalized uniform fuzzy partitions with the normal case to obtain better approximation solutions for solving Cauchy problems were presented by [16].
FzT is a soft computing method developed by Perfilieva [5] that has many applications, for example, in differential and integral equations.FzT for solving ordinary Cauchy problems with one variable was initiated by [4].The generalization of the Euler method has been discussed by [17] for solving ordinary Cauchy problems.The author has applied this technique to reef growth and sea level variations models.Further, FzT has been generalized from the case of constant components to the case of polynomial components by [8].Later, the first and second degree FzT based mid-point rule for solving the Cauchy problem and the uncertain initial value problem have been proposed by [18].Furthermore, an algorithm to obtain the approximate solutions of second order initial value problems was constructed by [19].From this idea, FzT for numerical solutions of two point boundary value problems was proposed by [20].
FzT of two variables based on finite differences method was used by [21] for solving a type of partial differential equations with Dirichlet boundary conditions and initial conditions.In addition, the first degree FzT of two variables was introduced by [22].By [23], the partial derivatives using the first FzT were approximated and modification of the Canny edge detector was proposed.Furthermore, the uniform stability result for the vibrations of a telegraph equation using FzT of two variables was proposed by [24].The composition of inverse and direct discrete FzT method was extended to numerical solution of Fredholm integral equations and Volterra Fredholm integral equations [25].The general form of the higher order FzT was constructed by [26] for solving differential and integral equations using any arbitrary basis functions.The FzT has investigated for solving the Volterra population growth model using the approximation for the Caputo derivative [27].A new numerical method based FzT was demonstrated to solve a class of delay differential equations by means of the Picard-like numerical scheme [28].FzT was considered to approximate the solution of boundary value problems by minimizing the integral squared error in 2-norm [29].In [30], the dynamical properties of a two neuron system with respect to FzT and a single delay have been investigated.The conditions under which quasi-consensus in a multi-agent system with sampled data based on FzT were proposed by [31].
NIM was proposed to solve nonlinear functional equations and the existence of solution for nonlinear Volterra integral equations [2].At the same time, NIM was introduced for solving nonlinear equations by using a different decomposition technique [32].From this conception, NIM was considered in terms up to fourth-order in Taylor series for solving nonlinear equations [33].Sufficiency conditions have been presented for convergence of the NIM [34].A new predictor-corrector approach was developed based on NIM for fractional differential equations [35].Classical methods are modified by [3] to derive numerous formulas for solving the differential equations.
The motivation of the proposed study comes from [16,36,37].In [16], new fuzzy numerical methods to solve the Cauchy problem was considered and the authors showed that the error can be reduced by FzT and NIM with respect to new generalized uniform fuzzy partitions, namely power of the triangular and raised cosine generalized uniform fuzzy partitions, where generating functions are normal (see also [37] for another approach).In addition, two basic approximation methods, modified Euler method and Trapezoidal rule, with help from FzT for solving SODEs are analyzed in detail by [36].For this purpose, more generally, new generalized uniform fuzzy partitions are proposed in this study, where a generating function is not normal.
The membership functions in underlying fuzzy partitions are often called basic functions.There has been a growing interest in investigating the properties of fuzzy partitions.However, the problem arises on how one can effectively construct the basic function of fuzzy partitions.
In this paper, new representations of basic functions are proposed.This is achieved by introducing new generalized uniform fuzzy partitions, where a generating function is not normal.Further, new fuzzy numerical methods based on NIM and FzT for solving SODEs are introduced and discussed.In particular, we consider functions of two variables with initial conditions.In accordance with the existing methods, Trapezoidal rule and Adams Moulton are improved using FzT and NIM.The methods are combined with one-step, two-step and three-step.As an application, all these methods are used to solve a general model of the dynamical system, i.e., Lotka-Volterra equation with derivatives and with variable coefficients.Furthermore, numerical examples are presented.It is observed that the new fuzzy numerical methods yield more accurate results than classical Trapezoidal rule and classical Adams Moulton methods (2 and 3-step).
The paper is organized as follows.The main part of the paper is Sections 3 and 4, which provides new representations for basic functions of FzT, followed by the modified one step, 2-step and 3-step based on NIM and FzT method with respect to new representations formulas for generalized uniform fuzzy partition of FzT.In Section 5, numerical examples are discussed.Finally, conclusions are given in Section 6.

Basic Concepts
In this section, we give some definitions and introduce the necessary notation following [38], which will be used throughout the paper.Throughout this section, we deal with an interval [a, b] ⊂ R of real numbers.

Definition 1. (generalized uniform fuzzy partition
We say that the fuzzy sets A i : [a, b] → [0, 1] constitute a generalized fuzzy partition of [a, b] if the following conditions are fulfilled: 1. (positivity and locality)- Fuzzy sets A 1 , . . ., A n are called basic functions.It is important to remark that by conditions of locality and continuity, b a A i (t)dt > 0. A generalized uniform fuzzy partition of [a, b] is defined for equidistant nodes, i.e., for all i = 1, . . ., n − 1, t i = t i+1 + h, where h = (b − a) / (n − 1) and two additional properties are satisfied,

4.
A then the fuzzy partition is called h-uniform generalized fuzzy partition.
The following definition recalls the concept of the generalized fuzzy partition which can be easily extended to the interval [a, b].We assume that [a, b] is partitioned by A 1 , . . . ,A n , according to Definition 1.

Definition 3.
A h-uniform generalized fuzzy partition of interval [a, b], determined by the triplet (K, h, a), can be defined using generating function K (Definition 2).Then, basic functions of a h-uniform generalized fuzzy partition are shifted copies of K defined by for all i = 1, . . ., n.The parameter h is called the bandwidth or the shift of the fuzzy partition and the nodes t i = a + ih are called the central point of the fuzzy sets A 1 , . . ., A n .

Remark 1. A h-uniform fuzzy partition is called Ruspini if the following condition
holds for any t ∈ [t i , t i+1 ].This condition is often called Ruspini condition.

New Iterative Method
NIM have proposed by [2] for solving linear and nonlinear functional equations of the form where f 1 is a known function and N a non linear operator.Solutions obtained by this method are in the form of rapidly converging infinite series which can be effectively approximated by calculating only the first few terms.In this method non linear operator N is decomposed as In [2], the authors were defined the recurrence relation: Hence u satisfies the functional (2).

New Representations for Basic Functions of FzT
Let us recall the basic facts of an FzT of a continuous real function f as presented by [5,17].The first step in the definition of the FzT of f involves the selection of a fuzzy partition of the domain [a, b] by a finite number n ≥ 2 of fuzzy sets B k (t), k = 1, . . ., n.In those papers, five axioms specified B k (t), k = 1, . . ., n, in the fuzzy partition: normality, locality, continuity, unimodality (monotonicity) and orthogonality (Ruspini condition).A fuzzy partition is called uniform if the fuzzy sets B k (t), k = 2, . . ., n − 1, are shifted copies of symmetrized B 1 (more details can be found in [17]).The membership functions B k (t), k = 1, . . ., n, in a fuzzy partition are called basic functions.Later, a generalized fuzzy partition appeared in connection with the notion of a higher-degree FzT [8].Furthermore, summarize both these notions in [38].Three axioms specify B k (t), k = 1, . . ., n, in the fuzzy partition: positivity and locality, continuity and covering.Recently, the different conditions for generalized uniform fuzzy partitions was proposed [13,38] while another approach was demonstrated by [37] where a function can be reconstructed from its F-transform components.In the following, we modify the definition h-uniform generalized fuzzy partition.

Generalized Uniform Fuzzy Partitions with the Generalized Normal Case
Let us recall the h-uniform generalized fuzzy partition of real line can be defined using generating function K.Then, basic functions of the h-uniform generalized fuzzy partition are shifted copies of K. On the basis of Definition 1 can be also defined using a generating function λβK(t) where β = 1/K(0), K(0) = 0, β > 0 and λ > 0 (in general, not necessarily satisfying normal and Ruspini condition) which is that K(t) assumed to be even, continuous and K (t) > 0 if t ∈ (−1, 1).Therefore, we will modify the basic functions of the h-uniform generalized fuzzy partition so that they are shifted copies of λβK defined by The parameter h is bandwidth of the fuzzy partition and t 0 + kh = t k .The concept of the h-uniform generalized fuzzy partition can be easily extended to the interval [a, b] as follows.Definition 4. Let t 1 < . . .< t n be fixed nodes within [a, b] ⊂ R, such that t 1 = a, t n = b and n ≥ 2. We consider nodes t 1 , . . ., t n are equidistant, with distance (shift) h = (b − a) / (n − 1).A system of fuzzy sets B 1 , . . ., B n [a, b] → [0, 1] be a generalized uniform fuzzy partitions of [a, b] if it is defined by where In this subsection, we present the main principles of FzT with respect to new representations of h-uniform generalized fuzzy partition.Further, we will show that FzT components with respect to new representations of h-uniform generalized fuzzy partition can be simplified and approximated of an original function, say f .Definition 5. Let f be a continuous function on [a, b] and B k (t), k = 1, . . ., n, be h-uniform generalized fuzzy partition of [a, b], n ≥ 2. A vector of real numbers F[ f ] = (F 1 , F 2 , . . . ,F n ) given by for k = 1, . . ., n is called the direct FzT of f with respect to B k .
In the following, we will simplify the representation (6).
) and according to Definition 4, fuzzy sets B k , k = 1, . . ., n, n ≥ 2, be a h-uniform generalized fuzzy partition of [a, b] with a generating function K, then representation (6) of direct FzT can be simplified for k = 1, . . ., n as follows Proof.By Definition 4, we get for k = 1, . . ., n, t 0 = t 1 , t n+1 = t n , and substituting u = t−t k h and then substituting t = s/h.Thus, we get and its corresponding results with representation (6).
If λ > 0, the Lemma 1 still hold by choosing suitable constant λ, satisfying λ = 1/ 1 −1 βK(t)dt , where 1 −1 βK(t)dt > 0. So, we will restrict ourselves to h-uniform generalized fuzzy partition with 0 In the following, we will simplify the above given expressions for the coefficients (6).This fact is very important for applications which are more flexible and consequently easier to use.

Corollary 1.
Let the assumptions of Lemma 2 be fulfilled and for k = 1, . . ., n, where interval [a, b] is partitioned by the h-uniform generalized fuzzy partition B 1 , . . . ,B n .
Proof.Let k ∈ {1, . . . ,n} and consider set of fuzzy sets B k (t) be the h-uniform generalized fuzzy partition of [a, b] defined by (5).Using the proof of Lemma 2, we get where 0 h is bandwidth of the fuzzy partition and t k = a + (k − 1)h and then its corresponding in the expression (6).

Corollary 2. Let the conditions of Lemma 3 be fulfilled. Then for each
Proof.According to [5,39] The following theorem estimates the difference between the original function and its direct FzT with respect to the h-uniform generalized fuzzy partition.
and the conditions of Lemma 2 be fulfilled.Then for k = 1, . . ., n Proof.By locality condition of definition of h-uniform generalized fuzzy partition, Corollary 1, Lemma 1, and according to [17], using the trapezoid formula with nodes t k−1 , t k , t k+1 to the numerical computation of the integral, we get for k = 1, . . ., n and 0 ] and the conditions of Lemma 2 be fulfilled.Let moreover, f be Lipschitz continuous with respect to t, i.e., there exists a constant L ∈ R, such that for all t ∈ [a, b] and t, t ∈ R, Then for k = 1, . . ., n where Proof.By the assumption f has continuous second order derivatives on [a, b] and is Lipschitz continuous with respect to t.Therefore, using the trapezoid rule and let us choose a value of k where is called the inverse FzT of f .
The following lemma estimates the difference between the original function and its inverse FzT.
Lemma 4. Let the assumptions of Theorem 1 and let f (t) be the inverse FzT of f with respect to the fuzzy partition of [a, b] is given by Definition 4.Then, the following estimation holds for t ∈ [a, b] and k = 1, . . ., n where . Thus, for any ε > 0 there exist n ε ∈ N and B 1 , . . ., B n ε be the h-uniform generalized fuzzy partition of [a, b] defined by (5).Then, the following estimations hold f (t) Proof.From the proof of Lemma 4 and then using Lemma (3) in the sense that for all k = 1, . . ., n, Remark 3.According to Definition (4), it is easy to see that the inverse FzT f (t k ) = F k for all k = 1, . . ., n.
On the basis of Definition 4, necessary steps of a new method to construct generalized uniform fuzzy partitions of [−1, 1] for solve case K is not normal in the following.

1.
Select the generating function K which is assumed to be even, continuous and Specify the value β = 1/K(0), where K(0) = 0 to get the normal generating function K and then If conditions β > 0 and λ > 0 holds, then construct generalized uniform fuzzy partitions of [−1, 1] by λβK (t).
Example 1.Let K : R → [0, 1] be defined by One can see in The following remark is for modified Trapezoidal rule based on FzT and NIM to solve SODEs.

Remark 4.
In view of Equation (9), Important property of the direct FzT as well as inverse FzT is their linearity, namely, given

New Fuzzy Numerical Methods for Solving SODEs
Consider the initial value problem (IVP) for the SODEs: where Lipschitz condition on D in the variable x (y), then the initial-value problem ( 16) has a unique solution x(t) (y(t)) for a ≤ t ≤ b.In many cases, the problem ( 16) cannot be solved analytically so that numerical solutions are required.In [16], new representations of basic function based on the FzT are constructed for solving generalized Cauchy problems with help of NIM, FzT and classical methods (one-step, two-step and three-step) have presented while Euler method and Mid-point rule, based on FzT to solve Cauchy problem proposed by [17,18].Further, NIM has been proposed for solving ODEs and delay differential equations [3].Moreover, Adams-Bashforth methods and Adams-Moulton methods are noted as two families of multistep methods in literature.Multistep methods refer to using several previous values from the previous steps.The Adams-Bashforth methods were presented by John Couch Adams to solve a differential equation modelling capillary action due to Francis Bashforth and it follows that the Adams-Moulton method was developed improved multistep methods for solving ballistic equations by Forest Ray Moulton.In particular, the Adams-Moulton method is similar to the Adams-Bashforth method and the Adams-Moulton method was used Newton's method to solve the implicit equation.
Clearly, Adams-Bashforth methods are explicit methods and the Adams-Moulton methods are implicit methods, for example, see ( [40], p. 111).Necessary steps of construction of the generalized uniform fuzzy partitions can be summarized as follows.

1.
Specify the number n of components and compute the step h = (t n − t 1 ) / (n − 1).If we want to obtain as best approximation of f as possible, then n should be large.

2.
Construct the nodes t 1 < . . .< t n , where Select the shape of basic functions.This is achieved by selecting the shape of generating function.

4.
Construct a h-uniform generalized fuzzy partition of [t 1 , t n ] by new representations of basic functions are defined by Definition 4.
To begin the derivation of a modified Trapezoidal rule (1-step) and Adams Moulton method (2 and 3-step), integrate (16) Consider the following integral t k g (s, x(s), y (s)) ds.
However, we cannot integrate f (s, x(s), y (s)) and g (s, x(s), y (s)) without knowing x(s) and y(s).So, the above integral (18) can be approximated by the following approach where f 1 and g 1 are the approximation of f and g on the interval [t k , t k+1 ].Choosing different f 1 (g 1 ) leads to different schemes.In particular, we choose f 1 (g 1 ) which contributes to the one, two and three-step methods based on FzT.Later, modification of these methods based on FzT and NIM.
In this section, we present three new schemes to solve SODEs ( 16) that use the F-transform and NIM and suppose that the functions f and g on [a, b] are sufficiently smooth.The first scheme uses 1-step method while the second one uses the 2-step method and the last uses the 3-step method.

Numerical Scheme I: Modified Trapezoidal Rule Based on FzT and NIM for SODEs
According to necessary steps of construction of the generalized uniform fuzzy partitions in Section 4, we contributed to approximation methods of SODEs ( 16) by scheme provides formulas for the FzT components, X k (Y k ), k = 2, . . ., n − 1, of the unknown function x (t) (y (t)) with respect to choose some of the h-uniform generalized fuzzy partition, B 1 , . . ., B n , of interval [a, b] with parameter h to approximate solution of SODEs (16).As initial step, choose the number n ≥ 2 and compute h = (b − a) / (n − 1), then construct the h-uniform generalized fuzzy partition of [a, b] using Definition 4. Let X 1 = x 1 and Y 1 = y 1 .In the following, we apply the FzT and NIM to the SODEs (16) for obtaining the numerical Scheme I, where k = 1, . . ., n − 1.
First, let f 1 (g 1 ) in the Equation ( 19) is chosen as where and B k represents the generalized uniform fuzzy partitions that are defined by Definition 4.Then, substituting (20) into (19) for k = 1, . . ., n − 1 By Remark 4 in the interval [t k , t k+1 ], we have Hence, the one step method based on FzT for ( 17) is derived as follows, where k = 1, . . ., n − 1.
where F k and G k are defined by (21).
This method computes the approximate coordinates [X 1 , . . ., X n ] and [Y 1 , . . ., Y n ] of the FzT for the functions x(t) and y(t).The problem with the previous scheme (22) is that the unknown quantities F k+1 and G k+1 which means that X k+1 (Y k+1 ) appears on both sides and an implicit method.Therefore, one solution to this problem would be to use an explicit method such as another fuzzy approach, namely Scheme I.For this purpose, the scheme ( 22) is of the form and can be solved by NIM (2), where The three term approximation of the NIM (3) gives the following formulas for solving SODEs ( 16): Hence, the three term approximate solution is which leads to the following formulas. where In the sequel, the approximate solution of SODEs ( 16) can be obtained using the inverse FzT as follows:

Numerical Scheme II: Modified 2-Step Adams Moulton Method Based on FzT and NIM for SODEs
The Scheme I uses 1-step method for solving SODEs (16).In this subsection, we improve 2-step Adams Moulton method using FzT and NIM for solving SODEs (16).Let us recall that the modified 2-step Adams Moulton method proposed by [16].From this idea, the modified 2-step Adams Moulton method can be extended to approximate the solution of ( 16) by necessary steps of construction of the generalized uniform fuzzy partitions in Section 4. It is worth noting that three terms of NIM were used in [16], while four terms of NIM are used in this study.Let FzT components, X k (Y k ), k = 2, . . ., n − 1, of the unknown function x (t) (y (t)) with respect to choose some of the h-uniform generalized fuzzy partition (5) and let X 1 = x 1 , Y 1 = y 1 , X 2 = x 2 and Y 2 = y 2 if possible; otherwise, we can compute FzT components X 2 and Y 2 from the numerical Scheme I.In the following, we apply the F-transform and NIM to the SODEs (16) for obtaining the numerical Scheme II, where k = 2, . . ., n − 1.First, if f 1 in the Equation ( 19) is approximated by where 26) into ( 19), then for k = 1, . . ., n − 1 Similarity, Thus, the two step method based FzT for ( 17) is given for k = 1, . . ., n − 1 as where The problem with the previous scheme (27) is that the unknown quantities F k+1 and G k+1 .Therefore, one solution to this problem would be to use an explicit method.For this purpose, the scheme ( 27) is of the form and can be solved by NIM (2), where The four term approximation of the NIM (3) gives the following formulas for solving SODEs (16): Hence, the four term approximate solution is which leads to the following formulas. where Then, obtain the desired approximation for x and y by the inverse FzT (25) applied to [X 1 , . . . ,X n ] and [Y 1 , . . . ,Y n ].

Numerical Scheme III: Modified 3-Step Adams Moulton Method Based on FzT and NIM for SODEs
In this subsection, we improve 3-step Adams Moulton method using FzT and NIM for solving SODEs (16).The modified 3-step Adams Moulton method proposed by [16] for solving Cauchy problems.From this idea, we can propose to approximate the solution of ( 16) by NIM and FzT components, X k (Y k ), k = 2, . . ., n − 1, of the unknown function x(t) (y(t)) with respect to choose some of the h-uniform generalized fuzzy partition (see Definition 4) and Y 3 = y 3 if possible; otherwise, we can compute FzT components X 2 , Y 2 , X 3 and Y 3 from the numerical Scheme I. Now, we apply the F -transform and NIM to the SODEs ( 16) and obtain the following numerical Scheme III for k = 3, . . ., n − 1: According to steps of deriving Equation ( 27) and then steps of NIM in previous Subsection 4.2, we get the four term approximation of the NIM as follows.
Hence, the four term approximate solution is which leads to the following formulas. where In the sequel, the inverse FzT (25) approximates the solution x(t) (y(t)) of the problem (16).

Error Analysis of Numerical Scheme I for SODEs
In this subsection, we present error analysis for numerical Scheme I and consider the Formula ( 23).If x(t k ) = x k and y(t k ) = y k denote the exact solution and X k , Y k denote the numerical solution.Then, substituting the exact solution in the Formula (23), we get x * * k+1 = x * k+1 + hF e * k+1 /2, y * * k+1 = y * k+1 + hG e * k+1 /2, where and the truncation error Tx k and Ty k of the scheme I are given by Rearranging ( 23), we get 0 = Let e k+1 = X k+1 − x k+1 and d k+1 = Y k+1 − y k+1 , then subtracting ( 33) from (32), we get , Similarly to Lemma 8 and Theorem 2 by [16], we have the following results.
Lemma 5. Let f , g are assumed to be sufficiently smooth functions of its arguments on [a, b] and be Lipschitz continuous with respect to x and y, i.e., there exists a constant L ∈ R, such that for all t ∈ [a, b] and x, x , y, y ∈ R, , G e k , G e * * k+1 are determined by Formula (31) and M 2 , M 3 are upper bound for f and g respectively on [a, b].The technique of error analysis for rest schemes can be obtained analogously to numerical Scheme I.

Applications
A general model for the dynamical system may be written as dx dt = xg(x, y), dy dt = yh(x, y), where g and h are arbitrary functions of the prey and predator species whose populations are x(t) and y(t) at time t.However, the following problem of Lotka-Volterra equation with derivatives and with variable coefficients α (t) , β (t) , δ (t) , γ (t) as functions of time t have not yet been solved by any fuzzy numerical method.The new differential equations are represented by a non-autonomous ordinary differential equation system.The model, incorporating the above functions is as follows [41][42][43]: Two examples are discussed in order to prove that the results obtained by Scheme I (23), II (28) and III (29) for the numerical solution of the model (36).
Example 3. Consider the problem of Lotka-Volterra-prey-predator model (36) with The exact solution for these coefficients is proposed by [41,42].
The results are listed in Tables 2-6 by the proposed fuzzy approximation methods with respect to case K C 201 2 is defined by Example 1.The proposed fuzzy approximation methods are generated by Algorithms A1-A3 (Appendix A).The mean square error (MSE) defined as MSE = 1 n ( Y k − y(t k ) 2 ) 2 .This is an easily computable quantity for a particular sample.From the numerical tests, the results are summarized as follows: 1.
In view of Tables 2-5, a comparison between the three new proposed schemes (( 23), ( 28) and ( 29)) and the classical Trapezoidal rule (1-step), classical 2-step Adams Moulton Method and classical 3-step Adams Moulton Method based on Euler method for Examples 2 and 3.

2.
Moreover, comparison of MSE for Examples 2 and 3 shown in Table 6.It is observed that the new fuzzy approximation methods yield more accurate results in comparison with the classical Trapezoidal rule (one step) and classical Adams Moulton method (two and three steps).Hence, the new fuzzy approximation methods provide alternative techniques for solving SODEs with better results.     .In view of Figures 2 and 3, the graphical results of Examples 2 and 3 show a comparison between numerical Schemes (I, II and III) and exact solutions are shown separated from each other for clarity while a comparison between three proposed fuzzy numerical methods (Schemes I, II and III) and exact solution are shown in Figure 1.Furthermore, in view of Figures 2 and 3, a comparison between the numerical results and exact solutions for h = 0.01.All the graphs are plotted using MATLAB software.

Conclusions
Three approximation methods are used the new generalized uniform fuzzy partitions for solving SODEs.In accordance with the three approximation methods for Cauchy problem by [16], Trapezoidal rule (one step) and Adams Moulton method (two and three steps) are improved using FzT and NIM.The results proved that the first approximation method converged to the exact solution.As an application, a predator-prey model is solved by using three proposed approximation methods.From the numerical results, it is observed that the new fuzzy approximation methods yield more accurate results in comparison with the classical Trapezoidal rule (one step) and classical Adams Moulton method (two and three steps).So, it is recommended to use the proposed methods to solve differential equations.
In this regard, it is well-known that FzT has a certain advantage to cope with problems affected by noise.This is because the FzT components of original and noisy functions are very similar to each other.In addition, we can reduce a higher-order differential equation into a system of first-order differential equations by relabeling the variables.Thus, the proposed methods can also be applied to a higher-order differential equation in the case of non-noisy or noisy right-hand side.From Algorithms A1-A3 (Appendix A), it is observed that the new fuzzy approximation methods are more time consuming in comparison with the considered Trapezoidal rule and the Adams Moulton methods.In the future research, we plan to give more details about running time of proposed methods.Further, we plan to solve a boundary value problem for a second order ordinary differential equation with fuzzy boundary conditions, see preliminary results in [44].
Step 3 for k = 1 to N do Steps 04-15.
Algorithm A2.Two-step algorithm for system of ODEs.
Algorithm A3.Three-step algorithm for system of ODEs.

Table 1 .
Table 1 the h-uniform generalized fuzzy partition of [a, b] determined by Definition 4. Example 1.

Table 2 .
Comparison of numerical results of x(t) for Example 2.

Table 3 .
Comparison of numerical results of y(t) for Example 2.

Table 4 .
Comparison of numerical results of x(t) for Example 3.

Table 5 .
Comparison of numerical results of y(t) for Example 3.

Table 6 .
The values of MSE for Examples 2 and 3. (a)

The Values of MSE of x(t) for SODEs Case Proposed Scheme for x(t) Classical Method for x(t)
Further, the results obtained using proposed fuzzy approximation methods for Examples 2 and 3 are shown in Figures1-3by using K C 201 2