# A Novel Differential Evolution Invasive Weed Optimization Algorithm for Solving Nonlinear Equations Systems

Authors: Yongquan Zhou; Qifang Luo; Huan Chen
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

In view of the traditional numerical method to solve the nonlinear equations exist is sensitive to initial value and the higher accuracy of defects. This paper presents an invasive weed optimization (IWO) algorithm which has population diversity with the heuristic global search of differential evolution (DE) algorithm. In the iterative process, the global exploration ability of invasive weed optimization algorithm provides effective search area for differential evolution; at the same time, the heuristic search ability of differential evolution algorithm provides a reliable guide for invasive weed optimization. Based on the test of several typical nonlinear equations and a circle packing problem, the results show that the differential evolution invasive weed optimization (DEIWO) algorithm has a higher accuracy and speed of convergence, which is an efficient and feasible algorithm for solving nonlinear systems of equations.


---

## Body

## 1. Introduction

Systems of nonlinear equations arise in many domains of practical importance such as engineering, mechanics, medicine, chemistry, and robotics. Solving such a system involves finding all the solutions (there are situations when more than one solution exists) of the polynomial equations contained in the mentioned system. The algorithms of solving nonlinear equations systems are worse than linear equations in convergence speed and ratio, especially solving nonconvex nonlinear equations. The traditional solutions of nonlinear equations include Newton-Raphson method, Quasi-Newton method, and homotopy method. Newton-Raphson method is a much more classical method, but it is sensitive to initial iteration value. In addition, not only it requires a large amount of calculation, but also sometimes it accompanied by difficulty calculation. Quasi-Newton method is to solve the difficult caused by Jacobi matrix. It has now become one of the most effective methods that solve nonlinear equations and optimization problems. While, its stability is poor and sometimes its iterative effect is not well. The basic idea of homotopy method is to start from easily solved equations and then gradually transit to the original equations and get the solution of problems. In recent years, with the rapid development of computational intelligence, computational intelligence techniques have also been used to solve nonlinear equations, such as genetic algorithm [1–3], particle swarm optimization algorithm [4], differential evolution algorithm [5], artificial fish-swarm algorithm [6], artificial bee colony algorithm [7] harmony search algorithm [8], and probabilistic-driven search (PDS) algorithm [9]. These swarm intelligent algorithm has several advantages when adopted for searching solutions for systems of nonlinear equations: their does not require of a “good” initial point to perform the search, and the search space can be bounded by lower and upper values for each decision variable. Additionally, no continuity or differentiability of the objective function is required. What can be considered as the main disadvantage of swarm intelligent algorithm in this sort of application is its relatively poor accuracy, which is caused by the coarse granularity of the search performed by the algorithm. This can, of course, be improved either by running the swarm intelligent algorithm for a larger number of iterations (although at a higher computational cost) or by postprocessing the solution produced by the swarm intelligent algorithm with a traditional numerical optimization technique.In 2006, a novel stochastic optimization model, invasive weed optimization (IWO) algorithm [10], was proposed by Mehrabian and Lucas, which is inspired from a common phenomenon in agriculture: colonization of invasive weeds. The algorithm was inspired from colonizing weeds, which is used to mimic the natural behavior of weeds in colonizing and occupying suitable places for growth and reproduction. It has the robustness, adaptation, and randomness and is simple but effective with an accurate global search ability. So far, it has been applied in many engineering fields.In this paper, a new hybrid algorithm based on the population diversity of IWO and the heuristic differential evolution (DE) is presented to solve nonlinear systems of equation. With the population diversity, which enhances the global searching ability of algorithm, and heuristic method, which improves the local mining capacity, the hybrid algorithm gets a higher optimization accuracy and faster convergence speed. The performance of the proposed approach is evaluated for several well-known benchmark problems from kinematics, chemistry, combustion, and medicine. Numerical results reveal the efficiency of the proposed approach and its flexibility to solve large-scale systems of equations.


## 2. Descriptions of Nonlinear Equation Systems

The general form of nonlinear equations systems can be described as follows:(1)f1(X)=0,f2(X)=0,⋮fn(X)=0,
where X=(x1,x2,…,xn)∈D∈Rn and D={(x1,x2,…,xn)∣ai≤xi≤bi,i=1,2,…,n}. The following fitness function is used to calculate conveniently
(2)f(X)=∑i=1nfi2(X),
where fi(X), i=1,2,…,n are nonlinear functions. Because there are many equality constraints, the system of equations usually has no solution X such that fi(X)=0, i=1,2,…,n. Thus, we find an approximate solution of simultaneous equations such that |f(X)|<ε, where ε is an arbitrary small positive number. In order to do so, we define the fitness function error:
(3)ε=|0-f(X)|.
So, the smaller the fitness error is, the higher the solution quality is. The fitness error of theoretical solutions is 0. Thus, nonlinear equation systems problems can be formulated as solving function optimization
(4)minf(X)s.t.0X=(x1,x2,…,xn)∈D.
In this way, the best value of formula (4) is the solution of (1).


## 3. Algorithm of Invasive Weed Optimization

In the basic IWO, weeds represent the feasible solutions of problems and population is the set of all weeds. A finite number of weeds are being dispread over the search area. Every weed produces new weeds depending on its fitness. The generated weeds are randomly distributed over the search space by normally distributed random numbers with a mean with a mean equal to zero equal to zero. This process continues until the maximum number of weeds is reached. Only the weeds with better fitness can survive and produce seed, others are being eliminated. The process continues until maximum iterations are reached or hopefully the weed with the best fitness is closest to the optimal solution. The process is addressed in detail as follows.Step 1 (initialize a population).
A population of initial solutions is being dispread over theD-dimensional search space with random positions.Step 2 (reproduction).
The higher the weed’s fitness is, the more seeds it produces. The formula of weeds producing seeds is(5)weedn=f-fminfmax-fmin(smax-smin)+smin,
where f is the current weed’s fitness. fmax and fmin, respectively, represent the maximum and the least fitness of the current population. smax and smin, respectively, represent the maximum and the least value of a weed.Step 3 (spatial dispersal).
The generated seeds are randomly distributed over theD-dimensional search space by normally distributed random numbers with a mean equal to zero but with a varying variance. This ensures that seeds will be randomly distributed so that they abide near to the parent plant. However, standard deviation (σ) of the random function will be reduced from a previously defined initial value (σinit) to a final value (σfinal) in every generation. In simulations, a nonlinear alteration has shown satisfactory performance given as follows:
(6)σcur=(itermax-iter)n(itermax)n(σinit-σfinal)+σfinal,
where itermax is the maximum number of iterations, σcur is the standard deviation at the present time step, and n is the nonlinear modulation index. Generally, n is set to 3.Step 4 (competitive exclusion).
After passing some iteration, the number of weeds in a colony will reach its maximum (P_MAX) by fast reproduction. At this time, each weed is allowed to produce seeds. The produced seeds are then allowed to spread over the search area. When all seeds have found their position in the search area, they are ranked together with their parents (as a colony of weeds). Next, weeds with lower fitness are eliminated to reach the maximum allowable population in a colony. In this way, weeds and seeds are ranked together and the ones with better fitness survive and are allowed to replicate. The population control mechanism also is applied to their offspring to the end of a given run, realizing competitive exclusion.


## 4. Differential Evolution Invasive Weed Optimization

