ANew Approach for Solving Dual Fuzzy Nonlinear Equations Using Broyden ’ s and Newton ’ s Methods

where F : Rn → Rn is a nonlinear mapping. The value of variable x is called a solution or root of the nonlinear equations. The most widest approach to solve such nonlinear systems is Newton’s initiative [1], yet it required to compute the Jacobian matrix in every iteration. However, in some cases, the coefficients of the nonlinear systems are given in fuzzy numbers instead of crisp numbers. Therefore, there is a need to explore some possible numerical methods for solving fuzzy nonlinear equations. It is vital to mention that the basic concept of fuzzy numbers were first presented in [2–4], and the famous application of fuzzy number arithmetic is systems of nonlinear equations in which its coefficients are given as fuzzy numbers [5, 6]. Moreover, the standard analytical technique presented by [7, 8] cannot be suitable for handing the fuzzy nonlinear equations such as


Introduction
Solving systems of nonlinear equations is becoming more essential state in analysis and handling complex problems in many research areas (e.g.robotics, radiative transfer, chemistry, economics, etc.).Consider the nonlinear systems where F : R n → R n is a nonlinear mapping.The value of variable x is called a solution or root of the nonlinear equations.The most widest approach to solve such nonlinear systems is Newton's initiative [1], yet it required to compute the Jacobian matrix in every iteration.However, in some cases, the coefficients of the nonlinear systems are given in fuzzy numbers instead of crisp numbers.Therefore, there is a need to explore some possible numerical methods for solving fuzzy nonlinear equations.It is vital to mention that the basic concept of fuzzy numbers were first presented in [2][3][4], and the famous application of fuzzy number arithmetic is systems of nonlinear equations in which its coefficients are given as fuzzy numbers [5,6].Moreover, the standard analytical technique presented by [7,8] cannot be suitable for handing the fuzzy nonlinear equations such as where, a, b, c, d, e, f , g, h, i are fuzzy numbers.In general, we consider these equations as To tackle these situations, some numerical methods have been introduced [8][9][10][11][12].For example, [11] applied Newton's method while [9] employed Broyden's method and [10] uses steepest descent method to solve fuzzy nonlinear equations, respectively.Nevertheless, the weakness of Newton's method arises from the need to compute and invert the Jacobian Matrix in every iteration.It worth to mention that [12] has extended the approach of [11] to solve dual fuzzy nonlinear systems.Nevertheless, Advances in Fuzzy Systems their approach required to compute and store the Jacobian matrix in every iteration.In this paper, a new approach via Newton's and Broyden's methods is proposed to solve dual fuzzy nonlinear equations.The anticipation has been to reduce the computational burden of the Jacobian matrix in every iteration.This paper is arranged as follows: in the next section, we present brief overview and some basic definitions of the fuzzy nonlinear equations, and description of our method is given in Section 3. Section 4 presents the alternative appoach for solving fuzzy nonlinear systems.Numerical results are reported in Section 5, and finally conclusion is given in Section 6.

Preliminaries
This section presents some vital definitions of fuzzy numbers.
Definition 1.A fuzzy number is a set like u : R → I = [0, 1] which satisfies the following conditions [13]: (1) u is upper semicontinuous, (2) u(x) = 0 outside some interval [c, d], (3) there are real numbers a, b such that c ≤ a ≤ b ≤ d and The set all these fuzzy numbers is denoted by E. An equivalent parametric is as also given in [14].

The Classical Broyden's Method
It is well known that it is not always feasible to compute the full elements of the Jacobian matrix of the given nonlinear function or it may be very expensive, we often have to approximate the Jacobian matrix by some other approaches; the famous method of doing so is quasi-Newton's method [15] and Broyden's method in particular.The basic idea underlining this method has been to reduce the evaluation cost of the Jacobian matrix.Moreover, in some situations analytic derivatives could not be done precisely, or not available to obtain, and the promising method designed to embark upon this situation is Broyden's scheme.
Broyden's method is an iterative procedure that generates a sequence of points {x k } from a given initial guess x 0 via the following form: where B k is an approximation to the Jacobian which can be updated at each iteration using a rank-one matrix for k = 0, 1, 2 and so forth.The appealing feature of this method is that it requires only one function evaluation per iteration.
If F(x) is an affine function, and there exist holds.By (4), the updated matrix B k+1 is chosen in such a way that it satisfies the secant equation, that is, where From ( 5), the update formulae for the Broyden matrix B k is given as [15] Hence, the number of scalar function evaluations is reduced from n 2 + n to n.In this study, we use classical approach, that is, by using direct inverse B −1 instead of approximate B −1 used by [9].In the following, we state the convergence theorems of the Broyden's method, we referred to the proof in [16].
Theorem 3. Let the standard assumptions hold and let r ∈ (0, 1) be given.Then there are δ and δ B such that if x 0 ∈ B(δ) and E 0 2 < δ B the Broyden sequence for the data (F, x 0 , B 0 ) exists and x k → x * q-linearly with q-factor at most r.
Theorem 4. Let the standard assumptions hold.Then there are δ and δ B such that if x 0 ∈ B(δ) and E 0 2 < δ B the Broyden sequence for the data (F, x 0 , B 0 ) exists and x k → x * q-superlinearly.