Reviewing the IWO, some initial weeds are dispread over the search space randomly and then produce new individuals (seeds). We select the better plants from the population consisting of weeds and seeds. The process continues until maximum number of plants is reached. In order to speed up the optimization process, differential evolution (DE) [11] is used to cooperate with IWO so that every weed in iteration can move towards the best individual of the current iteration. In this way, the algorithm designated as DEIWO not only ensures the individual diversity by IWO, but also improves the optimization accuracy and the speed by DE (see Pseudocode 1).Pseudocode 1:The pseudocode for DEIWO algorithm.
BEGIN
Input: randomly initialize a population of DEIWO individuals.
While (the termination criterion is not satisfied)
Calculate every individual’s fitnessf(X);
Calculate the number of every individual’s offspring
weedn=f-fminfmax-fmin(smax-smin)+smin;
Updateσcur by
σcur=(itermax-iter)n(itermax)n(σinit-σfinal)+σfinal;
Generate seeds over the search space byN(0,σcur),
Add the generated seeds into the solution set;
If (the number of all weeds and seeds equal the maximum populationP_MAX)
Sort the population in descending order by fitness;
Truncate population of weeds with higher fitness until the number of set equal the
maximum population;
End If
Mutate the individual by
Vit=Xit+F×(Xbestt-Xit)+F×(Xr1t-Xr2t);
Cross the individual by
Ui,jt={Vi,jtif(randj[0,1]≤CRorj=jrand)Xi,jtotherwise, j=1,2,…,n;
Select individual betweenUit and Xit
Xit+1={Uitiff(Uit)<f(Xit)Xitotherwise;
End While
Output the best individual and its fitness;
END


### 4.1. Differential Evolution

Differential evolution (DE), introduced by Storn and Price in 1997, is one of the most prominent new generation EAs for solving real-valued optimization problems. Using only a few control parameters, DE offers exclusive advantages, such as a simple and easy-to-understand concept, ease of implementation, powerful search capability, and robustness. The main procedure of DE includes mutation, crossover, and selection. These operators are based on natural evolution principles in order to keep the population diversity, as well as to avoid premature convergence. In this way, mutation and crossover operators are used to generate new vectors. Then, the selection operator determines which vectors will survive in the next generation. This procedure is repeated until a stopping condition is reached. The basic principle of DE is given as follows.Suppose the objective function can be expressed as follow.(7)Minimizef(X).
The problem space is S={X∣X=(x1,x2,…,xi,…,xn)}, where ai≤xi≤bi, i=1,2,…,n; for each parent vector Xit a mutant vector is generated according to
(8)Vit=Xit+F×(Xbestt-Xit)+F×(Xr1t-Xr2t),
where parent vector Xit, Xr1t, Xr2t are mutually exclusive individuals. The parent vector Xbestt is the best individual of the current iteration. The scaling factor F is a positive constant which controls the amplification of the difference vector, and F∈[0,2]. After mutation, the trial vector is generated using the parent and mutated vectors as follows:
(9)Ui,jt={Vi,jtif(randj[0,1]≤CRorj=jrand)Xi,jtotherwise,j=1,2,…,n,
where CR∈[0,1] is a real-valued crossover rate constant and jrand∈[1,n] is a random integer.The individual of the next generation is generated as the follows:(10)Xit+1={Uit,iff(Uit)<f(Xit),Xit,otherwise,
where f(x) denotes the fitness function.


### 4.2. The Hybrid Strategy of DE and IWO

The hybrid strategy of DE and IWO can be summarized as follows. Every initial weed produces new seeds depending on its fitness. The generated seeds are randomly distributed over the search space by normally distributed random numbers. This process continues until the maximum number of weeds is reached. Only the weeds with better fitness can survive and produce seed, when the others are being eliminated. We can see that IWO ensures the diversity of individuals and the worse individuals also have the opportunity to reproduce, while individuals of DE move toward the best individual, using the information of distance and direction obtained from the current population. We can conclude that DE improves the convergence speed and searching precision. So, combining DE and IWO is useful to enrich the search behavior of optimization process and get high quality solutions.


## 4.1. Differential Evolution

Differential evolution (DE), introduced by Storn and Price in 1997, is one of the most prominent new generation EAs for solving real-valued optimization problems. Using only a few control parameters, DE offers exclusive advantages, such as a simple and easy-to-understand concept, ease of implementation, powerful search capability, and robustness. The main procedure of DE includes mutation, crossover, and selection. These operators are based on natural evolution principles in order to keep the population diversity, as well as to avoid premature convergence. In this way, mutation and crossover operators are used to generate new vectors. Then, the selection operator determines which vectors will survive in the next generation. This procedure is repeated until a stopping condition is reached. The basic principle of DE is given as follows.Suppose the objective function can be expressed as follow.(7)Minimizef(X).
The problem space is S={X∣X=(x1,x2,…,xi,…,xn)}, where ai≤xi≤bi, i=1,2,…,n; for each parent vector Xit a mutant vector is generated according to
(8)Vit=Xit+F×(Xbestt-Xit)+F×(Xr1t-Xr2t),
where parent vector Xit, Xr1t, Xr2t are mutually exclusive individuals. The parent vector Xbestt is the best individual of the current iteration. The scaling factor F is a positive constant which controls the amplification of the difference vector, and F∈[0,2]. After mutation, the trial vector is generated using the parent and mutated vectors as follows:
(9)Ui,jt={Vi,jtif(randj[0,1]≤CRorj=jrand)Xi,jtotherwise,j=1,2,…,n,
where CR∈[0,1] is a real-valued crossover rate constant and jrand∈[1,n] is a random integer.The individual of the next generation is generated as the follows:(10)Xit+1={Uit,iff(Uit)<f(Xit),Xit,otherwise,
where f(x) denotes the fitness function.


## 4.2. The Hybrid Strategy of DE and IWO

The hybrid strategy of DE and IWO can be summarized as follows. Every initial weed produces new seeds depending on its fitness. The generated seeds are randomly distributed over the search space by normally distributed random numbers. This process continues until the maximum number of weeds is reached. Only the weeds with better fitness can survive and produce seed, when the others are being eliminated. We can see that IWO ensures the diversity of individuals and the worse individuals also have the opportunity to reproduce, while individuals of DE move toward the best individual, using the information of distance and direction obtained from the current population. We can conclude that DE improves the convergence speed and searching precision. So, combining DE and IWO is useful to enrich the search behavior of optimization process and get high quality solutions.


## 5. Experiments and Results

This section reports several experiments and comparisons using the proposed approach. Some well-known applications are also considered in the subsequent section.


### 5.1. Testing Platform

The experimental program testing platform included processor: CPU Intel Core i3-370, frequency: 2.40 GHz, memory: 4 GB, operating system: Windows 7, and run software: MATLAB 7.6.


### 5.2. Testing Nonlinear Equation Systems

In order to test the performance of DEIWO for solving nonlinear equation systems, 8 nonlinear equation systems in the literature are used and the testing results are compared with the literature [12]. Table 2 shows the experiment parameters. The search area and the number of iterations are as shown in Table 3. Table 4 to Table 11 is the results of 8 equation systems.Example 1.
Consider the following nonlinear system:(11)f1(x1,x2)=cos(2x1)-cos(2x2)-0.4=0,f2(x1,x2)=2(x2-x1)+sin(2x2)-sin(2x1)-1.2=0.
The results obtained by applying Newton, Secant, Broyden, and Effati methods; evolutionary approach (EA); Probabilistic-driven search; IWO; and the proposed DEIWO method are presented in Table4.Example 2.
Consider the following nonlinear system:(12)f1(x1,x2)=ex1+x1x2-1=0,f2(x1,x2)=sin(x1x2)+x1+x2-1=0.
The results obtained by Effati, evolutionary approach (EA), probabilistic-driven search, IWO, and the proposed DEIWO method are given in Table5.Example 3 (interval arithmetic benchmark).
We consider one benchmark problem proposed from interval arithmetic [13]. The benchmark consists of the following system of equations:(13)0=x1-0.25428722-0.18324757x4x3x9,0=x2-0.37842197-0.16275449x1x10x6,0=x3-0.27162577-0.16955071x1x2x10,0=x4-0.19807914-0.15585316x7x1x6,0=x5-0.44166728-0.19950920x7x6x3,0=x6-0.14654113-0.18922793x8x5x10,0=x7-0.42937161-0.21180486x2x5x8,0=x8-0.07056438-0.17081208x1x7x6,0=x9-0.34504906-0.19612740x10x6x8,0=x10-0.42651102-0.21466544x4x8x1.
Parameters used by the DEIWO approach are listed in Tables2 and 3. The results obtained by evolutionary approach (EA) and the proposed DEIWO method are given in Table 6.Example 4 (neurophysiology application).
We considered the example proposed in [14], which consisted of the following equations:
(14)x12+x32=1,x22+x42=1,x5x33+x6x43=c1,x5x13+x6x23=c2,x5x1x32+x6x42x2=c3,x5x12x3+x6x22x4=c4,
where the constant ci can be randomly chosen. In our experiments, we considered ci=0, i=1,2,3,4, as in the literature [12]. In [15], this problem is used to show the limitations of Newton’s method for which the running time exponentially increases with the size of the initial intervals. We considered the following values for the parameters used by the DEIWO as given in Tables 2 and 3. Some of the solutions obtained by our approach as well as the values of the objective functions for these values are presented in Table 7.Example 5 (chemical equilibrium application).
We consider the chemical equilibrium system as given by the following [16]:
(15)x1x2+x1-3x5=0,2x1x2+x1+x2x32+R8x2-Rx5+2R10x22+R7x2x3+R9x2x4=0,2x2x32+2R5x32-8x5+R6x3+R7x2x3=0,R9x2x4+2x42-4Rx5=0,x1(x2+1)+R10x22+x2x32+R8x2+R5x32+x42-1+R6x3+R7x2x3+R9x2x4=0,
where
(16)R=10,R5=0.193,R6=0.00259740,R7=0.00344840,R8=0.0000179940,R9=0.000215540,R10=0.0000384640.
The parameters used by the DEIWO approach are presented in Tables 2 and 3. Some of the solutions obtained by the DEIWO approach for the chemical equilibrium application are depicted in Table 8.Example 6 (kinematic application).
We consider the kinematic application kin2 as introduced in [17], which describes the inverse position problem for a six-revolute-joint problem in mechanics. The equations describe a denser constraint system and are given as follows:
(17)xi2+xi+12-1=0,a1ix1x3+a2ix1x4+a3ix2x3+a4ix2x4+a5ix2x7+a6ix5x8+a7ix6x7+a8ix6x8+a9ix1+a10ix2+a11ix3+a12ix4+a13ix5+a14ix6+a15ix7+a16ix8+a17i=0,000000000001≤i≤4.
The coefficients, aki, 1≤k≤17, 1≤i≤4, are given in Table 1. The parameters used by the DEIWO approach are presented in Tables 2 and 3. Some of the solutions obtained by the DEIWO approach for the kinematic example kin2 are presented in Table 9.Table 1
Coefficientsaki for the kinematic application.
i
=
1
i
=
2
i
=
3
i
=
4
k
=
1
−0.249150680
+0.125016350
−0.635550077
+1.48947730
k
=
2
+1.609135400
−0.686607360
−0.115719920
+0.23062341
k
=
3
+0.279423430
−0.119228120
−0.666404480
+1.32810730
k
=
4
+1.434801600
−0.719940470
+0.110362110
−0.25864503
k
=
5
+0.000000000
−0.432419270
+0.290702030
+1.16517200
k
=
6
+0.400263840
+0.000000000
+1.258776700
−0.26908494
k
=
7
−0.800527680
+0.000000000
−0.629388360
+0.53816987
k
=
8
+0.000000000
−0.864838550
+0.581404060
+0.58258598
k
=
9
+0.074052388
−0.037157270
+0.195946620
−0.20816985
k
=
10
−0.083050031
+0.035436896
−1.228034200
+2.68683200
k
=
11
−0.386159610
+0.085383482
+0.000000000
−0.69910317
k
=
12
−0.755266030
+0.000000000
−0.079034221
+0.35744413
k
=
13
+0.504201680
−0.039251967
+0.026387877
+1.24991170
k
=
1
4
−1.091628700
+0.000000000
−0.057131430
+1.46773600
k
=
15
+0.000000000
−0.432419270
−1.162808100
+1.16517200
k
=
16
+0.049207290
+0.000000000
+1.258776700
+1.07633970
k
=
17
+0.049207290
+0.013873010
+2.162575000
−0.69686809Table 2
Parameters values used by the DEIWO.
G
_
SIZE
P
_
MAX
σ
final
s
max
s
min
n
F
CR
10
15
0.01
15
1
4
2
0.9Table 3
Benchmarks used in the experiments.
Benchmark
Number of variables
Variables range
Number of iterations
(1) Interval arithmetic
10
[
-
2,2
]
2000
(2) Neurophysiology application
6
[
-
10,10
]
800
(3) Chemical equilibrium application
5
[
-
10,10
]
800
(4) Kinematic application
8
[
-
10,10
]
800
(5) Combustion application
10
[
-
10,10
]
800
(6) Economics modeling application
20
[
-
10,10
]
800Table 4
Comparison of results for Example1.
Method
Solution
Functions values
Error
Newton
(
0.15,0.49
)
(
-
0.00168,0.01497
)
0.000
22156
Secant
(
0.15,0.49
)
(
-
0.00168,0.01497
)
0.000
22156
Broyden
(
0.15,0.49
)
(
-
0.00168,0.01497
)
0.000
22156
Effati
(
0.1575,0.4970
)
(
-
0.005455,0.00739
)
0.0000
84369
EA [12]
(
0.157772,0.49458
)
(
0.001264,0.000969
)
0.00000
25366
PDS [9]
(
0.156520
,
0.493376
)
(
-
0.0000005815
,
-
0.0000008892
)
0.0000003477
IWO
500 iterations
(
0.1564457783
,
0.493379536
6
)
(
0.0000
51022
,
0.000
29976
)
0.0000000
92
56
800 iterations
(
0.1565388795
,
0.493391007
2
)
(
0.0000
12829
,
-
0.0000
2800
8
)
0.0000000000
9
DEIWO
500 iterations
(
0.156520069
7
,
0.493376374
2
)
(
-
0.0000000000
,
-
0.000000000
)
0.0000000000
800 iterations
(
0.156520069
7
,
0.4933763742
)
(
-
0.000000000
,
-
0.000000000
)
0.0000000000Table 5
Comparison of results for Example2.
Method
Solution
Functions values
Error
Effati
(
0.0096,0.9976
)
(
-
0.019223,0.016776
)
0.000
65095
EA [12]
(
-
0.00138,1.0027
)
(
-
0.00276
,
-
6.37
e
-
005
)
0.00000
76216
PDS [9]
(
0.0
,
1.0
)
(
0
,
0
)
0
IWO
500 iterations
(
0.0000353012
,
1.0000971972
)
(
0.0000
70606
,
0.000
1678
)
0.00000000
33
800 iterations
(
0.000053562
8
,
0.999977272
7
)
(
0.000
10712
,
0.000
84396
)
0.0000000
18
6
DEIWO
500 iterations
(
0.0000000000
,
1
.
0000000000
)
(
0.0000000000
,
0.00000000000
)
0.0000000000
800 iterations
(
0,1
)
(
0,0
)
0Table 6
Comparison of results for interval arithmetic benchmark.
Algorithm
EA
DEIWO
Solution
Functions values
Error
Solution
Functions values
Error
Solution 1
0.04649051150.10135683570.0840577820−0.13884603090.4943905739−0.07606851630.2475819110−0.01707481560.00036675350.1481119311
0.20779592400.27697988460.18768632120.33678871140.05303913210.22237305350.18160847520.0878963860.34472003660.2784227489
0.5578712619
0.266267198585396   0.378860657321047   0.279521817014353   0.201184236934515   0.443621074720278   0.147345385930806   0.4335602098777670.076147905678193   0.347936915389233   0.426164139956603
0.0083944895882410.0022825354155090.0006069506301620.000454038870950.0016087850411080.0019199019521440.0014778742457500.0026780167873950.0019500565783850.001222532385146
0.00009717937
Solution 2
0.12248197610.18262006850.2356779803−0.03711504700.37481818560.22133113410.06978130430.0768058043−0.03121538670.1452667120
0.13185527900.19644283610.03649870690.23548901550.06757530640.07399865880.36070382920.00591829790.37674877630.2811693568
0.4739431843
0.2590708561100570.382754529706839   0.274163620083770   0.194046880920087   0.446234352737533   0.147250545938791   0.4309567896938030.070959777303814   0.348016260950456   0.428727381578941
0.0013908674071930.0016706723471030.0046702321651840.0065945286183550.0010960025611250.0018594515638650.0009818536199010.0024128010630460.0020886061068820.001450589369406
0.0000879354
Solution 3
0.06339443990.1017426933−0.1051842285−0.04770599430.41498583260.12151953210.25397771590.0843972823−0.05341329920.0880998746
0.19084366530.27678973670.37690634360.24609007020.02603377510.02560547600.17614864010.13498698510.39863956910.3383563536
0.6396343329
0.2581779734823070.379453934641213   0.282106464391719   0.203187486961515   0.440886347480303   0.151149001651968   0.4380736052048880.075807440998570   0.343828213649516   0.426228284414650
0.0002792359113930.0016751067868830.0034009084520080.0024440221469600.0045076660001450.0019122062055020.0060158253638390.0023230120026440.0021786966880970.001136405771158
0.0000920236
Solution 4
0.19398201990.01521144000.16186543450.00569858090.1904538879−0.16236040330.1864448178−0.04493027060.1675935311−0.0274959004
0.06033352800.36335147260.10974657920.91146537680.25023582290.30894605610.24289922220.11449162850.17741618960.4539962587
1.4462409770
0.2567239775508220.384592485498424   0.276941817209508   0.197489408658780   0.440994124940054   0.147051336865180   0.4311368547949640.072704222159548   0.340803111821165   0.427549623883685
0.0009788984305040.0035435466836580.0018413211436140.0041761183898890.0020837586812790.0008464904602560.0006403136778010.0051424546098790.000247320573489
0.0000663200
Solution 5
0.1169663983−0.0360324410−0.0517944631−0.18259074480.07419020560.250360462900.20430198030.01206070750.187993760800.09312965555
0.13764661610.41489827280.32335368400.38160231220.36694852620.10380356430.22506280070.05952559500.15711045160.3333260708
0.7766780224
0.2559035380827300.382537175786251   0.280411159251412   0.195875663823431   0.443780436204940   0.143389617571427   0.4298459432405990.075841793689244   0.346507927398647   0.427574364325872
0.0018712859668010.0015616894881070.0016886114282790.0046617057634130.0013350145077740.0058746754791760.0022526741607970.0025832411177900.000546906470420.000247274165109
0.0000789255
Solution 6
0.06006249220.06650344530.1163378165−0.04569937750.1649150798−0.12237710450.06665599530.07328665930.07459618230.0578573421
0.19415205260.31184931040.15532713710.24370215880.27656286840.26905055560.36288585610.00280596690.27035110710.3686104960
0.7079949322
0.2557991630284610.377253755510825   0.277018794733188   0.200978421031471   0.438492258899619   0.148958902639442   0.4284636903048680.075053838513460   0.348037273986552   0.422616629486848
0.0020388278383890.0037890824633560.0015217635231320.0003548234238250.0067024043708550.000214110851400.0035376083463510.0017007816817940.0020615458590330.004722682745247
0.0001078856
Solution 7
0.20775003020.0299198492−0.0339491324−0.20279503170.21317717070.05684580670.2267650517−0.0977041236−0.03399212000.2532921324
0.04649430500.34898896960.30584184740.40129155130.22840279880.08869702440.20247456580.16872594370.37876526750.1741025236
0.6818033772
0.2559303137991560.375799746775034   0.277228800096808   0.197664668158518   0.443579221165137   0.149296154125337   0.4345092712348090.078260058512172   0.343897367370457   0.428288957805820
0.0018102074713030.0052856465920940.0013811383057120.0030019995514350.0016760299688360.0000583888002630.0023745156573730.0048597976728310.0021331330430300.000928066676778
0.00007961456
Solution 8
−0.03642604440.1232874096−0.03499267860.09592066800.24747761350.08777905340.2453311373−0.1234286095−0.07675431000.0837953112
0.29076047400.25509095340.30655464430.10203621560.19403932320.05827772940.18324283360.19385899900.42162531070.3428082855
0.6614924105
0.2583123777227120.373209790270427   0.278669128269362   0.197734446091387   0.445854845752568   0.149636829778090   0.4297188518939490.070228985270279   0.347004190872835   0.430286378019452
0.0005213187823600.0079190926910690.0000101179821790.0029334058951190.0006125769030440.0005462132554050.0021278927758120.0031725735134940.0010682796677290.003005330807344
0.0000970288Table 7
Comparison of results for neurophysiology application benchmark.
Algorithm
EA
DEIWO
Solutions
Functions values
Error
Solutions
Functions values
Error
Solution 1
−0.82821929960.5446434961−0.00944376590.76336762300.01993259830.1466452805
0.31396360690.12063333410.06523327570.01236817930.04654083230.0330776356
0.1207940789
−0.513317629102986  0.494147884422256   0.857821264907119  −0.869276357234471   0.035896065698353   0.035302984144262
6.476891252397143
E
-
004 1.764830741770584E-004 5.303409416579433E-004 5.954489839492686E-004 3.768837806613808E-004 6.201719809721604E-004
0.0000016131
Solution 2
−0.6512719807−0.6858609598−0.4637572369−0.64508537480.1535909562−0.0036883801
0.36077403230.11345960290.01432913970.04123803430.02046071540.0290928705
0.1462019196
0.474689101876457−0.473073184862286  −0.879846849365024   0.880629205401128   0.021597922714450   0.021393481167399
5.397782221636671
E
-
004 6.939643588305122E-004 1.003432819985314E-004 4.515693669908784E-005 8.794875317431894E-005 6.561341485605571E-005
0.0000007971
Solution 3
0.0425943625−0.1626952821−0.92153247860.9841530788−0.6789794019−0.9070329917
0.14896361100.00497296250.33323206900.00385367110.11836989360.0224932754
0.1477907295
0.5094507817241020.513500001903592   0.860888906816676   0.858125508396048  −0.150516100248646   0.151302300254598
6.698088793095636
E
-
004 6.164011496423427E-005 4.250176942812017E-004 5.847818530452624E-004 3.818324226455561E-004 6.049508942429435E-004
0.0000014868
Solution 4
0.32699111980.02664251620.78868438350.9866658030−0.24032840170.2613854687
0.27105375070.02578076950.13316790030.00839764290.04210235240.2008350695
0.1340463271
−0.8666265987848980.864802076273302  −0.499291890387587   0.502959292343981   0.028130597707313   0.027625633057264
3.340535282903012
E
-
004 8.506808817729805E-004 1.347578366112058E-005 4.419943585350222E-004 3.384546104441034E-005 1.571583020721235E-004
0.0000010566
Solution 5
0.8625703877−0.7176375053−0.22719128010.5169409578−0.13052901290.1532817352
0.20435644830.21776845690.02270517520.14042118570.03520672330.0628718325
0.1151261990
−0.999219226854328−0.999969649944819   0.030375755917241   0.022170611106836   0.036181758812701  −0.036046852107358
6.382501370955085
E
-
004 4.308368076140212E-004 6.212509562937959E-007 5.350559335882965E-005 1.564056883091377E-005 2.982004995065353E-004
0.0000006850
Solution 6
0.76187115760.67753367960.5086028850−0.67138920350.25635430630.0555642759
0.16087544440.09018465030.01691110460.13064839850.06749165580.0585551925
0.0596389890
0.8553934092715790.865828873424174   0.518752115726210   0.500755833107956  −0.003420786702082   0.002831763536266
8.016421956744679
E
-
004 4.160424466184143E-004 1.219573478839314E-004 3.030021758793175E-004 1.726192288883470E-004 2.353934715277509E-004
0.0000010225
Solution 7
0.6609930931−0.48210433120.8042915766−0.8729660781−0.8987020407−0.1909288931
0.08379680940.00550564020.34056283190.23814717380.31412611090.2770687343
0.4711746758
0.7576331596669540.737153307517733  −0.652140850599093  −0.675824989376039   0.015414698314102  −0.015418973671849
7.043063529591498
E
-
004 1.344150494564378E-004 4.842380876615825E-004 5.273519108908182E-004 2.245659170263979E-004 1.077840862562709E-004
0.0000013232
Solution 8
0.27391870540.10775413360.96567343960.9240784300−0.3143660356−0.0314940456
0.00755664880.13446810180.30794290490.00650041710.08319830950.0231155825
0.2152949830
0.643020583613510−0.662772222883848  −0.766232339637761   0.748885036051881  −0.015577794653013  −0.016554865087646
5.874692574154850
E
-
004 9.581664882496455E-005 5.491014584886340E-005 6.779650279244136E-004 2.724638197823467E-004 5.105665891994952E-004
0.0000011549
Solution 9
−0.0838634907−0.14372226500.88472214850.84776454970.1777227339−0.0455327341
0.21023363480.26063918180.09533050600.00003035050.00696295570.0003085023
0.1303553523
−0.943535194739874−0.856135906960048   0.330665259088240   0.516998990062110   0.003385954254385  −0.003000958247795
4.018227192950086
E
-
004 2.566469115463654E-004 2.922768011053473E-004 9.610052402960732E-004 3.374094591558359E-004 1.404462148891126E-004
0.0000014553
Solution 10
0.46123590640.27835846870.9360523694−0.9009125260−0.1421082154−0.2759388163
0.08893259970.11087318320.08521992280.01989553210.11977275080.0090365095
0.0495497431
−0.594047159308489−0.642823012010247  −0.804411318752719  −0.765987085952347   0.004090936206997  −0.004733324210377
3.040278002663932
E
-
005 4.235938430474651E-005 2.090792562180113E-006 3.997031908528915E-004 2.127225008508227E-004 3.369070905144755E-004
0.0000003212
Solution 11
−0.69077417580.8565963646−0.54284005280.5465986672−0.2327716625−0.0607828078
0.22815571500.03252743490.02730824140.03852121040.03182576350.0359158033
0.0583912523
0.818188591099612−0.819888086423616  −0.574896979975124   0.572533090358964   0.007674059020755   0.007253138354942
6.089180991386556
E
-
005 1.061381536349160E-005 9.690837191759410E-005 2.057404583100777E-004 1.258799064935345E-004 1.619105484420356E-004
0.0000001070
Solution 12
−0.8078668904−0.95605627260.5850998782−0.22194390270.0620152964−0.0057942792
0.00500921970.03669730760.01248527080.02763429070.01687848490.0248569233
0.0033499526
0.7285163758984630.710260686224578   0.684731209007270   0.704073008291695   0.022870698246558  −0.022612099542858
4.070614592122102
E
-
004 1.890434011242448E-004 5.496957312890105E-004 7.409173419877374E-004 1.495429352398985E-004 2.800232682831274E-004
0.0000014555Table 8
Comparison of results for chemical equilibrium application benchmark.
Algorithm
EA
DEIWO
Solutions
Functions values
Error
Solutions
Functions
Error
Solution 1
−0.01630874550.26136047090.59815592240.86069838830.0440020125
0.15257724470.37124835490.02655352800.27846940660.1168649339
0.2537179573
0.0470390508171011.869415636872731  −0.265064210870296  −0.8493945482387170.036068394869839
0.0267694033487260.0067471698986298.805250733959802E-004 1.522979575456773E-004 9.201328624012376E-004
0.0007637704
Solution 2
0.33573112850.10159723840.19598077150.52981495840.0069016628
0.34913549530.33884815910.03249191990.28534309850.3380473798
0.4345218367
0.0264747191059153.527094917732416   0.199080066736952  −0.855305091111645   0.036701709037720
0.0097484391996210.0136902448772390.0017269819820170.0050775548385637.862731788029286E-004
0.0003118371
Solution 3
0.32733186760.03965529070.5208586308−0.44427298600.0065409250
0.32068953280.29865241010.07413357150.13311937030.3989667326
0.3799218328
0.0437013691031262.049069927756282   0.252703525295020  −0.852671096585522   0.036246028365058
0.0245104452389430.0085811293139180.0032291564239910.0041953305995700.003804521422204
0.0007169003
Solution 4
0.1626252165−0.3017126041−0.17838890660.71372752510.0278993324
0.02986114150.22407949880.23015760330.09716667580.3805446921
0.3113039788
0.0416082851349302.045381329230565   0.256879114658442  −0.848553325138083   0.036215168144720
0.0180675902570580.0151291750039490.0060781777079390.0085803731376080.005202632121528
0.0006.929642
Solution 5
0.32210512150.2525051386−0.4189363818−0.63316042640.0174557111
0.35107118650.35446811890.01650377310.10355035860.1177149436
0.2740228596
0.0419946729121872.105373855163727  −0.254482330009030   0.850700424653915   0.036482968084062
0.0209602550654710.0098809763622580.0054308209034610.0118752711906180.002615595117843
0.0007143232
Solution 6
0.04291583540.0811055468−0.2662203512−0.77116700690.0303442027
0.04463606060.24783040560.20402119380.02437313770.3396013173
0.2625846864
0.0391426140220702.315919558093960  −0.242200440749196   0.851329829612683   0.036288481660450
0.0209283144093560.0068128941196970.0036386240727270.0019471290577430.001400229470722
0.0005034014
Solution 7
0.7276812579−0.4551167619−0.2113909448−0.48959995650.0075452271
0.37386563860.03040806050.08382213690.17761474140.3755303113
0.3273224789
0.0297771449838893.130317377778047  −0.210038553008973  −0.855206644040595   0.036744343061275
0.0127560302037470.0135741817138280.0011752109805810.0071081316245044.544958713665148E-004
0.0003990879
Solution 8
0.12960883990.2275206857−0.1061140447−0.81249751780.0310264084
0.06601830670.11913519310.23879771640.07924179660.1760779901
0.1698830307
0.0345686113945932.701986140890615   0.225725393195049  −0.853903291996822   0.036519499172778
0.0184140227742850.0058779454601860.0032795977209410.0025569184085240.004983364513138
0.0004157540
Solution 9
0.33672080300.1420207287−0.1427721429−0.84356185340.0349599086
0.27966241090.08564428480.26609090070.02479277410.1028940828
0.2383566164
0.0297002905685783.233972423385453   0.207655335416109  −0.849852570852197   0.036240046204135
0.0170300726214898.550186286035441E-004 0.0060781746882980.0051967115119190.003857570695716
0.0003695853
Solution 10
0.4224008806−0.18890790920.3561384679−0.73900693080.0237712845
0.27129215970.00110630240.18902267870.14141586270.1106294577
0.1772971294
0.0390764195226262.195407554156290  −0.246825556242421   0.857523650896189   0.036753497883680
0.0146045924809380.0233517481262320.0034075415828406.178556688816173E-004 0.005393500809564
0.0007996812
Solution 11
0.23957062530.46375678810.0154833612−0.58613949610.0177082676
0.29754832610.28479940960.14134097420.021220949950.3056084889
0.3034471293
0.0299895773999483.042439401620619   0.212865462882654  −0.854930903089161   0.036718962042137
0.0110741631930890.0165753681319181.046629449772713E-004 0.0070334114253498.967043275154037E-004
0.0004476638
Solution 12
0.16716621050.2358035522−0.1121274369−0.79531860410.0273318802
0.12458895600.02437175630.20793317990.17178176390.1555592488
0.1562964598
0.0262645504971743.679702397787328   0.193832890273980   0.850241727008394   0.036012370926070
0.0148731671602210.0017937620233530.0033735405897900.0054337553768940.008086975223007
0.0003307343
Solution 13
0.4586478321−0.1011456067−0.0115635220−0.758910346030.02454099949
0.33863462040.12044711330.19630754190.17025246240.0117867162
0.2353790429
0.0397637932764412.261109646335777   0.245179433647380   0.851709222649360   0.036466833083980
0.0202735905767620.0087832614367820.0037158089786560.0077905043699550.003080933278725
0.0005721555
Solution 14
0.40648106860.00732467010.08469535600.67268437410.0234023812
0.33925126460.17846498740.18430995450.03108654630.1365650475
0.2344978766
0.0424945752525382.060377057920538  −0.253929806014199   0.846954259610168   0.035785075844889
0.0226941956542790.0076093426038150.0039269589819180.0033194618160770.007645558185311
0.0006578230Table 9
Comparison of results for kinematic application benchmark.
Algorithm
EA
DEIWO
Solutions
Functions values
Error
Solutions
Functions values
Error
Solution 1
−0.06258203370.7777446281−0.05037258280.3805368959−0.5592587603−0.69883388650.39639276750.0861763643
0.39119678240.39257589630.85265427370.54242130970.77421162240.15371057180.91160199770.1519175234
2.8055186567
0.0455073717064930.998085408972999   0.024366076545306  −0.998881362713880  −0.008164213640613  −0.773486615357711   0.001737164451486−1.189751083986727
0.0017545955155670.0032318107089890.0016423175366500.0021693688384920.0042820271043260.0012141332819590.0039822422232590.006623094167720
0.0001004601
Solution 2
−0.15643535250.45071223200.4622139796−0.8818348503−0.65228242840.40828262350.4718261386−0.5070478474
0.77238646430.58321672090.00872553370.20310506970.60569294030.36636824930.35323598020.4646334692
1.8197999485
−0.2099084023445010.978141679315801   0.203396084482184   0.981576540651021  −0.195397812001838   0.225356569439925  −0.346358155642600  −1.054154642146118
0.0008226821895560.0018688880025800.0048624723391100.0016728100915320.0066698560894770.0042745951048250.0033928270533530.006142247670554
0.0001426091
Solution 3
−0.56181778140.84738135170.2226897354−0.0846064846−0.7914861841−0.41110141660.30560983140.4290046184
0.03369437480.23235412640.94325102440.36639136290.00813273330.22066884530.6249670410.4690397335
1.7384318314
−0.2022827508753490.980451072049731   0.204486405102704   0.980955265324429  −0.200955394943197   0.226603999519539  −0.346215357036744  −1.055887922168507
0.002202615985165 0.003098994555294   0.004087922439548   0.002656303324497   0.007030249990859  0.000351392074988  0.002940170368879   0.004663375717434
0.0001181619
Solution 4
0.06083632940.50533987700.3301025811−0.3441350935−0.2611454909−0.13354394410.7518856650−0.3711959678
0.74093054970.63566389460.77260332330.81337406990.01761091300.09294076530.25656636330.260158013
2.3539965598
0.0894226211977120.997785665917926   0.068431332656990  −0.997200694366805   0.035163660058698  −0.842053211908446  −0.019839283580881  −1.151839767619862
0.0035726402931490.000259082400491  0.000907927865151  0.004354292165639  0.003188494681972  0.001119794486023   0.003566152462141   0.004730361789871
0.0000791293
Solution 5
−0.74617426470.53659856980.49070941980.1144124666−0.54332906100.34438416800.6758924483−0.3341809018
0.15528594130.47126624010.74611405270.69170331880.99099174700.02714734960.09003554950.20247922410
2.3132511111
0.0056186233942850.993946653150363   0.118580049821602  −0.993497493618748  −0.114812848523787  −0.699491795013951   0.036947466879746  −1.208913136968972
0.0120384817623460.001991177524501   0.001098498042427   0.000219260012879   0.001036835150341   0.001125049094542  0.000179375216701   0.002011941060148
0.0001565654
Solution 6
−0.3064809352−0.0056167467−0.5007294639−0.0944531990−0.71612653760.63710961000.4792262814−0.4680908930
0.90603788840.74923845610.74034859710.47824137520.88826658770.20810007850.76905628640.4979341279
3.8308029184
0.0777799878539340.997317971174144   0.056964437816277  −0.998437109323217   0.037128587012228  −0.830589842164460  −0.026418733632420  −1.159254624161276
0.0006928621374700.002111917197364   0.000121608449425  0.001744806752775  0.001103067513220   0.003814260073128  0.005523418787933  0.006023917232275
0.0000905605
Solution 7
−0.57392758150.37671420360.1299295442−0.8025240903−0.30267617560.12267795940.7855891471−0.1240877403
0.52869353990.84120472220.33907339790.26434221700.51625064930.26594490700.40299861270.5347985046
1.9576485695
0.0624516941477460.999289247478883  −0.065940504171922  −0.998657897816642   0.031995659386056  −0.777915484740598   0.011189710814883  −1.188813754784150
0.0024792142288370.002927150217360   0.001665746962001  0.001658680908898  0.000472802926307  0.000768554472374   0.004544515610755   0.003312159443232
0.0000526779
Solution 8
0.7822939914−0.49358656570.80296537520.08043856740.92234373730.02969192510.7980078255−0.8324923146
0.14438841300.11161890820.34877624290.14281166700.38029213430.53114144880.03380665820.0083023855
0.6032920001
0.1297363746264300.991078327417070   0.114651095272193  −0.993508195985984   0.094700498150336  −0.914637593030690  −0.052275570550158  −1.115505051433812
0.0009322220229730.004618875277069   0.000203409138437  0.003973280158754   0.001204313453836  0.001429540419579  0.003053228353648   0.004263676032891
0.0000690265
Solution 9
−0.74617426470.60579263810.92714823760.2664085959−0.4794066217−0.33021910100.7692157072−0.1642217317
0.07623924600.22658857490.06942260550.69919575090.02907254660.19216906870.44592738530.7592540224
1.3773369662
0.0284244421965980.997354294976930   0.049494813895825  −0.999159498545526  −0.049400903337364  −0.740573888955816   0.016727719945304  −1.200998766591111
0.0044764613768830.002834673688489   0.000769440136329   0.000760152784294   0.000030652634633   0.001888111028099  0.003940336679685  0.002887323459874
0.0000566727
Solution 10
−0.62053990280.6846519932−0.9933817672−0.0541621937−0.07481527300.37913950333−0.4826393335−0.7830381952
0.14618187710.45555687280.01025912120.99146913160.29204143560.59092707720.83668102150.7918662138
2.5250914998
0.2702763891193420.966742987729601   0.251572811188323  −0.968707517259438   0.233925117931725  −1.122930726803439  −0.122906457137246  −1.014959836413780
0.0076413308397460.002119116346448   0.001683133324140  0.006884785205684  0.002237943072133   0.002719402992536   0.004568280905804   0.000607655156908
0.0001467558Example 7 (combustion application).
We consider the combustion problem for a temperature of 3000°C as proposed in [18]. The problem is described by the following sparse system of equations:
(18)x2+2x6+x9+2x10=10-5,x3+x8=3×10-5,x1+x3+2x5+2x8+x9+x10=5×10-5,x4+2x7=10-5,0.5140437×10-7x5=x12,0.1006932×10-6x6=2x22,0.7816278×10-15x7=x42,0.1496236×10-6x8=x1x3,0.619441×10-7x9=x1x2,0.2089296×10-14x10=x1x22.
The parameters used by the DEIWO approach are presented in Tables 2 and 3. Some of the solutions obtained by the DEIWO approach are presented in Table 10.Table 10
Comparison of results for combustion application benchmark.
Algorithm
EA
DEIWO
Solutions
Functions values
Error
Solutions
Functions values
Error
Solution 1
−0.0552429896−0.00233775330.0455880930−0.12870294720.0539771728−0.01510360790.10631590190.0386267592−0.11449051350.0872294353
0.02741338780.08418485220.14824188930.08391885670.00305178510.00001093170.01656444860.00251842830.00012915150.0000003019
0.0371466690
0.0036729295858850.025839710713681   0.521599141528195   0.004984881426735   0.628058027576599   1.337360448010958  −0.001636439489080  −0.5257287071743031.237228898173834  −1.967867906657906
0.0020436915936200.004159565646108  0.000758296565400   0.001702002448575  0.000013458126816  0.001335246636430  0.000024849042839  0.0019158755803130.000094830798928  0.000002452380746
0.0000304137
Solution 2
−0.03383785580.01856693330.0534924988−0.03927845170.01838822470.0005245892−0.10242696290.0500461848−0.10133611020.0404252678
0.00087946260.10350868370.09556261970.24414237770.00114499950.00068946190.00154279670.00181007890.00062825890.0000116649
0.0794603736
0.000170916648864−0.007197346940591  −0.063594314222215  −0.012650838988592   1.737695584563462  −0.347952826255259   0.006032977400039   0.064865246780411−7.783877291359830   4.243918883408473
0.0008474760060040.001240932558196   0.001348023865310  0.000594884188514   0.000000060112646  0.000103638642450  0.000160043727115  0.0000108596217020.000001712311771   0.000000008853799
0.0000044657
Solution 3
−0.01032093330.00212011080.1207182825−0.02630266790.0044219824−0.08508385790.0053645992−0.04803333240.07322690610.1059498141
0.11706892950.07265495010.20230136960.01558346940.00010652140.00000899830.00069183030.00124591810.00002188600.0000000463
0.0601546074
0.012819319246173−0.019071587101257   0.081305679904288   0.029533158589909   0.443752747560672  −1.815950119305152  −0.016526743634385  −0.082102054652280−5.292468955772288   4.472259351235464
0.0010679209870790.000826374747992  0.002833219569578  0.003530328678861  0.000164312135105  0.000727633722950  0.000872207456297  0.0010422957516250.000244156926303  0.000004662712454
0.0000247769
Solution 4
0.01771987470.00301004240.0676669725−0.04080399030.08525655980.05360566600.1635419218−0.0031889394−0.13907942760.0275601661
0.02625227910.06444803310.13795282660.2862698530.00031398950.00001811530.00166496560.00119905070.00005334610.0000001605
0.1058284540
0.015469554753745−0.027948524788676  −0.031112943145748   0.000991240275462  −3.984894922917517  −4.439885695760647  −0.000857889955385   0.0306339245091366.948027176097301   0.978107982563108
0.0034867750864510.000509018636612   0.001919773451644  0.000734539635309  0.000239511965292  0.001562687142025  0.000000982557284   0.0004813079611010.000432781623864  0.000012083579193
0.0000195602
Solution 5
0.03483577000.10923861080.1250085306−0.0107958218−0.1399310251−0.0177642059−0.07879416050.0917803902−0.0422169082−0.0302349392
0.02898658770.21675892090.00895881650.16839414280.00121353800.02386614990.00011654970.00435475460.00380541370.0004156978
0.0768661970
−0.0012694924827520.013623874571987  −1.175677942247783  −0.030597939250755  −2.528288495278615  −4.876718138697484   0.016294345284874   1.173144422478730−1.968788171623771   5.854744353118305
0.0008781317898570.002563519769053  0.001329398835772   0.001980751318993  0.000001741576241  0.000371710969061  0.000936233886393  0.0014923387797290.000017173451524   0.000000235630459
0.0000162755
Solution 6
0.0172348545−0.0049839785−0.0036835674−0.04016477610.0334303826−0.00495890410.0505724112−0.0076509738−0.11416787240.0544069796
0.02026569990.01136454120.00529921190.06097004620.00029703840.00004968050.00161320920.00006348460.00008589100.0000004281
0.0042879841
−0.0018583951845130.023721236556167   0.653216998921955  −0.018603147142958   3.776633185434082   4.088647518759656   0.008694690681967  −0.657906149925889−5.577580694917199  −1.310802573641980
0.0018204318743190.004719151003933   0.000379406194652  0.001223765779024  0.000003259497212  0.001124982428505  0.000346077083623   0.0012138368869520.000043737933515   0.000001045713511
0.0000300866
Solution 7
−0.14434753550.01371247490.0532523778−0.04075933150.00533406720.03908412610.01965930750.03960940250.0463257617−0.0921334590
0.04607042910.09283178030.04706591550.00145071640.02083621070.00037606000.00166132310.00768685540.00197936480.0000271419
0.0134575843
−0.002441213851907−0.007212781373596  −0.632390997817911   0.021860340733110  −1.135708922430422   3.176958936133236  −0.009642875855651   0.6361762253063249.619896251014470  −7.984037980824944
0.0014846197425420.003755227488414   0.001910664271513   0.002564589021807  0.000006017905472  0.000103728532125  0.000477874496968  0.0015437064767170.000017012045888   0.000000127002218
0.0000291561
Solution 8
0.16120544720.1001108591−0.03035257580.00155415910.04641697090.0906816701−0.0263359667−0.05404778390.0577884947−0.121281367
0.09668995900.08443035980.05204837240.05112777420.02598719380.02004435910.00000241540.00489299240.01613841220.0016156306
0.0231646270
−0.0138806825095570.015121949074709  −0.275566456809759  −0.040530104159630  −1.600723894610612  −2.981776673353555   0.021263311619477   0.275465137078247−0.071979725029653   3.011085635104662
0.0017501475472700.000131319731513  0.000908744309037   0.001986519079324  0.000192755631134  0.000457646932271  0.001642689343190  0.0038250092811740.000209898515312   0.000003174142090
0.0000254721Table 11
Comparison of results for economics modeling application benchmark.
Algorithm
EA
DEIWO
Solutions
Functions values
Error
Solutions
Functions values
Error
Solution 1
−0.1639324−0.38132090.2242448−0.07550940.11710980.0174083−0.0594358−0.22182840.1856304−0.2653962−0.3712114−0.3440810−0.10601680.0218564−0.20287480.0533728−0.05871110.0057098−0.0149290−0.0004102
1.94
E
-
005 9.73E-005 0.00012012.39E-005 5.61E-005 3.89E-005 3.90E-005 9.31E-005 0.00012945.01E-005 0.0001009 0.0001601 6.32E-006 7.96E-006 7.66E-006 2.35E-005 2.21E-005 3.34E-006 6.12E-006 0.4110599
0.4094911811
−6.219110636449961−2.286125269050309  −0.190906907089769   3.293847780216458   0.226175997267840   6.431110110640911   0.873513075442108  −1.765644489553301−6.288552449032001  −1.366427157490679   3.124344370134706  −7.620830976635871   3.649665491921485  −4.964622122726966   8.520732370320570   4.694681906366633−3.567582193815149   1.165936019206677   1.289481831200988   0.000003335851933
0.0002066212882370.000055179798466  0.000179369758855  0.000153340171517   0.000012764866830  0.000366739260435   0.000027744616667   0.0000256556368320.000211906268785  0.000005300945657   0.000160285427757   0.000008959715481   0.000036544006572  0.000213079172405  0.000028339070859   0.0000799611163690.000045923286818  0.000022862241718   0.000004301520459  0.000313249125625
0.0000004629
Solution 2
−0.2071340−0.2251718−0.0910972−0.0028412−0.21103370.4501557−0.02638000.0086212−0.20657000.1663536−0.1450036−0.0743482−0.2007066−0.1451752−0.2078702−0.2750080−0.0422618−0.06021860.0643765−0.0327867
0.00504960.0003158 0.0035983 0.0044357 0.0117217 0.0153466 0.00050130.00010270.00766980.00456370.00096528.36E-006 0.00186780.00163140.00446180.00847730.00145190.00241150.00211060.2107177
0.1866532809
−8.2581782015204603.820226326345229   0.896792437864452  −1.378246977599881  −1.977980138918008   9.317994920873035   7.625886190069895  −9.854325674137856−8.451314609323713  −1.288302214965556  −4.806635040632184   8.714568113773245  −0.240230209985321  −2.242967577834238  −3.026321149977592   7.848935266531751−0.856054140099465   4.362368405116176  −1.207387470433794  −0.000003389324678
0.0001124484212090.000508863338639   0.000443065499734  0.000118951963880   0.000139208990183  0.000243214789175   0.000139948558457   0.0002533038259530.000004024795816  0.000387682069229   0.000182576533214  0.000073658228323   0.000016739957752  0.000163845735074   0.000222130681599  0.0001033772550530.000140635835620  0.000048579832185   0.000004092228149  0.001171744854285
0.0000023157
Solution 3
0.0936580−0.1113572−0.0652960−0.0274100−0.0515078−0.05257120.1674281−0.0284058−0.1587341−0.1284569−0.1326800−0.1138290−0.14305440.0521726−0.2608338−0.1602811−0.1141750−0.1677992−0.11597210.0020995
0.00049091.66e-005 0.00014130.00013871.50e-005 0.00016870.00062897.11e-006 0.00024280.00019810.00028170.00024420.00016960.00013350.00052270.00030380.00024550.00037510.00024540.3975329
0.2694715119
0.410175655605525−2.945267796294842  −3.835549781983572  −2.443847819769871  −5.748058015027992  −1.590111925088804   2.813985490565141  −4.171269624265078−0.630245452070986   8.113507792718970   4.439827399023531   4.827513156671341   7.390298377478004  −8.776390912646111  −3.473157113136982   0.2312514602670285.048928765222839  −2.003672367032949   1.344442609750887  −0.000002463056953
0.0001740672141600.000070546833656  0.000015906476143   0.000141821985185   0.000226460376758  0.000111802212012   0.000123909026492   0.0002706538734010.000001910012918  0.000074219143303  0.000015270345745  0.000060127527268  0.000025057749856   0.000081474905001   0.000034111315945  0.0000075046248430.000000658443143   0.000003576887608  0.000003311438718   0.002359899986076
0.0000057950
Solution 4
−0.26862920.33913400.07325620.0797120−0.11093620.01778940.4220681−0.0583526−0.2610232−0.2838340−0.35798280.0214270−0.62825580.2185146−0.0897853−0.0178795−0.2514783−0.05464660.02750840.0048114
0.00156730.00281172.58e-005 0.00139860.00039910.00141870.00230954.67e-006 0.00176550.00148910.00286330.00141850.00352290.00101390.00082790.00015950.00109440.00029840.00013230.2261284
0.0336841112
0.309460292219735−9.464274256657685   3.672220916230188   3.501423117721803   4.263588177101166   3.919464724708721   5.299860346560947   4.570456826901647−9.419089944307890   4.270196782723423  −9.794235483541998   3.575029940105042  −9.670595944106060   4.067150477104734  −9.508380909500426   2.6383823975987152.824986452927806   2.617827805414225   1.326365492565140  −0.000002314853582
0.0005262413252410.000697399646506   0.000524135889717  0.000291514335086   0.000511322143838   0.000111237908905   0.000002321170209   0.0003570739343130.000119999377139   0.000161166520658  0.000247081697216   0.000079293991687  0.000272037260975  0.000003123855339   0.000049007635690   0.0000379463060220.000020643833165  0.000007010036976  0.000003070341910  0.000162788230767
0.0000017369Example 8 (economics modeling application).
The following modeling problem is considered as difficult and can be scaled up to arbitrary dimensions [18]. The problem is given by the following system of equations:
(19)(xk+∑i=1n-k-1xixi+k)xn-ck=0,1≤k≤n-1,∑l=1n-1xl+1=0.
The  constantck can be randomly chosen. We considered the value 0 for the constants in our experiments and the case of 20 equations as in literature [12].
The parameters used by the DEIWO approach are presented in Tables2 and 3. Some of the solutions obtained by the DEIWO approach are presented in Table 11.
Table2 gives the parameters values used by the DEIWO.


### 5.3. Experimental Results and Discussion

For Example1, we can see from Table 4 that the error of DEIWO is very small. When iteration number is 500, the error of DEIWO is E-020, compared to the traditional numerical method has high precision 16, compared to 14 EA high accuracy, precision is higher than IWO 12.  What is more, all the accuracy for functions values is very high simultaneously. Thus, the solutions acquired can be approximately equal to the theoretical value. Figures 1 and 2 show that the convergence speed of DEIWO is faster than that of IWO, and the optimization precision higher than that of IWO. For Example 2, Table 5 shows the error of DEIWO is E-024,  compared with traditional methods, and high precision of 20, 18 accuracy is higher than the EA, and precision is higher than IWO 16; almost equal to 0 when iteration number is 500. When iteration number reaches 800, the error is 0. Figures 3 and 4 show that IWO easy to fall into local optimum, while DEIWO can effectively jump out of local optimal solution. In addition,  DEIWO can effectively jump out of local optimal, fewer than 200 generations has been largely convergence.  All of this indicates that DEIWO is more effective than traditional numerical method and IWO.Figure 1
Results of 500 iterations.Figure 2
Results of 800 iterations.Figure 3
Results of 500 iterations.Figure 4
Results of 800 iterations.For interval arithmetic benchmark, we can see the errors of all solutions obtained by DEIWO are almost in  E-005. Moreover, the given 8 solutions are stable. Phenomenon of fluctuation does not appear. Thus, DEIWO can effectively resolve high-dimensional equation systems.For neurophysiology application benchmark, 12 solutions are given in Table7. The solutions’ errors of DEIWO are between E-006 and E-007, EA solving accuracy in  E-001 and E-003.  This also indicates that DEIWO is more effective than EA.Table8 shows 14 solutions of DEIWO and EA for chemical equilibrium application. The errors of DEIWO are all E-004, while those of EA are E-001.For kinematic application benchmark, we can see that the errors of EA are very large, such as Solution 6; the error is 3.8, while the errors of DEIWO are between E-004 and E-005. Solutions for combustion application are given in Table10. It is obvious that DEIWO can effectively resolve the nonlinear equation systems, and the error is much smaller than that of EA.The number of equations for economics modeling application is variable. In our experiment, the number is 20. As can be seen from the Table11, even if the equations of variable number is 20, DEIWO by solving the error precision can achieve the fitness of E-006, and even E-007. Than the error of high precision of the solutions of EA for five orders of magnitude. At the same time, precision of functions values obtained by DEIWO is between E-004 and E-006. No functions values are too large or too small.The running time required for our algorithm to converge is presented in Table12. It is measured in seconds using processor: CPU Intel Core i3-370 using frequency: 2.40 GHz, memory: 4 GB, operating system: Windows 7, run software: MATLAB 7.6. The comparison of the number of function evaluations required by the DEIWO for all the examples is shown in Table 13.Table 12
Comparison of results for CPU time required by the DEIWO for all the examples.
Example1
Example2
Example3
Example4
Example5
Example6
Example7
Example8
EA [12]
5.14
5.09
39.07
28.90
32.71
221.09
151.12
640.92
PDS [9]
5.00
5.00
30
30
30
30
30
20
DEIWO
2.27
2.27
15.05
13.11
12.25
16.05
17.25
18.14Table 13
Comparison of the number of function evaluations required by the DEIWO for all the examples.
Example1
Example2
Example3
Example4
Example5
Example6
Example7
Example8
EA [12]
1
1
8
12
14
10
8
4
PDS [9]
1
1
1
2
2
2
2
3
DEIWO
2
2
8
12
14
10
8
4


### 5.4. A Practical Example for Geometric Constraint

The geometric constraint problem is generally divided into three categories: integrity constraint, underconstraint, and overconstraint [19]. Solving geometric constraint problems generally uses the divide and conquer method; the problem is decomposed into a number of small problems, we use numerical methods to solve. This paper uses the example in literature [20]: a circle packing problem. First, transform the geometric constraints into equations, and then the equations are transformed into an optimization problem; finally solve optimization equation. In this way, we do not need to consider if the number of variables and the number of equations are equal; thus, this problem can be effectively solved.Single circle packing problem model is as shown in Figure5, where, the Point  O1 at  the origin of coordinates, the line segment  O1O3 at the x-axis of positive direction, the length of line segment O1O2 is L1, the length of line segment O1O3 is L2, and the length of line segment O2O3 is L3. There is an inscribed circle in triangle O1O2O3, center is O4, and radius is r. The geometric constraint problem can be  transformed and simplified, can obtain the following equations:
(20)f1=x22+y22-L12=0,f2=(L2-x2)2+y22-L32=0,f3=(y2x4-x2r)2(x22+y22)-r2=0,f4=((-y2)(x4-x2)+(L2-x2)(y2-r))2((L2-x2)2+y22)-r2=0.
Solve the equations, and then get the coordinate values of points O1(x1,y1), O2(x2,y2), O3(x3,y3), O4(x4,y4), and r. Parameters are set as G_SIZE=100, P_MAX=150, F=1.06, and CR=0.1. The remaining parameters are the same as in Table 2. Unknown quantity needed to be solved is x2, y2, x4, r, with range is (0,40). Input parameters are L1, L2,  and L3. For the first set of parameters, the fitness function value obtained by DEIWO is 1 precision higher than that of PSO. For the second set of parameters, accuracy is the same. From the third and the fourth set of parameters, we can see that there is conflict. The results of solving geometric constraint are shown in Table 14.Table 14
Results of solving geometric constraint.
L
1
,
L
2
,
L
3
Algorithm
x
2
y
2
x
4
r
F
(
x
)
20,20,20
PSO
10.0000
17.3205
10.0000
5.7735
1.446
E
-
5
DEIWO
10.0000000215
17.3205081063
10.0000000126
5.7735026797
1.4937
E
-
6
10,20,15
PSO
6.8750
7.2618
7.5000
3.2275
9.482
E
-
6
DEIWO
6.8750000954
7.2618437417
7.5000003760
3.2274862446
6.4021
E
-
6
13.6,8.25,16.3
PSO
−0.7677
13.5783
2.7750
2.9363
8.626
E
-
5
DEIWO
0
13.7395290033
2.9816440198
2.9817236110
12.66797
10,20,40
PSO
−20.0294
−0.0169
−6.2739
1.9062
303.5985
DEIWO
0
10.9711605462
4.0800150911
4.0797379306
1.1000
E
+
3Figure 5
Single circle packing problem model.The DEIWO algorithm is very efficient for solving equations systems. The algorithm has the abilities to overcome local optimal solutions and to obtain global optimal solutions.


## 5.1. Testing Platform

The experimental program testing platform included processor: CPU Intel Core i3-370, frequency: 2.40 GHz, memory: 4 GB, operating system: Windows 7, and run software: MATLAB 7.6.


## 5.2. Testing Nonlinear Equation Systems

In order to test the performance of DEIWO for solving nonlinear equation systems, 8 nonlinear equation systems in the literature are used and the testing results are compared with the literature [12]. Table 2 shows the experiment parameters. The search area and the number of iterations are as shown in Table 3. Table 4 to Table 11 is the results of 8 equation systems.Example 1.
Consider the following nonlinear system:(11)f1(x1,x2)=cos(2x1)-cos(2x2)-0.4=0,f2(x1,x2)=2(x2-x1)+sin(2x2)-sin(2x1)-1.2=0.
The results obtained by applying Newton, Secant, Broyden, and Effati methods; evolutionary approach (EA); Probabilistic-driven search; IWO; and the proposed DEIWO method are presented in Table4.Example 2.
Consider the following nonlinear system:(12)f1(x1,x2)=ex1+x1x2-1=0,f2(x1,x2)=sin(x1x2)+x1+x2-1=0.
The results obtained by Effati, evolutionary approach (EA), probabilistic-driven search, IWO, and the proposed DEIWO method are given in Table5.Example 3 (interval arithmetic benchmark).
We consider one benchmark problem proposed from interval arithmetic [13]. The benchmark consists of the following system of equations:(13)0=x1-0.25428722-0.18324757x4x3x9,0=x2-0.37842197-0.16275449x1x10x6,0=x3-0.27162577-0.16955071x1x2x10,0=x4-0.19807914-0.15585316x7x1x6,0=x5-0.44166728-0.19950920x7x6x3,0=x6-0.14654113-0.18922793x8x5x10,0=x7-0.42937161-0.21180486x2x5x8,0=x8-0.07056438-0.17081208x1x7x6,0=x9-0.34504906-0.19612740x10x6x8,0=x10-0.42651102-0.21466544x4x8x1.
Parameters used by the DEIWO approach are listed in Tables2 and 3. The results obtained by evolutionary approach (EA) and the proposed DEIWO method are given in Table 6.Example 4 (neurophysiology application).
We considered the example proposed in [14], which consisted of the following equations:
(14)x12+x32=1,x22+x42=1,x5x33+x6x43=c1,x5x13+x6x23=c2,x5x1x32+x6x42x2=c3,x5x12x3+x6x22x4=c4,
where the constant ci can be randomly chosen. In our experiments, we considered ci=0, i=1,2,3,4, as in the literature [12]. In [15], this problem is used to show the limitations of Newton’s method for which the running time exponentially increases with the size of the initial intervals. We considered the following values for the parameters used by the DEIWO as given in Tables 2 and 3. Some of the solutions obtained by our approach as well as the values of the objective functions for these values are presented in Table 7.Example 5 (chemical equilibrium application).
We consider the chemical equilibrium system as given by the following [16]:
(15)x1x2+x1-3x5=0,2x1x2+x1+x2x32+R8x2-Rx5+2R10x22+R7x2x3+R9x2x4=0,2x2x32+2R5x32-8x5+R6x3+R7x2x3=0,R9x2x4+2x42-4Rx5=0,x1(x2+1)+R10x22+x2x32+R8x2+R5x32+x42-1+R6x3+R7x2x3+R9x2x4=0,
where
(16)R=10,R5=0.193,R6=0.00259740,R7=0.00344840,R8=0.0000179940,R9=0.000215540,R10=0.0000384640.
The parameters used by the DEIWO approach are presented in Tables 2 and 3. Some of the solutions obtained by the DEIWO approach for the chemical equilibrium application are depicted in Table 8.Example 6 (kinematic application).
We consider the kinematic application kin2 as introduced in [17], which describes the inverse position problem for a six-revolute-joint problem in mechanics. The equations describe a denser constraint system and are given as follows:
(17)xi2+xi+12-1=0,a1ix1x3+a2ix1x4+a3ix2x3+a4ix2x4+a5ix2x7+a6ix5x8+a7ix6x7+a8ix6x8+a9ix1+a10ix2+a11ix3+a12ix4+a13ix5+a14ix6+a15ix7+a16ix8+a17i=0,000000000001≤i≤4.
The coefficients, aki, 1≤k≤17, 1≤i≤4, are given in Table 1. The parameters used by the DEIWO approach are presented in Tables 2 and 3. Some of the solutions obtained by the DEIWO approach for the kinematic example kin2 are presented in Table 9.Table 1
Coefficientsaki for the kinematic application.
i
=
1
i
=
2
i
=
3
i
=
4
k
=
1
−0.249150680
+0.125016350
−0.635550077
+1.48947730
k
=
2
+1.609135400
−0.686607360
−0.115719920
+0.23062341
k
=
3
+0.279423430
−0.119228120
−0.666404480
+1.32810730
k
=
4
+1.434801600
−0.719940470
+0.110362110
−0.25864503
k
=
5
+0.000000000
−0.432419270
+0.290702030
+1.16517200
k
=
6
+0.400263840
+0.000000000
+1.258776700
−0.26908494
k
=
7
−0.800527680
+0.000000000
−0.629388360
+0.53816987
k
=
8
+0.000000000
−0.864838550
+0.581404060
+0.58258598
k
=
9
+0.074052388
−0.037157270
+0.195946620
−0.20816985
k
=
10
−0.083050031
+0.035436896
−1.228034200
+2.68683200
k
=
11
−0.386159610
+0.085383482
+0.000000000
−0.69910317
k
=
12
−0.755266030
+0.000000000
−0.079034221
+0.35744413
k
=
13
+0.504201680
−0.039251967
+0.026387877
+1.24991170
k
=
1
4
−1.091628700
+0.000000000
−0.057131430
+1.46773600
k
=
15
+0.000000000
−0.432419270
−1.162808100
+1.16517200
k
=
16
+0.049207290
+0.000000000
+1.258776700
+1.07633970
k
=
17
+0.049207290
+0.013873010
+2.162575000
−0.69686809Table 2
Parameters values used by the DEIWO.
G
_
SIZE
P
_
MAX
σ
final
s
max
s
min
n
F
CR
10
15
0.01
15
1
4
2
0.9Table 3
Benchmarks used in the experiments.
Benchmark
Number of variables
Variables range
Number of iterations
(1) Interval arithmetic
10
[
-
2,2
]
2000
(2) Neurophysiology application
6
[
-
10,10
]
800
(3) Chemical equilibrium application
5
[
-
10,10
]
800
(4) Kinematic application
8
[
-
10,10
]
800
(5) Combustion application
10
[
-
10,10
]
800
(6) Economics modeling application
20
[
-
10,10
]
800Table 4
Comparison of results for Example1.
Method
Solution
Functions values
Error
Newton
(
0.15,0.49
)
(
-
0.00168,0.01497
)
0.000
22156
Secant
(
0.15,0.49
)
(
-
0.00168,0.01497
)
0.000
22156
Broyden
(
0.15,0.49
)
(
-
0.00168,0.01497
)
0.000
22156
Effati
(
0.1575,0.4970
)
(
-
0.005455,0.00739
)
0.0000
84369
EA [12]
(
0.157772,0.49458
)
(
0.001264,0.000969
)
0.00000
25366
PDS [9]
(
0.156520
,
0.493376
)
(
-
0.0000005815
,
-
0.0000008892
)
0.0000003477
IWO
500 iterations
(
0.1564457783
,
0.493379536
6
)
(
0.0000
51022
,
0.000
29976
)
0.0000000
92
56
800 iterations
(
0.1565388795
,
0.493391007
2
)
(
0.0000
12829
,
-
0.0000
2800
8
)
0.0000000000
9
DEIWO
500 iterations
(
0.156520069
7
,
0.493376374
2
)
(
-
0.0000000000
,
-
0.000000000
)
0.0000000000
800 iterations
(
0.156520069
7
,
0.4933763742
)
(
-
0.000000000
,
-
0.000000000
)
0.0000000000Table 5
Comparison of results for Example2.
Method
Solution
Functions values
Error
Effati
(
0.0096,0.9976
)
(
-
0.019223,0.016776
)
0.000
65095
EA [12]
(
-
0.00138,1.0027
)
(
-
0.00276
,
-
6.37
e
-
005
)
0.00000
76216
PDS [9]
(
0.0
,
1.0
)
(
0
,
0
)
0
IWO
500 iterations
(
0.0000353012
,
1.0000971972
)
(
0.0000
70606
,
0.000
1678
)
0.00000000
33
800 iterations
(
0.000053562
8
,
0.999977272
7
)
(
0.000
10712
,
0.000
84396
)
0.0000000
18
6
DEIWO
500 iterations
(
0.0000000000
,
1
.
0000000000
)
(
0.0000000000
,
0.00000000000
)
0.0000000000
800 iterations
(
0,1
)
(
0,0
)
0Table 6
Comparison of results for interval arithmetic benchmark.
Algorithm
EA
DEIWO
Solution
Functions values
Error
Solution
Functions values
Error
Solution 1
0.04649051150.10135683570.0840577820−0.13884603090.4943905739−0.07606851630.2475819110−0.01707481560.00036675350.1481119311
0.20779592400.27697988460.18768632120.33678871140.05303913210.22237305350.18160847520.0878963860.34472003660.2784227489
0.5578712619
0.266267198585396   0.378860657321047   0.279521817014353   0.201184236934515   0.443621074720278   0.147345385930806   0.4335602098777670.076147905678193   0.347936915389233   0.426164139956603
0.0083944895882410.0022825354155090.0006069506301620.000454038870950.0016087850411080.0019199019521440.0014778742457500.0026780167873950.0019500565783850.001222532385146
0.00009717937
Solution 2
0.12248197610.18262006850.2356779803−0.03711504700.37481818560.22133113410.06978130430.0768058043−0.03121538670.1452667120
0.13185527900.19644283610.03649870690.23548901550.06757530640.07399865880.36070382920.00591829790.37674877630.2811693568
0.4739431843
0.2590708561100570.382754529706839   0.274163620083770   0.194046880920087   0.446234352737533   0.147250545938791   0.4309567896938030.070959777303814   0.348016260950456   0.428727381578941
0.0013908674071930.0016706723471030.0046702321651840.0065945286183550.0010960025611250.0018594515638650.0009818536199010.0024128010630460.0020886061068820.001450589369406
0.0000879354
Solution 3
0.06339443990.1017426933−0.1051842285−0.04770599430.41498583260.12151953210.25397771590.0843972823−0.05341329920.0880998746
0.19084366530.27678973670.37690634360.24609007020.02603377510.02560547600.17614864010.13498698510.39863956910.3383563536
0.6396343329
0.2581779734823070.379453934641213   0.282106464391719   0.203187486961515   0.440886347480303   0.151149001651968   0.4380736052048880.075807440998570   0.343828213649516   0.426228284414650
0.0002792359113930.0016751067868830.0034009084520080.0024440221469600.0045076660001450.0019122062055020.0060158253638390.0023230120026440.0021786966880970.001136405771158
0.0000920236
Solution 4
0.19398201990.01521144000.16186543450.00569858090.1904538879−0.16236040330.1864448178−0.04493027060.1675935311−0.0274959004
0.06033352800.36335147260.10974657920.91146537680.25023582290.30894605610.24289922220.11449162850.17741618960.4539962587
1.4462409770
0.2567239775508220.384592485498424   0.276941817209508   0.197489408658780   0.440994124940054   0.147051336865180   0.4311368547949640.072704222159548   0.340803111821165   0.427549623883685
0.0009788984305040.0035435466836580.0018413211436140.0041761183898890.0020837586812790.0008464904602560.0006403136778010.0051424546098790.000247320573489
0.0000663200
Solution 5
0.1169663983−0.0360324410−0.0517944631−0.18259074480.07419020560.250360462900.20430198030.01206070750.187993760800.09312965555
0.13764661610.41489827280.32335368400.38160231220.36694852620.10380356430.22506280070.05952559500.15711045160.3333260708
0.7766780224
0.2559035380827300.382537175786251   0.280411159251412   0.195875663823431   0.443780436204940   0.143389617571427   0.4298459432405990.075841793689244   0.346507927398647   0.427574364325872
0.0018712859668010.0015616894881070.0016886114282790.0046617057634130.0013350145077740.0058746754791760.0022526741607970.0025832411177900.000546906470420.000247274165109
0.0000789255
Solution 6
0.06006249220.06650344530.1163378165−0.04569937750.1649150798−0.12237710450.06665599530.07328665930.07459618230.0578573421
0.19415205260.31184931040.15532713710.24370215880.27656286840.26905055560.36288585610.00280596690.27035110710.3686104960
0.7079949322
0.2557991630284610.377253755510825   0.277018794733188   0.200978421031471   0.438492258899619   0.148958902639442   0.4284636903048680.075053838513460   0.348037273986552   0.422616629486848
0.0020388278383890.0037890824633560.0015217635231320.0003548234238250.0067024043708550.000214110851400.0035376083463510.0017007816817940.0020615458590330.004722682745247
0.0001078856
Solution 7
0.20775003020.0299198492−0.0339491324−0.20279503170.21317717070.05684580670.2267650517−0.0977041236−0.03399212000.2532921324
0.04649430500.34898896960.30584184740.40129155130.22840279880.08869702440.20247456580.16872594370.37876526750.1741025236
0.6818033772
0.2559303137991560.375799746775034   0.277228800096808   0.197664668158518   0.443579221165137   0.149296154125337   0.4345092712348090.078260058512172   0.343897367370457   0.428288957805820
0.0018102074713030.0052856465920940.0013811383057120.0030019995514350.0016760299688360.0000583888002630.0023745156573730.0048597976728310.0021331330430300.000928066676778
0.00007961456
Solution 8
−0.03642604440.1232874096−0.03499267860.09592066800.24747761350.08777905340.2453311373−0.1234286095−0.07675431000.0837953112
0.29076047400.25509095340.30655464430.10203621560.19403932320.05827772940.18324283360.19385899900.42162531070.3428082855
0.6614924105
0.2583123777227120.373209790270427   0.278669128269362   0.197734446091387   0.445854845752568   0.149636829778090   0.4297188518939490.070228985270279   0.347004190872835   0.430286378019452
0.0005213187823600.0079190926910690.0000101179821790.0029334058951190.0006125769030440.0005462132554050.0021278927758120.0031725735134940.0010682796677290.003005330807344
0.0000970288Table 7
Comparison of results for neurophysiology application benchmark.
Algorithm
EA
DEIWO
Solutions
Functions values
Error
Solutions
Functions values
Error
Solution 1
−0.82821929960.5446434961−0.00944376590.76336762300.01993259830.1466452805
0.31396360690.12063333410.06523327570.01236817930.04654083230.0330776356
0.1207940789
−0.513317629102986  0.494147884422256   0.857821264907119  −0.869276357234471   0.035896065698353   0.035302984144262
6.476891252397143
E
-
004 1.764830741770584E-004 5.303409416579433E-004 5.954489839492686E-004 3.768837806613808E-004 6.201719809721604E-004
0.0000016131
Solution 2
−0.6512719807−0.6858609598−0.4637572369−0.64508537480.1535909562−0.0036883801
0.36077403230.11345960290.01432913970.04123803430.02046071540.0290928705
0.1462019196
0.474689101876457−0.473073184862286  −0.879846849365024   0.880629205401128   0.021597922714450   0.021393481167399
5.397782221636671
E
-
004 6.939643588305122E-004 1.003432819985314E-004 4.515693669908784E-005 8.794875317431894E-005 6.561341485605571E-005
0.0000007971
Solution 3
0.0425943625−0.1626952821−0.92153247860.9841530788−0.6789794019−0.9070329917
0.14896361100.00497296250.33323206900.00385367110.11836989360.0224932754
0.1477907295
0.5094507817241020.513500001903592   0.860888906816676   0.858125508396048  −0.150516100248646   0.151302300254598
6.698088793095636
E
-
004 6.164011496423427E-005 4.250176942812017E-004 5.847818530452624E-004 3.818324226455561E-004 6.049508942429435E-004
0.0000014868
Solution 4
0.32699111980.02664251620.78868438350.9866658030−0.24032840170.2613854687
0.27105375070.02578076950.13316790030.00839764290.04210235240.2008350695
0.1340463271
−0.8666265987848980.864802076273302  −0.499291890387587   0.502959292343981   0.028130597707313   0.027625633057264
3.340535282903012
E
-
004 8.506808817729805E-004 1.347578366112058E-005 4.419943585350222E-004 3.384546104441034E-005 1.571583020721235E-004
0.0000010566
Solution 5
0.8625703877−0.7176375053−0.22719128010.5169409578−0.13052901290.1532817352
0.20435644830.21776845690.02270517520.14042118570.03520672330.0628718325
0.1151261990
−0.999219226854328−0.999969649944819   0.030375755917241   0.022170611106836   0.036181758812701  −0.036046852107358
6.382501370955085
E
-
004 4.308368076140212E-004 6.212509562937959E-007 5.350559335882965E-005 1.564056883091377E-005 2.982004995065353E-004
0.0000006850
Solution 6
0.76187115760.67753367960.5086028850−0.67138920350.25635430630.0555642759
0.16087544440.09018465030.01691110460.13064839850.06749165580.0585551925
0.0596389890
0.8553934092715790.865828873424174   0.518752115726210   0.500755833107956  −0.003420786702082   0.002831763536266
8.016421956744679
E
-
004 4.160424466184143E-004 1.219573478839314E-004 3.030021758793175E-004 1.726192288883470E-004 2.353934715277509E-004
0.0000010225
Solution 7
0.6609930931−0.48210433120.8042915766−0.8729660781−0.8987020407−0.1909288931
0.08379680940.00550564020.34056283190.23814717380.31412611090.2770687343
0.4711746758
0.7576331596669540.737153307517733  −0.652140850599093  −0.675824989376039   0.015414698314102  −0.015418973671849
7.043063529591498
E
-
004 1.344150494564378E-004 4.842380876615825E-004 5.273519108908182E-004 2.245659170263979E-004 1.077840862562709E-004
0.0000013232
Solution 8
0.27391870540.10775413360.96567343960.9240784300−0.3143660356−0.0314940456
0.00755664880.13446810180.30794290490.00650041710.08319830950.0231155825
0.2152949830
0.643020583613510−0.662772222883848  −0.766232339637761   0.748885036051881  −0.015577794653013  −0.016554865087646
5.874692574154850
E
-
004 9.581664882496455E-005 5.491014584886340E-005 6.779650279244136E-004 2.724638197823467E-004 5.105665891994952E-004
0.0000011549
Solution 9
−0.0838634907−0.14372226500.88472214850.84776454970.1777227339−0.0455327341
0.21023363480.26063918180.09533050600.00003035050.00696295570.0003085023
0.1303553523
−0.943535194739874−0.856135906960048   0.330665259088240   0.516998990062110   0.003385954254385  −0.003000958247795
4.018227192950086
E
-
004 2.566469115463654E-004 2.922768011053473E-004 9.610052402960732E-004 3.374094591558359E-004 1.404462148891126E-004
0.0000014553
Solution 10
0.46123590640.27835846870.9360523694−0.9009125260−0.1421082154−0.2759388163
0.08893259970.11087318320.08521992280.01989553210.11977275080.0090365095
0.0495497431
−0.594047159308489−0.642823012010247  −0.804411318752719  −0.765987085952347   0.004090936206997  −0.004733324210377
3.040278002663932
E
-
005 4.235938430474651E-005 2.090792562180113E-006 3.997031908528915E-004 2.127225008508227E-004 3.369070905144755E-004
0.0000003212
Solution 11
−0.69077417580.8565963646−0.54284005280.5465986672−0.2327716625−0.0607828078
0.22815571500.03252743490.02730824140.03852121040.03182576350.0359158033
0.0583912523
0.818188591099612−0.819888086423616  −0.574896979975124   0.572533090358964   0.007674059020755   0.007253138354942
6.089180991386556
E
-
005 1.061381536349160E-005 9.690837191759410E-005 2.057404583100777E-004 1.258799064935345E-004 1.619105484420356E-004
0.0000001070
Solution 12
−0.8078668904−0.95605627260.5850998782−0.22194390270.0620152964−0.0057942792
0.00500921970.03669730760.01248527080.02763429070.01687848490.0248569233
0.0033499526
0.7285163758984630.710260686224578   0.684731209007270   0.704073008291695   0.022870698246558  −0.022612099542858
4.070614592122102
E
-
004 1.890434011242448E-004 5.496957312890105E-004 7.409173419877374E-004 1.495429352398985E-004 2.800232682831274E-004
0.0000014555Table 8
Comparison of results for chemical equilibrium application benchmark.
Algorithm
EA
DEIWO
Solutions
Functions values
Error
Solutions
Functions
Error
Solution 1
−0.01630874550.26136047090.59815592240.86069838830.0440020125
0.15257724470.37124835490.02655352800.27846940660.1168649339
0.2537179573
0.0470390508171011.869415636872731  −0.265064210870296  −0.8493945482387170.036068394869839
0.0267694033487260.0067471698986298.805250733959802E-004 1.522979575456773E-004 9.201328624012376E-004
0.0007637704
Solution 2
0.33573112850.10159723840.19598077150.52981495840.0069016628
0.34913549530.33884815910.03249191990.28534309850.3380473798
0.4345218367
0.0264747191059153.527094917732416   0.199080066736952  −0.855305091111645   0.036701709037720
0.0097484391996210.0136902448772390.0017269819820170.0050775548385637.862731788029286E-004
0.0003118371
Solution 3
0.32733186760.03965529070.5208586308−0.44427298600.0065409250
0.32068953280.29865241010.07413357150.13311937030.3989667326
0.3799218328
0.0437013691031262.049069927756282   0.252703525295020  −0.852671096585522   0.036246028365058
0.0245104452389430.0085811293139180.0032291564239910.0041953305995700.003804521422204
0.0007169003
Solution 4
0.1626252165−0.3017126041−0.17838890660.71372752510.0278993324
0.02986114150.22407949880.23015760330.09716667580.3805446921
0.3113039788
0.0416082851349302.045381329230565   0.256879114658442  −0.848553325138083   0.036215168144720
0.0180675902570580.0151291750039490.0060781777079390.0085803731376080.005202632121528
0.0006.929642
Solution 5
0.32210512150.2525051386−0.4189363818−0.63316042640.0174557111
0.35107118650.35446811890.01650377310.10355035860.1177149436
0.2740228596
0.0419946729121872.105373855163727  −0.254482330009030   0.850700424653915   0.036482968084062
0.0209602550654710.0098809763622580.0054308209034610.0118752711906180.002615595117843
0.0007143232
Solution 6
0.04291583540.0811055468−0.2662203512−0.77116700690.0303442027
0.04463606060.24783040560.20402119380.02437313770.3396013173
0.2625846864
0.0391426140220702.315919558093960  −0.242200440749196   0.851329829612683   0.036288481660450
0.0209283144093560.0068128941196970.0036386240727270.0019471290577430.001400229470722
0.0005034014
Solution 7
0.7276812579−0.4551167619−0.2113909448−0.48959995650.0075452271
0.37386563860.03040806050.08382213690.17761474140.3755303113
0.3273224789
0.0297771449838893.130317377778047  −0.210038553008973  −0.855206644040595   0.036744343061275
0.0127560302037470.0135741817138280.0011752109805810.0071081316245044.544958713665148E-004
0.0003990879
Solution 8
0.12960883990.2275206857−0.1061140447−0.81249751780.0310264084
0.06601830670.11913519310.23879771640.07924179660.1760779901
0.1698830307
0.0345686113945932.701986140890615   0.225725393195049  −0.853903291996822   0.036519499172778
0.0184140227742850.0058779454601860.0032795977209410.0025569184085240.004983364513138
0.0004157540
Solution 9
0.33672080300.1420207287−0.1427721429−0.84356185340.0349599086
0.27966241090.08564428480.26609090070.02479277410.1028940828
0.2383566164
0.0297002905685783.233972423385453   0.207655335416109  −0.849852570852197   0.036240046204135
0.0170300726214898.550186286035441E-004 0.0060781746882980.0051967115119190.003857570695716
0.0003695853
Solution 10
0.4224008806−0.18890790920.3561384679−0.73900693080.0237712845
0.27129215970.00110630240.18902267870.14141586270.1106294577
0.1772971294
0.0390764195226262.195407554156290  −0.246825556242421   0.857523650896189   0.036753497883680
0.0146045924809380.0233517481262320.0034075415828406.178556688816173E-004 0.005393500809564
0.0007996812
Solution 11
0.23957062530.46375678810.0154833612−0.58613949610.0177082676
0.29754832610.28479940960.14134097420.021220949950.3056084889
0.3034471293
0.0299895773999483.042439401620619   0.212865462882654  −0.854930903089161   0.036718962042137
0.0110741631930890.0165753681319181.046629449772713E-004 0.0070334114253498.967043275154037E-004
0.0004476638
Solution 12
0.16716621050.2358035522−0.1121274369−0.79531860410.0273318802
0.12458895600.02437175630.20793317990.17178176390.1555592488
0.1562964598
0.0262645504971743.679702397787328   0.193832890273980   0.850241727008394   0.036012370926070
0.0148731671602210.0017937620233530.0033735405897900.0054337553768940.008086975223007
0.0003307343
Solution 13
0.4586478321−0.1011456067−0.0115635220−0.758910346030.02454099949
0.33863462040.12044711330.19630754190.17025246240.0117867162
0.2353790429
0.0397637932764412.261109646335777   0.245179433647380   0.851709222649360   0.036466833083980
0.0202735905767620.0087832614367820.0037158089786560.0077905043699550.003080933278725
0.0005721555
Solution 14
0.40648106860.00732467010.08469535600.67268437410.0234023812
0.33925126460.17846498740.18430995450.03108654630.1365650475
0.2344978766
0.0424945752525382.060377057920538  −0.253929806014199   0.846954259610168   0.035785075844889
0.0226941956542790.0076093426038150.0039269589819180.0033194618160770.007645558185311
0.0006578230Table 9
Comparison of results for kinematic application benchmark.
Algorithm
EA
DEIWO
Solutions
Functions values
Error
Solutions
Functions values
Error
Solution 1
−0.06258203370.7777446281−0.05037258280.3805368959−0.5592587603−0.69883388650.39639276750.0861763643
0.39119678240.39257589630.85265427370.54242130970.77421162240.15371057180.91160199770.1519175234
2.8055186567
0.0455073717064930.998085408972999   0.024366076545306  −0.998881362713880  −0.008164213640613  −0.773486615357711   0.001737164451486−1.189751083986727
0.0017545955155670.0032318107089890.0016423175366500.0021693688384920.0042820271043260.0012141332819590.0039822422232590.006623094167720
0.0001004601
Solution 2
−0.15643535250.45071223200.4622139796−0.8818348503−0.65228242840.40828262350.4718261386−0.5070478474
0.77238646430.58321672090.00872553370.20310506970.60569294030.36636824930.35323598020.4646334692
1.8197999485
−0.2099084023445010.978141679315801   0.203396084482184   0.981576540651021  −0.195397812001838   0.225356569439925  −0.346358155642600  −1.054154642146118
0.0008226821895560.0018688880025800.0048624723391100.0016728100915320.0066698560894770.0042745951048250.0033928270533530.006142247670554
0.0001426091
Solution 3
−0.56181778140.84738135170.2226897354−0.0846064846−0.7914861841−0.41110141660.30560983140.4290046184
0.03369437480.23235412640.94325102440.36639136290.00813273330.22066884530.6249670410.4690397335
1.7384318314
−0.2022827508753490.980451072049731   0.204486405102704   0.980955265324429  −0.200955394943197   0.226603999519539  −0.346215357036744  −1.055887922168507
0.002202615985165 0.003098994555294   0.004087922439548   0.002656303324497   0.007030249990859  0.000351392074988  0.002940170368879   0.004663375717434
0.0001181619
Solution 4
0.06083632940.50533987700.3301025811−0.3441350935−0.2611454909−0.13354394410.7518856650−0.3711959678
0.74093054970.63566389460.77260332330.81337406990.01761091300.09294076530.25656636330.260158013
2.3539965598
0.0894226211977120.997785665917926   0.068431332656990  −0.997200694366805   0.035163660058698  −0.842053211908446  −0.019839283580881  −1.151839767619862
0.0035726402931490.000259082400491  0.000907927865151  0.004354292165639  0.003188494681972  0.001119794486023   0.003566152462141   0.004730361789871
0.0000791293
Solution 5
−0.74617426470.53659856980.49070941980.1144124666−0.54332906100.34438416800.6758924483−0.3341809018
0.15528594130.47126624010.74611405270.69170331880.99099174700.02714734960.09003554950.20247922410
2.3132511111
0.0056186233942850.993946653150363   0.118580049821602  −0.993497493618748  −0.114812848523787  −0.699491795013951   0.036947466879746  −1.208913136968972
0.0120384817623460.001991177524501   0.001098498042427   0.000219260012879   0.001036835150341   0.001125049094542  0.000179375216701   0.002011941060148
0.0001565654
Solution 6
−0.3064809352−0.0056167467−0.5007294639−0.0944531990−0.71612653760.63710961000.4792262814−0.4680908930
0.90603788840.74923845610.74034859710.47824137520.88826658770.20810007850.76905628640.4979341279
3.8308029184
0.0777799878539340.997317971174144   0.056964437816277  −0.998437109323217   0.037128587012228  −0.830589842164460  −0.026418733632420  −1.159254624161276
0.0006928621374700.002111917197364   0.000121608449425  0.001744806752775  0.001103067513220   0.003814260073128  0.005523418787933  0.006023917232275
0.0000905605
Solution 7
−0.57392758150.37671420360.1299295442−0.8025240903−0.30267617560.12267795940.7855891471−0.1240877403
0.52869353990.84120472220.33907339790.26434221700.51625064930.26594490700.40299861270.5347985046
1.9576485695
0.0624516941477460.999289247478883  −0.065940504171922  −0.998657897816642   0.031995659386056  −0.777915484740598   0.011189710814883  −1.188813754784150
0.0024792142288370.002927150217360   0.001665746962001  0.001658680908898  0.000472802926307  0.000768554472374   0.004544515610755   0.003312159443232
0.0000526779
Solution 8
0.7822939914−0.49358656570.80296537520.08043856740.92234373730.02969192510.7980078255−0.8324923146
0.14438841300.11161890820.34877624290.14281166700.38029213430.53114144880.03380665820.0083023855
0.6032920001
0.1297363746264300.991078327417070   0.114651095272193  −0.993508195985984   0.094700498150336  −0.914637593030690  −0.052275570550158  −1.115505051433812
0.0009322220229730.004618875277069   0.000203409138437  0.003973280158754   0.001204313453836  0.001429540419579  0.003053228353648   0.004263676032891
0.0000690265
Solution 9
−0.74617426470.60579263810.92714823760.2664085959−0.4794066217−0.33021910100.7692157072−0.1642217317
0.07623924600.22658857490.06942260550.69919575090.02907254660.19216906870.44592738530.7592540224
1.3773369662
0.0284244421965980.997354294976930   0.049494813895825  −0.999159498545526  −0.049400903337364  −0.740573888955816   0.016727719945304  −1.200998766591111
0.0044764613768830.002834673688489   0.000769440136329   0.000760152784294   0.000030652634633   0.001888111028099  0.003940336679685  0.002887323459874
0.0000566727
Solution 10
−0.62053990280.6846519932−0.9933817672−0.0541621937−0.07481527300.37913950333−0.4826393335−0.7830381952
0.14618187710.45555687280.01025912120.99146913160.29204143560.59092707720.83668102150.7918662138
2.5250914998
0.2702763891193420.966742987729601   0.251572811188323  −0.968707517259438   0.233925117931725  −1.122930726803439  −0.122906457137246  −1.014959836413780
0.0076413308397460.002119116346448   0.001683133324140  0.006884785205684  0.002237943072133   0.002719402992536   0.004568280905804   0.000607655156908
0.0001467558Example 7 (combustion application).
We consider the combustion problem for a temperature of 3000°C as proposed in [18]. The problem is described by the following sparse system of equations:
(18)x2+2x6+x9+2x10=10-5,x3+x8=3×10-5,x1+x3+2x5+2x8+x9+x10=5×10-5,x4+2x7=10-5,0.5140437×10-7x5=x12,0.1006932×10-6x6=2x22,0.7816278×10-15x7=x42,0.1496236×10-6x8=x1x3,0.619441×10-7x9=x1x2,0.2089296×10-14x10=x1x22.
The parameters used by the DEIWO approach are presented in Tables 2 and 3. Some of the solutions obtained by the DEIWO approach are presented in Table 10.Table 10
Comparison of results for combustion application benchmark.
Algorithm
EA
DEIWO
Solutions
Functions values
Error
Solutions
Functions values
Error
Solution 1
−0.0552429896−0.00233775330.0455880930−0.12870294720.0539771728−0.01510360790.10631590190.0386267592−0.11449051350.0872294353
0.02741338780.08418485220.14824188930.08391885670.00305178510.00001093170.01656444860.00251842830.00012915150.0000003019
0.0371466690
0.0036729295858850.025839710713681   0.521599141528195   0.004984881426735   0.628058027576599   1.337360448010958  −0.001636439489080  −0.5257287071743031.237228898173834  −1.967867906657906
0.0020436915936200.004159565646108  0.000758296565400   0.001702002448575  0.000013458126816  0.001335246636430  0.000024849042839  0.0019158755803130.000094830798928  0.000002452380746
0.0000304137
Solution 2
−0.03383785580.01856693330.0534924988−0.03927845170.01838822470.0005245892−0.10242696290.0500461848−0.10133611020.0404252678
0.00087946260.10350868370.09556261970.24414237770.00114499950.00068946190.00154279670.00181007890.00062825890.0000116649
0.0794603736
0.000170916648864−0.007197346940591  −0.063594314222215  −0.012650838988592   1.737695584563462  −0.347952826255259   0.006032977400039   0.064865246780411−7.783877291359830   4.243918883408473
0.0008474760060040.001240932558196   0.001348023865310  0.000594884188514   0.000000060112646  0.000103638642450  0.000160043727115  0.0000108596217020.000001712311771   0.000000008853799
0.0000044657
Solution 3
−0.01032093330.00212011080.1207182825−0.02630266790.0044219824−0.08508385790.0053645992−0.04803333240.07322690610.1059498141
0.11706892950.07265495010.20230136960.01558346940.00010652140.00000899830.00069183030.00124591810.00002188600.0000000463
0.0601546074
0.012819319246173−0.019071587101257   0.081305679904288   0.029533158589909   0.443752747560672  −1.815950119305152  −0.016526743634385  −0.082102054652280−5.292468955772288   4.472259351235464
0.0010679209870790.000826374747992  0.002833219569578  0.003530328678861  0.000164312135105  0.000727633722950  0.000872207456297  0.0010422957516250.000244156926303  0.000004662712454
0.0000247769
Solution 4
0.01771987470.00301004240.0676669725−0.04080399030.08525655980.05360566600.1635419218−0.0031889394−0.13907942760.0275601661
0.02625227910.06444803310.13795282660.2862698530.00031398950.00001811530.00166496560.00119905070.00005334610.0000001605
0.1058284540
0.015469554753745−0.027948524788676  −0.031112943145748   0.000991240275462  −3.984894922917517  −4.439885695760647  −0.000857889955385   0.0306339245091366.948027176097301   0.978107982563108
0.0034867750864510.000509018636612   0.001919773451644  0.000734539635309  0.000239511965292  0.001562687142025  0.000000982557284   0.0004813079611010.000432781623864  0.000012083579193
0.0000195602
Solution 5
0.03483577000.10923861080.1250085306−0.0107958218−0.1399310251−0.0177642059−0.07879416050.0917803902−0.0422169082−0.0302349392
0.02898658770.21675892090.00895881650.16839414280.00121353800.02386614990.00011654970.00435475460.00380541370.0004156978
0.0768661970
−0.0012694924827520.013623874571987  −1.175677942247783  −0.030597939250755  −2.528288495278615  −4.876718138697484   0.016294345284874   1.173144422478730−1.968788171623771   5.854744353118305
0.0008781317898570.002563519769053  0.001329398835772   0.001980751318993  0.000001741576241  0.000371710969061  0.000936233886393  0.0014923387797290.000017173451524   0.000000235630459
0.0000162755
Solution 6
0.0172348545−0.0049839785−0.0036835674−0.04016477610.0334303826−0.00495890410.0505724112−0.0076509738−0.11416787240.0544069796
0.02026569990.01136454120.00529921190.06097004620.00029703840.00004968050.00161320920.00006348460.00008589100.0000004281
0.0042879841
−0.0018583951845130.023721236556167   0.653216998921955  −0.018603147142958   3.776633185434082   4.088647518759656   0.008694690681967  −0.657906149925889−5.577580694917199  −1.310802573641980
0.0018204318743190.004719151003933   0.000379406194652  0.001223765779024  0.000003259497212  0.001124982428505  0.000346077083623   0.0012138368869520.000043737933515   0.000001045713511
0.0000300866
Solution 7
−0.14434753550.01371247490.0532523778−0.04075933150.00533406720.03908412610.01965930750.03960940250.0463257617−0.0921334590
0.04607042910.09283178030.04706591550.00145071640.02083621070.00037606000.00166132310.00768685540.00197936480.0000271419
0.0134575843
−0.002441213851907−0.007212781373596  −0.632390997817911   0.021860340733110  −1.135708922430422   3.176958936133236  −0.009642875855651   0.6361762253063249.619896251014470  −7.984037980824944
0.0014846197425420.003755227488414   0.001910664271513   0.002564589021807  0.000006017905472  0.000103728532125  0.000477874496968  0.0015437064767170.000017012045888   0.000000127002218
0.0000291561
Solution 8
0.16120544720.1001108591−0.03035257580.00155415910.04641697090.0906816701−0.0263359667−0.05404778390.0577884947−0.121281367
0.09668995900.08443035980.05204837240.05112777420.02598719380.02004435910.00000241540.00489299240.01613841220.0016156306
0.0231646270
−0.0138806825095570.015121949074709  −0.275566456809759  −0.040530104159630  −1.600723894610612  −2.981776673353555   0.021263311619477   0.275465137078247−0.071979725029653   3.011085635104662
0.0017501475472700.000131319731513  0.000908744309037   0.001986519079324  0.000192755631134  0.000457646932271  0.001642689343190  0.0038250092811740.000209898515312   0.000003174142090
0.0000254721Table 11
Comparison of results for economics modeling application benchmark.
Algorithm
EA
DEIWO
Solutions
Functions values
Error
Solutions
Functions values
Error
Solution 1
−0.1639324−0.38132090.2242448−0.07550940.11710980.0174083−0.0594358−0.22182840.1856304−0.2653962−0.3712114−0.3440810−0.10601680.0218564−0.20287480.0533728−0.05871110.0057098−0.0149290−0.0004102
1.94
E
-
005 9.73E-005 0.00012012.39E-005 5.61E-005 3.89E-005 3.90E-005 9.31E-005 0.00012945.01E-005 0.0001009 0.0001601 6.32E-006 7.96E-006 7.66E-006 2.35E-005 2.21E-005 3.34E-006 6.12E-006 0.4110599
0.4094911811
−6.219110636449961−2.286125269050309  −0.190906907089769   3.293847780216458   0.226175997267840   6.431110110640911   0.873513075442108  −1.765644489553301−6.288552449032001  −1.366427157490679   3.124344370134706  −7.620830976635871   3.649665491921485  −4.964622122726966   8.520732370320570   4.694681906366633−3.567582193815149   1.165936019206677   1.289481831200988   0.000003335851933
0.0002066212882370.000055179798466  0.000179369758855  0.000153340171517   0.000012764866830  0.000366739260435   0.000027744616667   0.0000256556368320.000211906268785  0.000005300945657   0.000160285427757   0.000008959715481   0.000036544006572  0.000213079172405  0.000028339070859   0.0000799611163690.000045923286818  0.000022862241718   0.000004301520459  0.000313249125625
0.0000004629
Solution 2
−0.2071340−0.2251718−0.0910972−0.0028412−0.21103370.4501557−0.02638000.0086212−0.20657000.1663536−0.1450036−0.0743482−0.2007066−0.1451752−0.2078702−0.2750080−0.0422618−0.06021860.0643765−0.0327867
0.00504960.0003158 0.0035983 0.0044357 0.0117217 0.0153466 0.00050130.00010270.00766980.00456370.00096528.36E-006 0.00186780.00163140.00446180.00847730.00145190.00241150.00211060.2107177
0.1866532809
−8.2581782015204603.820226326345229   0.896792437864452  −1.378246977599881  −1.977980138918008   9.317994920873035   7.625886190069895  −9.854325674137856−8.451314609323713  −1.288302214965556  −4.806635040632184   8.714568113773245  −0.240230209985321  −2.242967577834238  −3.026321149977592   7.848935266531751−0.856054140099465   4.362368405116176  −1.207387470433794  −0.000003389324678
0.0001124484212090.000508863338639   0.000443065499734  0.000118951963880   0.000139208990183  0.000243214789175   0.000139948558457   0.0002533038259530.000004024795816  0.000387682069229   0.000182576533214  0.000073658228323   0.000016739957752  0.000163845735074   0.000222130681599  0.0001033772550530.000140635835620  0.000048579832185   0.000004092228149  0.001171744854285
0.0000023157
Solution 3
0.0936580−0.1113572−0.0652960−0.0274100−0.0515078−0.05257120.1674281−0.0284058−0.1587341−0.1284569−0.1326800−0.1138290−0.14305440.0521726−0.2608338−0.1602811−0.1141750−0.1677992−0.11597210.0020995
0.00049091.66e-005 0.00014130.00013871.50e-005 0.00016870.00062897.11e-006 0.00024280.00019810.00028170.00024420.00016960.00013350.00052270.00030380.00024550.00037510.00024540.3975329
0.2694715119
0.410175655605525−2.945267796294842  −3.835549781983572  −2.443847819769871  −5.748058015027992  −1.590111925088804   2.813985490565141  −4.171269624265078−0.630245452070986   8.113507792718970   4.439827399023531   4.827513156671341   7.390298377478004  −8.776390912646111  −3.473157113136982   0.2312514602670285.048928765222839  −2.003672367032949   1.344442609750887  −0.000002463056953
0.0001740672141600.000070546833656  0.000015906476143   0.000141821985185   0.000226460376758  0.000111802212012   0.000123909026492   0.0002706538734010.000001910012918  0.000074219143303  0.000015270345745  0.000060127527268  0.000025057749856   0.000081474905001   0.000034111315945  0.0000075046248430.000000658443143   0.000003576887608  0.000003311438718   0.002359899986076
0.0000057950
Solution 4
−0.26862920.33913400.07325620.0797120−0.11093620.01778940.4220681−0.0583526−0.2610232−0.2838340−0.35798280.0214270−0.62825580.2185146−0.0897853−0.0178795−0.2514783−0.05464660.02750840.0048114
0.00156730.00281172.58e-005 0.00139860.00039910.00141870.00230954.67e-006 0.00176550.00148910.00286330.00141850.00352290.00101390.00082790.00015950.00109440.00029840.00013230.2261284
0.0336841112
0.309460292219735−9.464274256657685   3.672220916230188   3.501423117721803   4.263588177101166   3.919464724708721   5.299860346560947   4.570456826901647−9.419089944307890   4.270196782723423  −9.794235483541998   3.575029940105042  −9.670595944106060   4.067150477104734  −9.508380909500426   2.6383823975987152.824986452927806   2.617827805414225   1.326365492565140  −0.000002314853582
0.0005262413252410.000697399646506   0.000524135889717  0.000291514335086   0.000511322143838   0.000111237908905   0.000002321170209   0.0003570739343130.000119999377139   0.000161166520658  0.000247081697216   0.000079293991687  0.000272037260975  0.000003123855339   0.000049007635690   0.0000379463060220.000020643833165  0.000007010036976  0.000003070341910  0.000162788230767
0.0000017369Example 8 (economics modeling application).
The following modeling problem is considered as difficult and can be scaled up to arbitrary dimensions [18]. The problem is given by the following system of equations:
(19)(xk+∑i=1n-k-1xixi+k)xn-ck=0,1≤k≤n-1,∑l=1n-1xl+1=0.
The  constantck can be randomly chosen. We considered the value 0 for the constants in our experiments and the case of 20 equations as in literature [12].
The parameters used by the DEIWO approach are presented in Tables2 and 3. Some of the solutions obtained by the DEIWO approach are presented in Table 11.
Table2 gives the parameters values used by the DEIWO.


## 5.3. Experimental Results and Discussion

For Example1, we can see from Table 4 that the error of DEIWO is very small. When iteration number is 500, the error of DEIWO is E-020, compared to the traditional numerical method has high precision 16, compared to 14 EA high accuracy, precision is higher than IWO 12.  What is more, all the accuracy for functions values is very high simultaneously. Thus, the solutions acquired can be approximately equal to the theoretical value. Figures 1 and 2 show that the convergence speed of DEIWO is faster than that of IWO, and the optimization precision higher than that of IWO. For Example 2, Table 5 shows the error of DEIWO is E-024,  compared with traditional methods, and high precision of 20, 18 accuracy is higher than the EA, and precision is higher than IWO 16; almost equal to 0 when iteration number is 500. When iteration number reaches 800, the error is 0. Figures 3 and 4 show that IWO easy to fall into local optimum, while DEIWO can effectively jump out of local optimal solution. In addition,  DEIWO can effectively jump out of local optimal, fewer than 200 generations has been largely convergence.  All of this indicates that DEIWO is more effective than traditional numerical method and IWO.Figure 1
Results of 500 iterations.Figure 2
Results of 800 iterations.Figure 3
Results of 500 iterations.Figure 4
Results of 800 iterations.For interval arithmetic benchmark, we can see the errors of all solutions obtained by DEIWO are almost in  E-005. Moreover, the given 8 solutions are stable. Phenomenon of fluctuation does not appear. Thus, DEIWO can effectively resolve high-dimensional equation systems.For neurophysiology application benchmark, 12 solutions are given in Table7. The solutions’ errors of DEIWO are between E-006 and E-007, EA solving accuracy in  E-001 and E-003.  This also indicates that DEIWO is more effective than EA.Table8 shows 14 solutions of DEIWO and EA for chemical equilibrium application. The errors of DEIWO are all E-004, while those of EA are E-001.For kinematic application benchmark, we can see that the errors of EA are very large, such as Solution 6; the error is 3.8, while the errors of DEIWO are between E-004 and E-005. Solutions for combustion application are given in Table10. It is obvious that DEIWO can effectively resolve the nonlinear equation systems, and the error is much smaller than that of EA.The number of equations for economics modeling application is variable. In our experiment, the number is 20. As can be seen from the Table11, even if the equations of variable number is 20, DEIWO by solving the error precision can achieve the fitness of E-006, and even E-007. Than the error of high precision of the solutions of EA for five orders of magnitude. At the same time, precision of functions values obtained by DEIWO is between E-004 and E-006. No functions values are too large or too small.The running time required for our algorithm to converge is presented in Table12. It is measured in seconds using processor: CPU Intel Core i3-370 using frequency: 2.40 GHz, memory: 4 GB, operating system: Windows 7, run software: MATLAB 7.6. The comparison of the number of function evaluations required by the DEIWO for all the examples is shown in Table 13.Table 12
Comparison of results for CPU time required by the DEIWO for all the examples.
Example1
Example2
Example3
Example4
Example5
Example6
Example7
Example8
EA [12]
5.14
5.09
39.07
28.90
32.71
221.09
151.12
640.92
PDS [9]
5.00
5.00
30
30
30
30
30
20
DEIWO
2.27
2.27
15.05
13.11
12.25
16.05
17.25
18.14Table 13
Comparison of the number of function evaluations required by the DEIWO for all the examples.
Example1
Example2
Example3
Example4
Example5
Example6
Example7
Example8
EA [12]
1
1
8
12
14
10
8
4
PDS [9]
1
1
1
2
2
2
2
3
DEIWO
2
2
8
12
14
10
8
4


## 5.4. A Practical Example for Geometric Constraint

The geometric constraint problem is generally divided into three categories: integrity constraint, underconstraint, and overconstraint [19]. Solving geometric constraint problems generally uses the divide and conquer method; the problem is decomposed into a number of small problems, we use numerical methods to solve. This paper uses the example in literature [20]: a circle packing problem. First, transform the geometric constraints into equations, and then the equations are transformed into an optimization problem; finally solve optimization equation. In this way, we do not need to consider if the number of variables and the number of equations are equal; thus, this problem can be effectively solved.Single circle packing problem model is as shown in Figure5, where, the Point  O1 at  the origin of coordinates, the line segment  O1O3 at the x-axis of positive direction, the length of line segment O1O2 is L1, the length of line segment O1O3 is L2, and the length of line segment O2O3 is L3. There is an inscribed circle in triangle O1O2O3, center is O4, and radius is r. The geometric constraint problem can be  transformed and simplified, can obtain the following equations:
(20)f1=x22+y22-L12=0,f2=(L2-x2)2+y22-L32=0,f3=(y2x4-x2r)2(x22+y22)-r2=0,f4=((-y2)(x4-x2)+(L2-x2)(y2-r))2((L2-x2)2+y22)-r2=0.
Solve the equations, and then get the coordinate values of points O1(x1,y1), O2(x2,y2), O3(x3,y3), O4(x4,y4), and r. Parameters are set as G_SIZE=100, P_MAX=150, F=1.06, and CR=0.1. The remaining parameters are the same as in Table 2. Unknown quantity needed to be solved is x2, y2, x4, r, with range is (0,40). Input parameters are L1, L2,  and L3. For the first set of parameters, the fitness function value obtained by DEIWO is 1 precision higher than that of PSO. For the second set of parameters, accuracy is the same. From the third and the fourth set of parameters, we can see that there is conflict. The results of solving geometric constraint are shown in Table 14.Table 14
Results of solving geometric constraint.
L
1
,
L
2
,
L
3
Algorithm
x
2
y
2
x
4
r
F
(
x
)
20,20,20
PSO
10.0000
17.3205
10.0000
5.7735
1.446
E
-
5
DEIWO
10.0000000215
17.3205081063
10.0000000126
5.7735026797
1.4937
E
-
6
10,20,15
PSO
6.8750
7.2618
7.5000
3.2275
9.482
E
-
6
DEIWO
6.8750000954
7.2618437417
7.5000003760
3.2274862446
6.4021
E
-
6
13.6,8.25,16.3
PSO
−0.7677
13.5783
2.7750
2.9363
8.626
E
-
5
DEIWO
0
13.7395290033
2.9816440198
2.9817236110
12.66797
10,20,40
PSO
−20.0294
−0.0169
−6.2739
1.9062
303.5985
DEIWO
0
10.9711605462
4.0800150911
4.0797379306
1.1000
E
+
3Figure 5
Single circle packing problem model.The DEIWO algorithm is very efficient for solving equations systems. The algorithm has the abilities to overcome local optimal solutions and to obtain global optimal solutions.


## 6. Conclusions

In this paper, based on the IWO algorithm and DE algorithm, we present a new hybrid DEIWO algorithm for solving nonlinear equation systems. The proposed DEIWO approach seems to be very efficient for solving equations systems. We analyzed the case of nonlinear equations systems. We first compared our approach for some simple equations systems having only two equations that were recently used for analyzing the performance of a newly proposed method. The results obtained using the proposed DEIWO optimization approach are very promising. In optimization process, useful information obtained by DE from population is used to guide evolution direction of weeds. Experiments show that the DEIWO is effective and the proposed method could be extended for more higher dimensional systems, although this will also involve increased computational complexity. In a similar manner, we can also solve inequality systems and systems of differential equations, which are part of our future research work.



---

*Source: 757391-2013-12-25.xml*