Classical Broyden's Method for Solving Dual Fuzzy Nonlinear Equations
Generally, there exists no inverse of any given fuzzy number, say x ∈ E such that where y ∈ E [17].In fact, for any nonscrip fuzzy number x ∈ E [12], it is true that Here, we consider the dual fuzzy nonlinear system as where all parameters are fuzzy numbers.The basic idea of this section is to obtain a solution for the above dual fuzzy nonlinear equations, whose parametric version is given as Assume that x = (λ, λ) is the solution to the above fuzzy nonlinear equation, then From ( 11), we have Letting F(λ, λ, r) = c(r) and F(λ, λ, r) = c(r) in ( 12), we obtained To solve (10), an initial guess is required and then generating a sequence of points {x k } k≥0 .Now, we can describe the algorithm for our proposed approach (Broyden's Newton) as follows.
Step 1. Transform the dual fuzzy nonlinear equations into parametric form.
Step 2. Determine the initial guess x 0 by solving the parametric equations for r = 0 and r = 1.
Step 5. Compute the next point via Broyden's method.
Step 6. Repeat steps from 3 to 5 and continue with the next k until the stopping criteria ( ≤ 10 −5 ) are satisfied.

Numerical Results
In this section, we consider two problems to illustrate the performance of our approach for solving dual fuzzy nonlinear equations.The computations are performed in Matlab 7.0 using double precision computer.The benchmark problems are from [12].
Problem 1 (see [12]).Consider (2, 1, 1) Without loss of generality, assume x is positive, then we have the parametric equation as [12] ( Equivlantly, We obtained the initial point by letting r = 0 in (18) as follows: For r = 1, we have We consider the initial guess x 0 = (0.8, 0.9, 0.9).By implementing Algorithm 5, the solution was obtained in three iterations with maximum error less than 10 −5 .The performance profile of the positive solution for r ∈ [0, 1] is given in Figure 1.
Problem 2 (see [12]).Consider (6, 2, 2) Without loss of generality, let x be positive; hence the parametric form of ( 22) is given as [12] (4 + 2r)x 2 (r) + (1 + r)x(r) (23) In other way, (23) can be written as and the Jacobian is (25) Hence, the Jacobian inverse is given as follows: To obtain the initial values, we let r = 0 and r = 1 in (23), respectively, therefore Thus, from (27), we have x(0) = 0.4343, x(0) = 0.5307 and x(1) = x(1) = 0.5000.Therefore, initial guess x 0 = (x(0), x(1), x(0)).From our own observation, the x 0 is very close to the solution.Therefore, in order to illustrate the performance of our approach, we consider x 0 = (0.4,0.5, 0.6).Via Algorithm 5 with x 0 = (0.4,0.5, 0.6) by repeating 3 to 5 until stopping criterion is satisfied.After three iterations, the solution was obtained with maximum error less than 10 −5 .It has been shown in [12] that the negative root of this dual fuzzy nonlinear systems does not exist, that is why we consider only the positive solutions.We present the details of the solution for all r ∈ [0, 1] in Figure 2. Figures 1 and  2 illustrates the efficiency of our approach on solving dual fuzzy nonlinear equations.Our approach converges in three iterations with maximum error less than 10 −5 .

Conclusion
A new approach for solving dual fuzzy nonlinear equations was presented.The approach reduces computational cost of the Jacobian matrix in every iteration.The fuzzy nonlinear equations are transformed into parametric and then solved via Newton's for initial iteration and Broyden's method for rest of the iterations.Numerical experiment shows that in all the benchmark problems, our approach is very promising.Finally, we can claim that our scheme is a good candidate for solving dual fuzzy nonlinear equations.

Figure 2 :
Figure 2: Positive solution of new method for Problem 2.
Figure 1: Positive solution of the new method for Problem 1.