I
nte
rna
t
io
na
l J
o
urna
l o
f
I
nfo
r
m
a
t
ics a
nd
Co
mm
u
n
ica
t
io
n T
ec
hn
o
lo
g
y
(
I
J
-
I
CT
)
Vo
l.
8
,
No
.
2
,
A
u
g
u
s
t
201
9
,
p
p
.
94
~
1
0
1
I
SS
N:
225
2
-
8776
,
DOI
:
1
0
.
1
1
5
9
1
/i
j
ict.
v
8
i2
.
p
p
9
4
-
101
94
J
o
ur
na
l ho
m
ep
a
g
e
:
h
ttp
:
//ia
e
s
co
r
e.
co
m/jo
u
r
n
a
ls
/in
d
ex
.
p
h
p
/
I
JI
C
T
Ta
ilo
red f
lo
w
er pollina
tion (TF
P)
a
lg
o
rith
m
f
o
r di
m
i
nution o
f
real
pow
er los
s
K
.
L
enin
De
p
a
rtme
n
t
o
f
EE
E
,
P
ra
sa
d
V
.
P
o
t
lu
ri
S
id
d
h
a
rth
a
I
n
stit
u
te o
f
T
e
c
h
n
o
lo
g
y
,
In
d
ia
Art
icle
I
nfo
AB
ST
RAC
T
A
r
ticle
his
to
r
y:
R
ec
eiv
ed
Oct
1
6
,
2
0
1
7
R
ev
i
s
ed
Ma
r
1
8
,
2
0
1
9
A
cc
ep
ted
A
p
r
2
0
,
2
0
1
9
In
th
is
p
a
p
e
r,
T
a
il
o
re
d
F
l
o
w
e
r
P
o
ll
in
a
ti
o
n
(T
F
P
)
a
lg
o
rit
h
m
is
p
ro
p
o
s
e
d
to
so
lv
e
th
e
o
p
ti
m
a
l
re
a
c
ti
v
e
p
o
w
e
r
p
ro
b
lem
.
Co
m
p
risin
g
o
f
th
e
e
le
m
e
n
t
s
o
f
c
h
a
o
s
th
e
o
ry
,
S
h
u
f
f
led
f
ro
g
le
a
p
in
g
se
a
r
c
h
a
n
d
L
e
v
y
F
li
g
h
t,
th
e
p
e
rf
o
rm
a
n
c
e
o
f
th
e
f
lo
we
r
p
o
ll
in
a
ti
o
n
a
lg
o
rit
h
m
h
a
s
b
e
e
n
im
p
ro
v
e
d
.
P
r
o
p
o
se
d
T
F
P
a
lg
o
rit
h
m
h
a
s
b
e
e
n
tes
ted
in
sta
n
d
a
rd
IEE
E
1
1
8
&
p
ra
c
ti
c
a
l
1
9
1
b
u
s
tes
t
sy
ste
m
s
a
n
d
sim
u
latio
n
re
su
lt
s
sh
o
w
c
lea
rl
y
th
e
b
e
tt
e
r
p
e
rf
o
rm
a
n
c
e
o
f
th
e
p
ro
p
o
se
d
a
lg
o
rit
h
m
in
re
d
u
c
i
n
g
th
e
re
a
l
p
o
w
e
r
lo
ss
.
K
ey
w
o
r
d
s
:
C
h
ao
s
t
h
eo
r
y
Flo
w
er
p
o
llin
a
tio
n
al
g
o
r
ith
m
Op
ti
m
al
r
ea
cti
v
e
p
o
w
er
Sh
u
f
f
led
f
r
o
g
leap
in
g
s
ea
r
ch
an
d
lev
y
f
l
ig
h
t
T
r
an
s
m
is
s
io
n
lo
s
s
Co
p
y
rig
h
t
©
2
0
1
9
In
stit
u
te o
f
A
d
v
a
n
c
e
d
E
n
g
i
n
e
e
rin
g
a
n
d
S
c
ien
c
e
.
Al
l
rig
h
ts re
se
rv
e
d
.
C
o
r
r
e
s
p
o
nd
ing
A
uth
o
r
:
K.
L
e
n
in
,
Dep
ar
t
m
en
t o
f
E
E
E
,
P
r
asad
V.
P
o
tlu
r
i Sid
d
h
ar
th
a
I
n
s
tit
u
te
o
f
T
ec
h
n
o
lo
g
y
,
Kan
u
r
u
,
Vij
a
y
a
w
ad
a,
An
d
h
r
a
P
r
ad
esh
-
5
2
0
0
0
7
,
I
n
d
ia.
E
m
ail:
g
k
len
i
n
@
g
m
ai
l.c
o
m
1.
I
NT
RO
D
UCT
I
O
N
Op
ti
m
al
r
ea
cti
v
e
p
o
w
er
d
is
p
at
ch
(
OR
P
D)
p
r
o
b
le
m
is
to
m
i
n
i
m
ize
t
h
e
r
ea
l
p
o
w
er
lo
s
s
a
n
d
b
u
s
v
o
lta
g
e
d
ev
iatio
n
.
V
ar
io
u
s
n
u
m
er
ical
m
et
h
o
d
s
lik
e
t
h
e
g
r
ad
ien
t
m
et
h
o
d
[
1
-
2
]
,
Ne
w
to
n
m
eth
o
d
[
3
]
an
d
lin
ea
r
p
r
o
g
r
am
m
i
n
g
[
4
-
7
]
h
av
e
b
ee
n
ad
o
p
ted
t
o
s
o
lv
e
th
e
o
p
ti
m
al
r
ea
ctiv
e
p
o
w
er
d
is
p
atch
p
r
o
b
le
m
.
B
o
th
th
e
g
r
ad
ien
t
an
d
Ne
w
to
n
m
et
h
o
d
s
h
a
v
e
t
h
e
co
m
p
le
x
it
y
i
n
m
an
a
g
i
n
g
in
eq
u
alit
y
co
n
s
tr
ai
n
ts
.
I
f
li
n
ea
r
p
r
o
g
r
a
m
m
in
g
i
s
ap
p
lied
th
en
th
e
i
n
p
u
t
-
o
u
tp
u
t
f
u
n
ct
io
n
h
a
s
to
b
e
u
tter
ed
as
a
s
et
o
f
l
in
ea
r
f
u
n
ctio
n
s
w
h
ic
h
m
o
s
tl
y
lead
to
lo
s
s
o
f
ac
cu
r
ac
y
.
T
h
e
p
r
o
b
lem
o
f
v
o
ltag
e
s
tab
ilit
y
an
d
co
llap
s
e
p
la
y
a
m
aj
o
r
r
o
le
in
p
o
w
er
s
y
s
te
m
p
la
n
n
in
g
a
n
d
o
p
er
atio
n
[
8
]
.
E
v
o
lu
tio
n
ar
y
al
g
o
r
ith
m
s
s
u
c
h
as
g
e
n
etic
alg
o
r
ith
m
h
av
e
b
ee
n
alr
ea
d
y
p
r
o
p
o
s
ed
to
s
o
lv
e
th
e
r
ea
ctiv
e
p
o
w
er
f
lo
w
p
r
o
b
le
m
[
9
-
1
1
]
.
E
v
o
lu
tio
n
ar
y
alg
o
r
it
h
m
i
s
a
h
eu
r
is
t
ic
ap
p
r
o
ac
h
u
s
ed
f
o
r
m
in
i
m
iza
tio
n
p
r
o
b
lem
s
b
y
u
ti
lizi
n
g
n
o
n
l
in
ea
r
an
d
n
o
n
-
d
if
f
er
en
tiab
le
co
n
tin
u
o
u
s
s
p
ac
e
f
u
n
ctio
n
s
.
I
n
[
1
2
]
,
H
y
b
r
id
d
if
f
er
en
t
ial
ev
o
lu
tio
n
alg
o
r
it
h
m
i
s
p
r
o
p
o
s
e
d
to
im
p
r
o
v
e
th
e
v
o
lta
g
e
s
tab
il
it
y
in
d
e
x
.
I
n
[
1
3
]
B
i
o
g
eo
g
r
ap
h
y
B
ased
al
g
o
r
ith
m
is
p
r
o
j
ec
ted
to
s
o
lv
e
th
e
r
ea
cti
v
e
p
o
w
er
d
is
p
atch
p
r
o
b
lem
.
I
n
[
1
4
]
,
a
f
u
zz
y
b
ased
m
et
h
o
d
is
u
s
ed
to
s
o
lv
e
th
e
o
p
tim
a
l
r
ea
ctiv
e
p
o
w
er
s
ch
ed
u
li
n
g
m
et
h
o
d
.
I
n
[
1
5
]
,
an
i
m
p
r
o
v
ed
ev
o
lu
tio
n
ar
y
p
r
o
g
r
a
m
m
in
g
is
u
s
ed
to
s
o
lv
e
th
e
o
p
ti
m
al
r
ea
cti
v
e
p
o
w
er
d
is
p
atch
p
r
o
b
lem
.
I
n
[
1
6
]
,
th
e
o
p
ti
m
al
r
ea
ctiv
e
p
o
w
er
f
lo
w
p
r
o
b
lem
i
s
s
o
lv
ed
b
y
in
te
g
r
atin
g
a
g
en
e
tic
alg
o
r
ith
m
w
it
h
a
n
o
n
li
n
ea
r
in
ter
io
r
p
o
in
t
m
et
h
o
d
.
I
n
[
1
7
]
,
a
p
atter
n
alg
o
r
ith
m
is
u
s
ed
to
s
o
lv
e
ac
-
dc
o
p
ti
m
a
l
r
ea
ctiv
e
p
o
w
er
f
lo
w
m
o
d
el
w
it
h
t
h
e
g
en
er
ato
r
ca
p
ab
ilit
y
li
m
i
ts
.
I
n
[
1
8
]
,
F
.
C
ap
itan
escu
p
r
o
p
o
s
es
a
t
w
o
-
s
tep
ap
p
r
o
ac
h
to
ev
alu
ate
R
ea
cti
v
e
p
o
w
er
r
eser
v
es
w
i
th
r
esp
ec
t
to
o
p
er
atin
g
co
n
s
tr
ai
n
ts
a
n
d
v
o
ltag
e
s
tab
ilit
y
.
I
n
[
1
9
]
,
a
p
r
o
g
r
a
m
m
i
n
g
b
ased
ap
p
r
o
ac
h
is
u
s
ed
to
s
o
lv
e
th
e
o
p
ti
m
a
l
r
ea
ctiv
e
p
o
w
er
d
is
p
atch
p
r
o
b
lem
.
I
n
[
2
0
]
,
A
.
Kar
g
ar
ia
n
et
al
p
r
ese
n
t
a
p
r
o
b
ab
ilis
tic
alg
o
r
ith
m
f
o
r
o
p
ti
m
al
r
ea
cti
v
e
p
o
w
er
p
r
o
v
is
io
n
i
n
h
y
b
r
id
elec
tr
icit
y
m
ar
k
et
s
w
i
th
u
n
ce
r
tai
n
lo
ad
s
.
T
h
is
p
ap
er
p
r
o
p
o
s
es
T
ailo
r
ed
Flo
w
er
P
o
llin
atio
n
(
T
FP
)
alg
o
r
ith
m
i
s
p
r
o
p
o
s
ed
to
s
o
lv
e
th
e
r
ea
cti
v
e
p
o
w
er
p
r
o
b
le
m
.
T
h
e
b
asic
id
ea
o
f
f
lo
w
er
p
o
lli
n
atio
n
p
r
o
ce
s
s
w
h
ic
h
lead
s
to
th
e
f
o
r
m
u
la
tio
n
o
f
f
lo
w
er
p
o
llin
at
io
n
al
g
o
r
ith
m
(
FP
A
)
[
2
1
]
is
f
ir
s
t
i
n
tr
o
d
u
ce
d
an
d
s
u
b
s
eq
u
en
t
l
y
,
c
h
ao
s
Evaluation Warning : The document was created with Spire.PDF for Python.
I
n
t J
I
n
f
&
C
o
m
m
u
n
T
ec
h
n
o
l
I
SS
N:
2252
-
8776
Ta
ilo
r
ed
flo
w
er p
o
llin
a
tio
n
(
TFP
)
a
lg
o
r
ith
m
fo
r
d
imin
u
tio
n
o
f rea
l p
o
w
er lo
s
s
(
K
.
Len
in
)
95
th
eo
r
y
,
S
h
u
f
f
led
f
r
o
g
leap
i
n
g
s
ea
r
ch
a
n
d
L
e
v
y
Fl
ig
h
t
ar
e
in
tr
o
d
u
ce
d
.
P
r
o
p
o
s
ed
T
F
P
a
lg
o
r
ith
m
h
as
b
ee
n
ev
alu
a
ted
in
s
tan
d
ar
d
I
E
E
E
1
1
8
&
p
r
ac
tical
1
9
1
b
u
s
test
s
y
s
te
m
s
.
Si
m
u
latio
n
r
e
s
u
l
ts
s
h
o
w
t
h
at
o
u
r
p
r
o
p
o
s
ed
ap
p
r
o
ac
h
o
u
tp
er
f
o
r
m
s
all
th
e
e
n
titl
ed
r
ep
o
r
ted
alg
o
r
ith
m
s
in
m
i
n
i
m
izatio
n
o
f
r
ea
l p
o
w
er
lo
s
s
.
2.
P
RO
B
L
E
M
F
O
R
M
UL
AT
I
O
N
T
h
e
o
p
tim
al
p
o
w
er
f
lo
w
p
r
o
b
le
m
is
tr
ea
ted
as
a
g
e
n
er
al
m
i
n
i
m
izatio
n
p
r
o
b
lem
w
it
h
co
n
s
tr
ain
ts
,
a
n
d
ca
n
b
e
m
at
h
e
m
a
ticall
y
w
r
itte
n
in
th
e
f
o
llo
w
i
n
g
f
o
r
m
:
Min
i
m
ize
f
(
x
,
u
)
(
1
)
s
u
b
j
ec
t to
g
(
x
,
u
)
=0
(
2
)
an
d
h
(
x
,
u
)
≤
0
(
3
)
w
h
er
e
f
(
x
,
u
)
i
s
th
e
o
b
j
ec
tiv
e
f
u
n
ctio
n
.
g
(
x
.
u
)
a
n
d
h
(
x
,
u
)
ar
e
r
esp
ec
tiv
el
y
th
e
s
et
o
f
eq
u
a
lit
y
a
n
d
in
eq
u
alit
y
co
n
s
tr
ain
ts
.
x
i
s
th
e
v
ec
to
r
o
f
s
tate
v
ar
iab
les,
an
d
u
i
s
t
h
e
v
ec
to
r
o
f
co
n
tr
o
l v
ar
iab
les.
T
h
e
s
tate
v
ar
iab
les ar
e
th
e
lo
ad
b
u
s
es (
P
Q
b
u
s
es)
v
o
ltag
e
s
,
an
g
le
s
,
th
e
g
en
er
ato
r
r
ea
ctiv
e
p
o
w
er
s
an
d
t
h
e
s
lac
k
ac
tiv
e
g
e
n
er
ato
r
p
o
w
er
:
x
=
(
P
g1
,
θ
2
,
.
.
,
θ
N
,
V
L1
,
.
,
V
LNL
,
Q
g1
,
.
.
,
Q
g
ng
)
T
(
4
)
T
h
e
co
n
tr
o
l
v
ar
iab
les
ar
e
th
e
g
en
er
ato
r
b
u
s
v
o
ltag
e
s
,
th
e
s
h
u
n
t
ca
p
ac
ito
r
s
/r
ea
cto
r
s
an
d
th
e
tr
an
s
f
o
r
m
er
s
tap
-
s
etti
n
g
s
:
u
=
(
V
g
,
T
,
Q
c
)
T
(
5
)
or
u
=
(
V
g1
,
…
,
V
g
ng
,
T
1
,
.
.
,
T
Nt
,
Q
c1
,
.
.
,
Q
cN
c
)
T
(
6
)
w
h
er
e
n
g
,
n
t
an
d
n
c
ar
e
th
e
n
u
m
b
er
o
f
g
en
er
ato
r
s
,
n
u
m
b
e
r
o
f
tap
t
r
an
s
f
o
r
m
er
s
an
d
th
e
n
u
m
b
er
o
f
s
h
u
n
t
co
m
p
e
n
s
ato
r
s
r
esp
ec
ti
v
el
y
.
3.
O
B
J
E
CT
I
V
E
F
UNC
T
I
O
N
3
.
1
.
Act
i
v
e
po
w
er
lo
s
s
T
h
e
o
b
j
ec
tiv
e
o
f
th
e
r
ea
ctiv
e
p
o
w
er
d
is
p
atch
is
to
m
in
i
m
ize
th
e
ac
tiv
e
p
o
w
er
lo
s
s
in
t
h
e
tr
an
s
m
i
s
s
io
n
n
et
w
o
r
k
,
w
h
ich
ca
n
b
e
d
escr
ib
ed
as f
o
llo
w
s
:
=
=
∑
∈
(
2
+
2
−
2
)
(
7
)
or
=
=
∑
−
=
+
∑
−
≠
∈
(
8
)
w
h
er
e
g
k
:
is
th
e
co
n
d
u
c
tan
ce
o
f
b
r
an
ch
b
et
w
ee
n
n
o
d
es
i
an
d
j
,
Nb
r
:
is
th
e
to
tal
n
u
m
b
er
o
f
tr
an
s
m
is
s
io
n
li
n
es
i
n
p
o
w
er
s
y
s
te
m
s
.
P
d
:
is
th
e
to
tal
ac
tiv
e
p
o
w
er
d
e
m
a
n
d
,
P
gi
:
is
t
h
e
g
e
n
er
ato
r
ac
tiv
e
p
o
w
er
o
f
u
n
i
t
i,
an
d
P
g
salc
k
:
i
s
th
e
g
e
n
er
ato
r
ac
tiv
e
p
o
w
er
o
f
s
lack
b
u
s
.
3
.
2
.
Vo
l
t
a
g
e
pro
f
ile
i
m
pro
v
e
m
e
nt
Fo
r
m
in
i
m
iz
in
g
t
h
e
v
o
lta
g
e
d
ev
iatio
n
i
n
P
Q
b
u
s
e
s
,
th
e
o
b
j
ec
t
iv
e
f
u
n
ctio
n
b
ec
o
m
e
s
:
=
+
×
(
9
)
Evaluation Warning : The document was created with Spire.PDF for Python.
I
SS
N
:
2252
-
8776
I
n
t J
I
n
f
&
C
o
m
m
u
n
T
ec
h
n
o
l
,
Vo
l.
8
,
No
.
2,
A
u
g
u
s
t
201
9
:
9
4
–
1
0
1
96
w
h
er
e
ω
v
: is a
w
ei
g
h
t
in
g
f
ac
to
r
o
f
v
o
ltag
e
d
ev
ia
tio
n
.
VD
i
s
t
h
e
v
o
lta
g
e
d
ev
iatio
n
g
i
v
e
n
b
y
:
=
∑
|
−
1
|
=
1
(
1
0
)
3
.
3
.
E
qu
a
lity
Co
ns
t
ra
int
T
h
e
eq
u
alit
y
co
n
s
tr
ai
n
t
g
(
x
,
u
)
o
f
t
h
e
O
R
P
D
p
r
o
b
lem
is
r
ep
r
esen
ted
b
y
t
h
e
p
o
w
er
b
ala
n
ce
eq
u
atio
n
,
w
h
er
e
th
e
to
tal
p
o
w
er
g
en
er
at
i
o
n
m
u
s
t c
o
v
er
t
h
e
to
tal
p
o
w
er
d
em
a
n
d
an
d
t
h
e
p
o
w
er
lo
s
s
e
s
:
=
+
(
1
1
)
T
h
is
eq
u
atio
n
is
s
o
lv
ed
b
y
r
u
n
n
i
n
g
Ne
w
to
n
R
ap
h
s
o
n
lo
ad
f
lo
w
m
et
h
o
d
,
b
y
ca
lcu
lati
n
g
th
e
ac
tiv
e
p
o
w
er
o
f
s
lac
k
b
u
s
to
d
eter
m
i
n
e
ac
tiv
e
p
o
w
er
lo
s
s
.
3
.
4
.
I
nequ
a
lity
Co
ns
t
ra
ints
T
h
e
in
eq
u
alit
y
co
n
s
tr
ai
n
ts
h
(
x
,
u
)
r
ef
lect
th
e
li
m
it
s
o
n
co
m
p
o
n
en
t
s
i
n
t
h
e
p
o
w
er
s
y
s
te
m
a
s
w
ell
as
th
e
li
m
it
s
cr
ea
ted
to
en
s
u
r
e
s
y
s
te
m
s
ec
u
r
it
y
.
Up
p
er
an
d
lo
w
er
b
o
u
n
d
s
o
n
t
h
e
ac
tiv
e
p
o
w
er
o
f
s
lac
k
b
u
s
,
a
n
d
r
ea
ctiv
e
p
o
w
er
o
f
g
e
n
er
ato
r
s
:
≤
≤
(
1
2
)
≤
≤
,
∈
(
1
3
)
Up
p
er
an
d
lo
w
er
b
o
u
n
d
s
o
n
th
e
b
u
s
v
o
lta
g
e
m
a
g
n
it
u
d
es:
≤
≤
,
∈
(
1
4
)
Up
p
er
an
d
lo
w
er
b
o
u
n
d
s
o
n
th
e
tr
an
s
f
o
r
m
er
s
tap
r
atio
s
:
≤
≤
,
∈
(
1
5
)
Up
p
er
an
d
lo
w
er
b
o
u
n
d
s
o
n
th
e
co
m
p
e
n
s
ato
r
s
r
ea
cti
v
e
p
o
w
e
r
s
:
≤
≤
,
∈
(
1
6
)
w
h
er
e
N
is
t
h
e
to
tal
n
u
m
b
er
o
f
b
u
s
e
s
,
N
T
is
t
h
e
to
tal
n
u
m
b
e
r
o
f
T
r
an
s
f
o
r
m
er
s
;
N
c
i
s
th
e
to
tal
n
u
m
b
er
o
f
s
h
u
n
t
r
ea
ctiv
e
co
m
p
en
s
ato
r
s
.
4.
F
L
O
WE
R
P
O
L
L
I
NA
T
I
O
N
AL
G
O
RI
T
H
M
Gen
er
all
y
w
e
u
s
e
th
e
f
o
llo
w
in
g
s
y
s
te
m
s
i
n
Flo
w
er
P
o
llin
atio
n
A
l
g
o
r
ith
m
(
FP
A
)
,
S
y
s
te
m
1
.
B
io
tic
an
d
cr
o
s
s
-
p
o
l
lin
atio
n
h
as
b
ee
n
tr
ea
ted
as
g
lo
b
al
p
o
llin
atio
n
p
r
o
ce
s
s
,
an
d
p
o
llen
-
ca
r
r
y
i
n
g
p
o
llin
ato
r
s
tr
av
el
i
n
a
w
a
y
w
h
i
ch
o
b
ey
s
L
ev
y
f
li
g
h
t
s
.
S
y
s
te
m
2
.
Fo
r
lo
c
al
p
o
llin
atio
n
,
A
-
b
io
tic
an
d
s
el
f
-
p
o
llin
a
tio
n
h
a
s
b
ee
n
u
t
ilized
.
S
y
s
te
m
3
.
P
o
llin
ato
r
s
s
u
ch
as
in
s
ec
t
s
ca
n
d
ev
elo
p
f
lo
w
er
r
eliab
ilit
y
,
w
h
ic
h
is
eq
u
iv
ale
n
t
to
a
r
ep
r
o
d
u
ctio
n
p
r
o
b
a
b
ilit
y
an
d
it is
p
r
o
p
o
r
tio
n
al
to
th
e
s
i
m
ilar
it
y
o
f
t
w
o
f
lo
w
er
s
i
m
p
licated
.
Sy
s
te
m
4
.
T
h
e
co
m
m
u
n
icatio
n
o
f
lo
ca
l
p
o
llin
atio
n
an
d
g
lo
b
al
p
o
llin
atio
n
ca
n
b
e
co
n
tr
o
lled
b
y
a
co
n
tr
o
l
p
r
o
b
a
b
ilit
y
p
∈
[0
,
1
]
,
w
it
h
a
s
lig
h
t b
ias to
w
ar
d
s
lo
ca
l p
o
llin
atio
n
.
S
y
s
te
m
1
an
d
f
lo
w
er
r
eliab
ilit
y
ca
n
b
e
r
ep
r
esen
ted
m
at
h
e
m
a
ticall
y
a
s
+
1
=
+
(
)
(
−
∗
)
(
1
7
)
w
h
er
e
is
th
e
p
o
llen
i
o
r
s
o
lu
ti
o
n
v
ec
to
r
x
i
at
iter
atio
n
t
,
an
d
∗
is
th
e
cu
r
r
en
t
b
est
s
o
lu
tio
n
f
o
u
n
d
a
m
o
n
g
all
s
o
lu
tio
n
s
at
t
h
e
cu
r
r
e
n
t
g
en
er
atio
n
/
iter
atio
n
.
Her
e
γ
is
a
s
c
alin
g
f
ac
to
r
to
co
n
tr
o
l
t
h
e
s
te
p
s
ize.
L
(
λ
)
i
s
t
h
e
p
ar
am
eter
t
h
at
co
r
r
esp
o
n
d
s
to
th
e
s
tr
en
g
t
h
o
f
t
h
e
p
o
llin
atio
n
,
w
h
ic
h
ess
e
n
tial
l
y
i
s
also
th
e
s
te
p
s
ize.
Sin
ce
i
n
s
ec
t
s
m
a
y
m
o
v
e
o
v
er
a
lo
n
g
d
is
ta
n
ce
w
it
h
v
ar
io
u
s
d
is
ta
n
ce
s
tep
s
,
we
ca
n
u
s
e
a
L
ev
y
f
lig
h
t
to
m
i
m
i
c
th
is
c
h
ar
ac
ter
is
tic
e
f
f
i
cien
tl
y
.
W
e
d
r
a
w
L
>
0
f
r
o
m
a
L
e
v
y
d
is
tr
ib
u
t
io
n
~
Г
(
(
/
2
)
)
1
1
+
,
(
≫
0
>
0
)
(
1
8
)
Evaluation Warning : The document was created with Spire.PDF for Python.
I
n
t J
I
n
f
&
C
o
m
m
u
n
T
ec
h
n
o
l
I
SS
N:
2252
-
8776
Ta
ilo
r
ed
flo
w
er p
o
llin
a
tio
n
(
TFP
)
a
lg
o
r
ith
m
fo
r
d
imin
u
tio
n
o
f rea
l p
o
w
er lo
s
s
(
K
.
Len
in
)
97
h
er
e,
Γ
(
λ
)
is
t
h
e
s
ta
n
d
ar
d
g
a
m
m
a
f
u
n
ctio
n
,
a
n
d
th
i
s
d
is
tr
ib
u
t
io
n
is
v
al
id
f
o
r
lar
g
e
s
tep
s
s
>
0
.
T
h
en
,
to
m
o
d
el
th
e
lo
ca
l p
o
llin
atio
n
,
f
o
r
b
o
th
s
y
s
te
m
2
a
n
d
s
y
s
te
m
3
ca
n
b
e
r
ep
r
esen
ted
as
+
1
=
+
∈
(
−
)
(
1
9
)
w
h
er
e
an
d
ar
e
p
o
llen
f
r
o
m
d
i
ff
er
en
t
f
lo
w
er
s
o
f
th
e
s
a
m
e
p
lan
t
s
p
ec
ies.
T
h
is
ess
en
tial
l
y
m
i
m
ics
t
h
e
f
lo
w
er
r
eliab
ilit
y
i
n
a
li
m
ited
n
e
ig
h
b
o
u
r
h
o
o
d
.
Ma
th
e
m
atica
ll
y
,
i
f
an
d
co
m
e
s
f
r
o
m
t
h
e
s
a
m
e
s
p
e
cies
o
r
s
elec
ted
f
r
o
m
th
e
s
a
m
e
p
o
p
u
latio
n
,
th
is
eq
u
iv
ale
n
tl
y
b
ec
o
m
es
a
lo
ca
l
r
an
d
o
m
w
al
k
if
w
e
d
r
a
w
∈
f
r
o
m
a
u
n
i
f
o
r
m
d
is
tr
ib
u
tio
n
i
n
[
0
,
1
]
.
T
h
o
u
g
h
Flo
w
er
p
o
llin
atio
n
p
er
f
o
r
m
a
n
ce
ca
n
o
cc
u
r
at
all
b
alan
ce
,
lo
ca
l
an
d
g
lo
b
al,
n
eig
h
b
o
u
r
i
n
g
f
lo
w
er
p
atch
o
r
f
lo
w
er
s
in
t
h
e
n
o
t
-
so
-
f
ar
-
a
w
a
y
n
ei
g
h
b
o
u
r
h
o
o
d
ar
e
m
o
r
e
lik
el
y
to
b
e
p
o
llin
ated
b
y
lo
ca
l
f
lo
w
er
p
o
llen
t
h
a
n
th
o
s
e
f
ar
a
w
a
y
.
I
n
o
r
d
er
to
m
i
m
i
c
th
is
,
w
e
ca
n
e
ff
ec
ti
v
el
y
u
s
e
a
co
n
tr
o
l p
r
o
b
a
b
ilit
y
(
s
y
s
te
m
4
)
o
r
p
r
o
x
i
m
it
y
p
r
o
b
ab
ilit
y
p
to
s
w
itc
h
b
et
w
ee
n
co
m
m
o
n
g
lo
b
al
p
o
llin
atio
n
to
in
ten
s
iv
e
lo
ca
l
p
o
llin
atio
n
.
T
o
s
tar
t
w
it
h
,
w
e
ca
n
u
s
e
a
r
a
w
v
alu
e
o
f
p
=
0
.
8
as a
n
i
n
itiall
y
v
alu
e.
T
h
e
s
i
m
p
lest
m
et
h
o
d
is
to
u
s
e
a
w
ei
g
h
ted
s
u
m
to
co
m
b
i
n
e
all
m
u
ltip
le
o
b
j
ec
tiv
es
i
n
to
a
co
m
p
o
s
i
te
s
in
g
le
o
b
j
ec
tiv
e
=
∑
=
1
∑
=
1
=
1
,
>
0
(
2
0
)
w
h
er
e
m
is
t
h
e
n
u
m
b
er
o
f
o
b
j
ec
tiv
es a
n
d
w
i
(
i
=
1
,
.
.
.
,
m
)
ar
e
n
o
n
-
n
eg
a
tiv
e
w
e
ig
h
ts
.
FP
A
l
g
o
r
ith
m
f
o
r
s
o
l
v
in
g
o
p
ti
m
al
r
ea
cti
v
e
p
o
w
er
o
p
ti
m
izati
o
n
S
tep
1
.
Ob
jective
min
o
f
(
x
)
,
x
=
(
x
1
,
x
2
, ...,
x
d
)
S
tep
2
.
I
n
itia
liz
e
a
p
o
p
u
la
tio
n
o
f n
flo
w
ers
S
tep
3
.
F
i
n
d
th
e
b
est s
o
lu
tio
n
∗
in
th
e
in
itia
l p
o
p
u
la
tio
n
S
tep
4
.
Defin
e
a
co
n
tr
o
l p
r
o
b
a
b
ilit
y
p
∈
[0
,
1]
S
tep
5
.
Defin
e
a
s
to
p
p
in
g
crit
erio
n
(
a
fixed
n
u
mb
er o
f g
e
n
era
tio
n
s
/
iter
a
tio
n
s
)
S
tep
6
.
w
h
ile
(
t
<
Ma
x
Gen
era
tio
n
)
S
tep
6
.
fo
r
i
=
1
:
n
(
a
ll n
flo
w
e
r
s
in
th
e
p
o
p
u
la
ti
o
n
)
S
tep
7
.
if
r
a
n
d
<
p,
S
tep
8
.
Dra
w
a
(
d
-
d
imen
s
io
n
a
l)
s
tep
ve
cto
r
L
w
h
ich
o
b
ey
s
a
Levy
d
is
tr
ib
u
tio
n
Glo
b
a
l
p
o
llin
a
tio
n
th
r
o
u
g
h
+
1
=
+
(
−
∗
)
else
s
tep
9
.
Dra
w
∈
fr
o
m
a
u
n
ifo
r
m
d
i
s
tr
ib
u
tio
n
in
[
0
,
1
]
s
tep
1
0
.
Do
lo
c
a
l p
o
llin
a
tio
n
th
r
o
u
g
h
+
1
=
+
∈
(
−
)
en
d
if
s
tep
1
0
.
E
va
lu
a
te
n
ew s
o
lu
tio
n
s
s
tep
1
1
.
I
f n
ew s
o
lu
tio
n
s
a
r
e
b
e
tter
,
u
p
d
a
te
th
em
in
t
h
e
p
o
p
u
l
a
tio
n
en
d
fo
r
s
tep
1
2
.
F
in
d
th
e
cu
r
r
en
t b
est s
o
lu
tio
n
∗
en
d
w
h
ile
O
u
tp
u
t
-
b
est s
o
lu
tio
n
h
a
s
b
ee
n
fo
u
n
d
5.
CH
AO
T
I
C
M
AP
S
C
h
ao
s
is
a
r
a
n
d
o
m
s
tate
f
o
u
n
d
i
n
t
h
e
n
o
n
-
lin
ea
r
d
y
n
a
m
ical
d
eter
m
in
is
tic
s
y
s
te
m
,
p
o
s
s
es
s
es
n
o
n
-
p
er
io
d
,
n
o
n
-
co
n
v
er
g
i
n
g
an
d
b
o
u
n
d
ed
p
r
o
p
er
ties
.
T
h
e
u
s
e
o
f
ch
ao
ti
c
s
eq
u
en
ce
s
i
s
m
o
r
e
b
en
ef
ic
i
al
th
an
t
h
e
r
an
d
o
m
s
eq
u
en
ce
s
d
u
e
to
its
n
o
n
-
r
ep
etitio
n
an
d
er
g
o
d
icit
y
p
r
o
p
er
tie
s
.
B
o
r
r
o
w
i
n
g
t
h
e
ad
v
a
n
tag
e
s
o
f
er
g
o
d
icit
y
,
n
o
n
-
r
ep
etitio
n
an
d
r
an
d
o
m
n
es
s
o
f
th
e
ch
ao
tic
s
eq
u
en
ce
s
,
th
e
c
h
ao
tic
m
ap
is
r
ep
lacin
g
t
h
e
r
an
d
o
m
s
eq
u
e
n
ce
s
i
n
g
en
er
ati
n
g
th
e
i
n
itial
p
o
p
u
lati
o
n
in
th
e
FP
A
in
th
is
s
t
u
d
y
.
T
h
is
is
to
en
s
u
r
e
th
at
t
h
e
d
iv
er
s
it
y
o
f
t
h
e
in
itia
l
p
o
p
u
latio
n
ca
n
b
e
im
p
r
o
v
ed
,
w
h
er
e
th
e
d
is
tr
ib
u
tio
n
o
f
th
e
in
itial
p
o
p
u
latio
n
is
m
o
r
e
u
n
i
f
o
r
m
.
T
en
d
if
f
er
en
t
ch
ao
tic
m
ap
s
ar
e
an
d
cir
cle
m
ap
is
s
elec
ted
f
o
r
th
e
i
n
teg
r
ati
o
n
w
ith
FP
A
.
+
1
=
(
+
0
.
2
−
(
0
.
5
2
⁄
)
(
2
)
)
×
(
1
)
(
2
1
)
6.
SH
UF
F
L
E
D
F
RO
G
L
E
AP
I
NG
A
L
G
O
R
I
T
H
M
Sh
u
f
f
led
f
r
o
g
leap
i
n
g
al
g
o
r
ith
m
is
a
b
io
lo
g
ical
ev
o
lu
tio
n
alg
o
r
ith
m
b
ased
o
n
s
w
ar
m
i
n
telli
g
en
ce
.
T
h
e
alg
o
r
ith
m
s
i
m
u
late
s
a
g
r
o
u
p
o
f
f
r
o
g
s
i
n
t
h
e
w
etla
n
d
p
ass
in
g
th
o
u
g
h
t
an
d
f
o
r
ag
i
n
g
b
y
cla
s
s
if
ica
tio
n
o
f
et
h
n
ic
Evaluation Warning : The document was created with Spire.PDF for Python.
I
SS
N
:
2252
-
8776
I
n
t J
I
n
f
&
C
o
m
m
u
n
T
ec
h
n
o
l
,
Vo
l.
8
,
No
.
2,
A
u
g
u
s
t
201
9
:
9
4
–
1
0
1
98
g
r
o
u
p
s
.
I
n
th
e
ex
ec
u
tio
n
o
f
th
e
alg
o
r
ith
m
,
F
f
r
o
g
s
ar
e
g
en
er
ated
at
f
ir
s
t
to
f
o
r
m
a
g
r
o
u
p
,
f
o
r
N
-
d
i
m
e
n
s
io
n
a
l
o
p
tim
izatio
n
p
r
o
b
lem
,
f
r
o
g
i
o
f
t
h
e
g
r
o
u
p
is
r
ep
r
ese
n
ted
as
=
(
1
,
2
,
.
.
,
)
th
e
n
i
n
d
i
v
id
u
al
f
r
o
g
s
i
n
t
h
e
g
r
o
u
p
ar
e
s
o
r
ted
in
d
escen
d
in
g
o
r
d
e
r
ac
co
r
d
in
g
to
f
it
n
ess
v
al
u
es,
to
f
i
n
d
th
e
g
lo
b
al
b
est s
o
lu
tio
n
P
x
.
T
h
e
g
r
o
u
p
is
d
iv
id
ed
in
to
m
eth
n
ic
g
r
o
u
p
s
,
ea
ch
eth
n
ic
g
r
o
u
p
in
c
lu
d
i
n
g
n
f
r
o
g
s
,
s
atis
f
y
in
g
t
h
e
r
elatio
n
F
=
m
×
n
.
T
h
e
r
u
le
o
f
eth
n
ic
g
r
o
u
p
d
i
v
is
io
n
i
s
:
t
h
e
f
ir
s
t
f
r
o
g
i
n
to
th
e
f
ir
s
t
s
u
b
-
g
r
o
u
p
,
th
e
s
ec
o
n
d
f
r
o
g
in
to
th
e
s
ec
o
n
d
s
u
b
-
g
r
o
u
p
,
f
r
o
g
m
i
n
to
s
u
b
-
g
r
o
u
p
m
,
f
r
o
g
m
+
1
in
to
th
e
f
ir
s
t
s
u
b
-
g
r
o
u
p
ag
ain
,
f
r
o
g
m
+
2
in
to
th
e
s
ec
o
n
d
s
u
b
-
g
r
o
u
p
,
an
d
s
o
o
n
,
u
n
til
all
t
h
e
f
r
o
g
s
ar
e
d
iv
id
ed
,
th
e
n
f
in
d
t
h
e
b
es
t
f
r
o
g
i
n
ea
c
h
s
u
b
-
g
r
o
u
p
,
d
en
o
ted
b
y
P
b
; g
et
a
w
o
r
s
t
f
r
o
g
co
r
r
esp
o
n
d
in
g
l
y
,
d
e
n
o
ted
b
y
P
w
.
I
ts
iter
ativ
e
f
o
r
m
u
la
ca
n
b
e
ex
p
r
ess
ed
as:
=
(
)
∗
(
−
)
(
2
2
)
_
−
=
+
,
−
≤
≤
(
2
3
)
w
h
er
e
r
an
d
(
)
r
ep
r
esen
ts
a
r
an
d
o
m
n
u
m
b
er
b
et
w
ee
n
0
an
d
1
,
P
b
r
ep
r
esen
ts
th
e
p
o
s
itio
n
o
f
th
e
b
est
f
r
o
g
,
P
w
r
ep
r
esen
ts
th
e
p
o
s
itio
n
o
f
th
e
w
o
r
s
t
f
r
o
g
,
D
r
ep
r
esen
ts
th
e
d
is
tan
ce
m
o
v
ed
b
y
th
e
w
o
r
s
t
f
r
o
g
,
_
−
is
th
e
i
m
p
r
o
v
ed
p
o
s
itio
n
o
f
t
h
e
f
r
o
g
,
D
m
ax
r
ep
r
ese
n
ts
t
h
e
s
tep
le
n
g
t
h
o
f
f
r
o
g
leap
in
g
.
I
n
th
e
ex
ec
u
t
io
n
o
f
th
e
alg
o
r
it
h
m
,
if
th
e
u
p
d
ated
_
−
is
in
th
e
f
ea
s
i
b
le
s
o
lu
tio
n
s
p
ac
e,
ca
lcu
late
th
e
co
r
r
esp
o
n
d
in
g
f
it
n
e
s
s
v
al
u
e
o
f
_
−
,
if
th
e
co
r
r
esp
o
n
d
in
g
f
it
n
ess
v
al
u
e
o
f
_
−
is
w
o
r
s
e
th
a
n
th
e
co
r
r
esp
o
n
d
in
g
f
it
n
es
s
v
a
lu
e
o
f
P
w
,
th
en
u
s
e
P
w
to
r
ep
lace
P
b
in
eq
u
at
io
n
(
2
2
)
an
d
r
e
-
u
p
d
ate
_
−
;
if
th
er
e
i
s
s
till
n
o
i
m
p
r
o
v
e
m
e
n
t,
th
e
n
r
an
d
o
m
l
y
g
e
n
er
ate
a
n
e
w
f
r
o
g
to
r
ep
lace
P
w
;
r
ep
ea
t
th
e
u
p
d
ate
p
r
o
ce
s
s
u
n
til
s
atis
f
y
in
g
s
to
p
co
n
d
i
tio
n
s
.
7.
L
E
V
Y
F
L
I
G
H
T
L
e
v
y
f
li
g
h
t
is
a
r
a
n
k
o
f
n
o
n
-
G
au
s
s
ian
r
a
n
d
o
m
p
r
o
ce
s
s
es
w
h
o
s
e
ar
b
itra
r
y
w
al
k
s
ar
e
d
r
a
w
n
f
r
o
m
L
e
v
y
s
tab
le
d
is
tr
ib
u
tio
n
.
T
h
is
allo
c
atio
n
is
a
s
i
m
p
le
p
o
w
er
-
la
w
f
o
r
m
u
la
L
(
s
)
~
|
s
|
-
1
-
β
w
h
er
e
0
<
ß
<
2
is
an
in
d
e
x
.
Ma
th
e
m
atica
ll
y
ex
c
la
m
a
tio
n
,
a
ea
s
y
v
er
s
io
n
o
f
L
e
v
y
d
i
s
tr
ib
u
tio
n
ca
n
b
e
d
ef
in
ed
as ,
(
,
,
)
=
{
√
2
0
≤
0
[
−
2
(
−
)
]
1
(
−
)
3
2
⁄
0
<
<
<
∞
(
2
4
)
w
h
er
e
>
0
p
ar
am
eter
i
s
s
ca
le
(
co
n
tr
o
ls
t
h
e
s
ca
le
o
f
d
is
tr
ib
u
tio
n
)
p
ar
am
eter
,
μ
p
ar
a
m
eter
i
s
lo
ca
tio
n
o
r
s
h
if
t
p
ar
am
eter
.
I
n
g
en
er
al,
L
ev
y
d
i
s
tr
ib
u
tio
n
s
h
o
u
ld
b
e
d
ef
in
ed
in
ter
m
s
o
f
Fo
u
r
ier
tr
an
s
f
o
r
m
(
)
=
[
−
|
|
]
,
0
<
≤
2
,
(
2
5
)
w
h
er
e
α
is
a
p
a
r
am
eter
w
it
h
i
n
[
-
1
,
1
]
in
ter
v
al
an
d
k
n
o
w
n
a
s
s
c
ale
f
ac
to
r
.
A
n
i
n
d
ex
o
f
o
s
tab
ilit
y
β
∈
[
0
,
2
]
is
also
r
ef
er
r
ed
to
as
L
ev
y
in
d
e
x
.
I
n
p
ar
ticu
lar
,
f
o
r
β
=
1
,
th
e
in
te
g
r
a
l
ca
n
b
e
ca
r
r
ied
o
u
t
an
al
y
tical
l
y
an
d
is
k
n
o
w
n
as
th
e
C
a
u
ch
y
p
r
o
b
ab
ilit
y
d
is
tr
ib
u
tio
n
.
On
e
m
o
r
e
s
p
ec
ial
ca
s
e
w
h
e
n
β
=
2
,
th
e
d
is
tr
ib
u
tio
n
co
r
r
esp
o
n
d
to
Gau
s
s
ia
n
d
is
tr
ib
u
tio
n
.
β
an
d
α
p
ar
am
ete
r
s
tak
e
a
k
e
y
p
ar
t
in
d
eter
m
in
a
tio
n
o
f
th
e
d
is
tr
ib
u
tio
n
.
T
h
e
p
ar
a
m
eter
β
co
n
tr
o
ls
th
e
s
il
h
o
u
e
tte
o
f
t
h
e
p
r
o
b
ab
ilit
y
d
i
s
tr
ib
u
tio
n
in
s
u
c
h
a
w
a
y
th
at
o
n
e
ca
n
ac
q
u
ir
e
d
if
f
er
en
t
s
h
ap
es
o
f
p
r
o
b
ab
ilit
y
d
is
tr
ib
u
tio
n
,
e
s
p
ec
iall
y
i
n
t
h
e
tail
r
eg
io
n
d
e
p
en
d
i
n
g
o
n
th
e
p
ar
a
m
eter
β
.
T
h
u
s
,
t
h
e
s
m
aller
β
p
ar
a
m
eter
ca
u
s
es
th
e
d
is
tr
ib
u
t
io
n
to
m
a
k
e
lo
n
g
e
r
j
u
m
p
s
s
in
ce
t
h
er
e
w
ill
b
e
lo
n
g
er
tail.
I
t
m
ak
e
s
lo
n
g
er
j
u
m
p
s
f
o
r
s
m
aller
v
a
lu
e
s
w
h
er
ea
s
it m
a
k
es
s
h
o
r
ter
j
u
m
p
s
f
o
r
b
ig
g
er
v
al
u
es.
B
y
L
e
v
y
f
l
ig
h
t,
n
e
w
-
f
an
g
led
s
ta
te
o
f
t
h
e
p
ar
ticle
is
d
esig
n
ed
as,
+
1
=
+
⊕
(
)
(
2
6
)
α
is
th
e
s
tep
s
ize
w
h
ic
h
m
u
s
t
b
e
r
elate
d
to
th
e
s
ca
les
o
f
t
h
e
p
r
o
b
lem
o
f
i
n
ter
est.
I
n
th
e
p
r
o
p
o
s
ed
m
et
h
o
d
α
is
r
an
d
o
m
n
u
m
b
er
f
o
r
all
d
i
m
en
s
io
n
s
o
f
p
ar
ticles.
+
1
=
+
(
(
)
)
⊕
(
)
(
2
7
)
t
h
e
p
r
o
d
u
ct
⊕
m
ea
n
s
en
tr
y
-
w
is
e
m
u
l
tip
licatio
n
s
.
A
n
o
n
-
tr
i
v
ial
s
c
h
e
m
e
o
f
g
e
n
er
atin
g
s
tep
s
ize
s
s
a
m
p
le
s
ar
e
s
u
m
m
ar
ized
as
f
o
llo
w
s
,
+
1
=
+
(
(
)
)
⊕
(
)
~
0
.
01
|
|
1
⁄
(
−
)
(
2
8
)
Evaluation Warning : The document was created with Spire.PDF for Python.
I
n
t J
I
n
f
&
C
o
m
m
u
n
T
ec
h
n
o
l
I
SS
N:
2252
-
8776
Ta
ilo
r
ed
flo
w
er p
o
llin
a
tio
n
(
TFP
)
a
lg
o
r
ith
m
fo
r
d
imin
u
tio
n
o
f rea
l p
o
w
er lo
s
s
(
K
.
Len
in
)
99
w
h
er
e
u
an
d
v
ar
e
d
r
aw
n
f
r
o
m
n
o
r
m
al
d
is
tr
ib
u
tio
n
s
.
T
h
at
is
~
(
0
,
2
)
~
(
0
,
2
)
(
2
9
)
w
it
h
=
{
Г
(
1
+
)
(
/
2
)
Г
[
(
1
+
)
/
2
]
2
(
−
1
)
/
2
}
1
⁄
,
=
1
(
3
0
)
h
er
e
Г
is
s
ta
n
d
ar
d
Ga
m
m
a
f
u
n
ct
io
n
.
On
e
o
f
th
e
i
m
p
o
r
tan
t
p
o
in
t
s
to
b
e
co
n
s
id
er
ed
w
h
ile
p
er
f
o
r
m
i
n
g
d
i
s
tr
ib
u
tio
n
b
y
L
e
v
y
f
li
g
h
ts
i
s
th
e
v
al
u
e
ta
k
en
b
y
t
h
e
β
p
ar
a
m
eter
an
d
it
s
u
b
s
ta
n
tial
l
y
a
f
f
ec
ts
d
is
tr
ib
u
ti
o
n
.
8.
T
AI
L
O
RE
D
F
L
O
W
E
R
P
O
L
L
I
NA
T
I
O
N
(
T
F
P
)
AL
G
O
RIT
H
M
I
n
th
e
T
FP
alg
o
r
ith
m
,
th
e
i
n
itia
l
p
o
p
u
latio
n
is
g
en
er
ated
u
s
i
n
g
th
e
cir
cle
m
ap
,
f
r
o
g
leap
in
g
lo
ca
l
s
ea
r
ch
is
p
er
f
o
r
m
ed
b
y
ea
ch
s
o
l
u
tio
n
an
d
w
h
en
r
a
n
d
>p
,
m
o
d
i
f
ied
L
e
v
y
f
li
g
h
t
w
it
h
in
te
g
r
atio
n
o
f
in
er
tia
w
e
ig
h
t
i
n
g
lo
b
al
p
o
llin
atio
n
i
s
p
er
f
o
r
m
e
d
o
n
th
at
p
ar
ticu
lar
s
o
l
u
tio
n
.
T
h
e
s
tep
s
i
n
v
o
lv
ed
in
t
h
e
T
FP
a
r
e
as f
o
llo
w
s
:
Step
1
.
P
ar
am
eter
I
n
itializatio
n
:
I
n
itia
lize
th
e
r
elev
a
n
t
p
ar
a
m
e
ter
s
o
f
p
o
p
u
latio
n
s
ize,
n
,
d
i
m
en
s
io
n
o
f
s
ea
r
ch
s
p
ac
e,
d
,
m
ax
i
m
u
m
iter
atio
n
,
ma
x_
iter
,
s
w
itc
h
p
r
o
b
ab
ilit
y
,
p
,
r
an
g
e
o
f
s
ea
r
c
h
s
p
ac
e
[
Lb
,
U
b
m
]
,
n
u
m
b
e
r
o
f
m
e
m
ep
lex
e
s
,
m
a
n
d
iter
atio
n
s
w
it
h
i
n
ea
ch
m
e
m
ep
lex
,
it
.
Step
2
.
cr
ea
te
I
n
itial P
o
p
u
latio
n
u
s
in
g
th
e
ch
ao
tic
m
ap
Step
3
.
Fin
d
t
h
e
B
est
So
l
u
tio
n
:
T
h
e
f
it
n
es
s
v
a
lu
e
o
f
ea
ch
s
o
lu
ti
o
n
is
ca
lc
u
lated
a
n
d
th
e
b
es
t
s
o
lu
tio
n
i
s
d
eter
m
in
ed
.
Step
4
.
P
er
f
o
r
m
t
h
e
Fro
g
L
ea
p
in
g
Sea
r
ch
-
Fo
r
ea
ch
s
o
l
u
tio
n
,
s
ea
r
ch
is
p
er
f
o
r
m
ed
.
Step
5
.
P
er
f
o
r
m
t
h
e
Glo
b
al
Sear
ch
o
f
Flo
w
er
P
o
l
lin
atio
n
A
l
g
o
r
ith
m
Step
6
.
Up
d
ate
th
e
So
lu
tio
n
:
T
h
e
f
itn
ess
v
al
u
e
o
f
ea
c
h
n
e
w
s
o
l
u
ti
o
n
is
ev
al
u
ated
.
T
h
e
h
is
to
r
ical
p
o
s
itio
n
i
s
u
p
d
ated
th
r
o
u
g
h
co
m
p
ar
is
o
n
w
it
h
t
h
e
n
e
w
s
o
lu
t
io
n
.
S
u
b
s
eq
u
en
t
l
y
,
t
h
e
b
est s
o
l
u
tio
n
i
s
u
p
d
ated
.
Step
7
.
C
h
ec
k
T
er
m
i
n
atio
n
C
o
n
d
it
io
n
o
r
else,
Step
4
is
r
ep
ea
ted
.
9.
SI
M
UL
AT
I
O
N
R
E
S
UL
T
S
A
t
f
ir
s
t
T
ailo
r
ed
Flo
w
er
P
o
llin
atio
n
(
T
FP
)
alg
o
r
ith
m
h
as
b
ee
n
test
ed
i
n
s
ta
n
d
ar
d
I
E
E
E
1
1
8
-
b
u
s
test
s
y
s
te
m
[
22
]
.
T
h
e
s
y
s
te
m
h
a
s
5
4
g
en
er
ato
r
b
u
s
es,
6
4
lo
ad
b
u
s
es,
1
8
6
b
r
an
ch
es
an
d
9
o
f
th
e
m
ar
e
w
i
th
t
h
e
tap
s
etti
n
g
t
r
an
s
f
o
r
m
er
s
.
T
h
e
li
m
i
ts
o
f
v
o
lta
g
e
o
n
g
e
n
er
ato
r
b
u
s
es
ar
e
0
.
9
5
-
1
.
1
p
er
-
u
n
it.,
a
n
d
o
n
lo
ad
b
u
s
es
ar
e
0
.
9
5
-
1
.
0
5
p
er
-
u
n
it.
T
h
e
li
m
it
o
f
tr
an
s
f
o
r
m
er
r
ate
is
0
.
9
-
1
.
1
,
w
it
h
th
e
c
h
an
g
es st
ep
o
f
0
.
0
2
5
.
T
h
e
li
m
itatio
n
s
o
f
r
ea
ctiv
e
p
o
w
er
s
o
u
r
ce
ar
e
lis
te
d
in
T
ab
l
e
1
,
w
i
th
t
h
e
ch
a
n
g
e
i
n
s
tep
o
f
0
.
0
1
.
T
ab
le
1
.
L
i
m
itatio
n
o
f
R
ea
c
tiv
e
P
o
w
er
So
u
r
ce
s
B
U
S
5
34
37
44
45
46
48
Q
C
M
A
X
0
14
0
10
10
10
15
Q
C
M
I
N
-
40
0
-
25
0
0
0
0
B
U
S
74
79
82
83
1
0
5
1
0
7
1
1
0
Q
C
M
A
X
12
20
20
10
20
6
6
Q
C
M
I
N
0
0
0
0
0
0
0
T
h
e
s
tatis
tical
co
m
p
ar
is
o
n
r
e
s
u
lts
o
f
5
0
tr
ial
r
u
n
s
h
a
v
e
b
ee
n
l
is
t
i
n
T
ab
le
2
an
d
th
e
r
esu
lts
cl
ea
r
l
y
s
h
o
w
th
e
b
etter
p
er
f
o
r
m
an
ce
o
f
p
r
o
p
o
s
ed
T
F
P
alg
o
r
ith
m
.
T
ab
le
2
.
C
o
m
p
ar
is
o
n
R
es
u
lt
s
A
c
t
i
v
e
p
o
w
e
r
l
o
ss (p
.
u
)
BBO
[
2
3
]
I
L
S
B
B
O
/
st
r
a
t
e
g
y
1
[
2
3
]
I
L
S
B
B
O
/
st
r
a
t
e
g
y
1
[
2
3
]
P
r
o
p
o
se
d
TF
P
M
i
n
1
2
8
.
7
7
1
2
6
.
9
8
1
2
4
.
7
8
1
1
6
.
0
2
M
a
x
1
3
2
.
6
4
1
3
7
.
3
4
1
3
2
.
3
9
1
1
9
.
1
2
A
v
e
r
a
g
e
1
3
0
.
2
1
1
3
0
.
3
7
1
2
9
.
2
2
1
1
7
.
3
2
T
h
en
th
e
T
ailo
r
ed
Flo
w
er
P
o
llin
atio
n
(
T
FP
)
alg
o
r
ith
m
h
as
b
ee
n
te
s
ted
in
p
r
ac
tical
1
9
1
test
s
y
s
te
m
a
n
d
th
e
f
o
llo
w
in
g
r
es
u
lt
s
h
av
e
b
e
en
o
b
tain
ed
.
I
n
P
r
ac
tical
1
9
1
test
b
u
s
s
y
s
te
m
–
Nu
m
b
er
o
f
Ge
n
er
ato
r
s
=
2
0
,
Nu
m
b
er
o
f
li
n
es
=
2
0
0
,
Nu
m
b
er
o
f
b
u
s
es
=
1
9
1
Nu
m
b
er
o
f
t
r
an
s
m
is
s
io
n
l
in
e
s
=
5
5
.
T
a
b
le
3
s
h
o
w
s
t
h
e
o
p
ti
m
a
l
Evaluation Warning : The document was created with Spire.PDF for Python.
I
SS
N
:
2252
-
8776
I
n
t J
I
n
f
&
C
o
m
m
u
n
T
ec
h
n
o
l
,
Vo
l.
8
,
No
.
2,
A
u
g
u
s
t
201
9
:
9
4
–
1
0
1
100
co
n
tr
o
l
v
al
u
es
o
f
p
r
ac
tical
1
9
1
test
s
y
s
te
m
o
b
tain
ed
b
y
T
FP
m
e
th
o
d
.
An
d
tab
le
4
s
h
o
w
s
th
e
r
es
u
lt
s
ab
o
u
t
th
e
v
alu
e
o
f
th
e
r
ea
l p
o
w
er
lo
s
s
b
y
o
b
tain
ed
b
y
T
ailo
r
ed
Flo
w
er
P
o
llin
atio
n
(
T
FP
)
alg
o
r
ith
m
.
T
ab
le
3
.
Op
tim
al
C
o
n
tr
o
l V
al
u
es
o
f
P
r
ac
tical
1
9
1
Utilit
y
(
I
n
d
ian
)
S
y
s
te
m
b
y
T
FP
Me
th
o
d
V
G
1
1
.
1
0
V
G
1
1
0
.
9
0
V
G
2
0
.
7
8
V
G
1
2
1
.
0
0
V
G
3
1
.
0
1
V
G
1
3
1
.
0
0
V
G
4
1
.
0
1
V
G
1
4
0
.
9
0
V
G
5
1
.
1
0
V
G
1
5
1
.
0
0
V
G
6
1
.
1
0
V
G
1
6
1
.
0
0
V
G
7
1
.
1
0
V
G
1
7
0
.
9
0
V
G
8
1
.
0
1
V
G
1
8
1
.
0
0
V
G
9
1
.
1
0
V
G
1
9
1
.
1
0
V
G
1
0
1
.
0
1
V
G
2
0
1
.
1
0
T1
1
.
0
0
T
2
1
0
.
9
0
T
4
1
0
.
9
0
T2
1
.
0
0
T
2
2
0
.
9
0
T
4
2
0
.
9
0
T3
1
.
0
0
T
2
3
0
.
9
0
T
4
3
0
.
9
1
T4
1
.
1
0
T
2
4
0
.
9
0
T
4
4
0
.
9
1
T5
1
.
0
0
T
2
5
0
.
9
0
T
4
5
0
.
9
1
T6
1
.
0
0
T
2
6
1
.
0
0
T
4
6
0
.
9
0
T7
1
.
0
0
T
2
7
0
.
9
0
T
4
7
0
.
9
1
T8
1
.
0
1
T
2
8
0
.
9
0
T
4
8
1
.
0
0
T9
1
.
0
0
T
2
9
1
.
0
1
T
4
9
0
.
9
0
T
1
0
1
.
0
0
T
3
0
0
.
9
0
T
5
0
0
.
9
0
T
1
1
0
.
9
0
T
3
1
0
.
9
0
T
5
1
0
.
9
0
T
1
2
1
.
0
0
T
3
2
0
.
9
0
T
5
2
0
.
9
0
T
1
3
1
.
0
1
T
3
3
1
.
0
1
T
5
3
1
.
0
0
T
1
4
1
.
0
1
T
3
4
0
.
9
0
T
5
4
0
.
9
0
T
1
5
1
.
0
1
T
3
5
0
.
9
0
T
5
5
0
.
9
0
T
1
9
1
.
0
2
T
3
9
0
.
9
0
T
2
0
1
.
0
1
T
4
0
0
.
9
0
T
ab
le
4
.
Op
tim
u
m
R
ea
l P
o
w
er
L
o
s
s
Val
u
e
s
Ob
tain
ed
f
o
r
P
r
ac
tical
1
9
1
Utilit
y
(
I
n
d
ian
)
S
y
s
t
e
m
b
y
T
FP
Me
th
o
d
R
e
a
l
p
o
w
e
r
L
o
ss (M
W
)
TF
P
M
i
n
1
4
4
.
9
8
2
M
a
x
1
4
7
.
9
9
2
A
v
e
r
a
g
e
1
4
5
.
8
3
2
10.
CO
NCLU
SI
O
N
T
ailo
r
e
d
Flo
w
er
P
o
llin
atio
n
(
T
FP
)
alg
o
r
ith
m
h
a
s
b
ee
n
ef
f
ec
tiv
el
y
ap
p
lied
f
o
r
s
o
lv
i
n
g
r
ea
ctiv
e
p
o
w
e
r
p
r
o
b
lem
.
An
d
it
h
as
b
ee
n
te
s
ted
in
s
ta
n
d
ar
d
I
E
E
E
1
1
8
&
p
r
ac
tical
1
9
1
b
u
s
te
s
t
s
y
s
te
m
s
.
P
er
f
o
r
m
a
n
ce
co
m
p
ar
is
o
n
s
w
i
th
w
el
l
-
k
n
o
w
n
p
o
p
u
latio
n
-
b
ased
al
g
o
r
ith
m
s
g
iv
e
i
m
p
r
o
v
ed
r
es
u
lt
s
.
T
ailo
r
e
d
Flo
w
er
P
o
llin
atio
n
(
T
F
P
)
alg
o
r
ith
m
em
er
g
e
s
to
f
in
d
g
o
o
d
s
o
lu
tio
n
s
w
h
e
n
co
m
p
ar
ed
to
th
at
o
f
o
th
er
r
ep
o
r
ted
alg
o
r
ith
m
s
.
T
h
e
s
i
m
u
lat
io
n
r
es
u
lt
s
p
r
ese
n
ted
in
p
r
ev
io
u
s
s
ec
tio
n
p
r
o
v
e
t
h
e
ca
p
ab
ilit
y
o
f
T
FP
ap
p
r
o
ac
h
to
ar
r
i
v
e
at
n
ea
r
to
g
lo
b
al
o
p
tim
a
l so
lu
tio
n
.
RE
F
E
R
E
NC
E
S
[
1
]
O.
A
lsa
c
a
n
d
B
.
S
to
t
t,
"
Op
ti
m
a
l
L
o
a
d
F
lo
w
w
it
h
S
tea
d
y
-
S
tate
S
e
c
u
rit
y
,
"
in
IE
EE
T
ra
n
sa
c
ti
o
n
s
o
n
P
o
we
r A
p
p
a
r
a
tu
s
a
n
d
S
y
ste
ms
,
v
o
l.
P
A
S
-
9
3
,
n
o
.
3
,
p
p
.
7
4
5
-
7
5
1
,
M
a
y
1
9
7
4
.
[
2
]
K.
Y
.
L
e
e
,
Y
.
M
.
P
a
rk
a
n
d
J.
L
.
O
rti
z
,
"
A
Un
it
e
d
A
p
p
ro
a
c
h
t
o
O
p
ti
m
a
l
Re
a
l
a
n
d
Re
a
c
ti
v
e
P
o
w
e
r
Disp
a
tch
,
"
i
n
IE
EE
T
ra
n
s
a
c
ti
o
n
s o
n
P
o
we
r
Ap
p
a
r
a
tu
s a
n
d
S
y
ste
ms
,
v
o
l.
P
A
S
-
1
0
4
,
n
o
.
5
,
p
p
.
1
1
4
7
-
1
1
5
3
,
M
a
y
1
9
8
5
.
[
3
]
A
.
M
o
n
ti
c
e
ll
i,
M
.
V
.
F
.
P
e
re
ira an
d
S
.
G
ra
n
v
il
le,
"
S
e
c
u
rit
y
-
Co
n
stra
in
e
d
Op
ti
m
a
l
P
o
w
e
r
F
lo
w
w
it
h
P
o
st
-
Co
n
ti
n
g
e
n
c
y
Co
rre
c
ti
v
e
Re
sc
h
e
d
u
li
n
g
,
"
in
IE
E
E
T
ra
n
sa
c
ti
o
n
s
o
n
P
o
we
r S
y
ste
ms
,
v
o
l.
2
,
n
o
.
1
,
p
p
.
1
7
5
-
1
8
0
,
F
e
b
.
1
9
8
7
.
[
4
]
N.
De
e
b
a
n
d
S
.
M
.
S
h
a
h
id
e
h
p
o
u
r
,
"
L
in
e
a
r
re
a
c
ti
v
e
p
o
w
e
r
o
p
ti
m
iz
a
ti
o
n
in
a
lar
g
e
p
o
w
e
r
n
e
tw
o
rk
u
sin
g
th
e
d
e
c
o
m
p
o
siti
o
n
a
p
p
r
o
a
c
h
,
"
in
IEE
E
T
ra
n
sa
c
ti
o
n
s
o
n
P
o
we
r S
y
ste
ms
,
v
o
l.
5
,
n
o
.
2
,
p
p
.
4
2
8
-
4
3
8
,
M
a
y
1
9
9
0
.
[
5
]
E.
Hc
b
so
n
,
"
Ne
tw
o
rk
Co
n
stra
in
e
d
Re
a
c
ti
v
e
P
o
w
e
r
Co
n
tro
l
Us
in
g
L
in
e
a
r
P
ro
g
ra
m
m
in
g
,
"
in
IEE
E
T
ra
n
sa
c
ti
o
n
s
o
n
Po
we
r
Ap
p
a
r
a
tu
s
a
n
d
S
y
ste
ms
,
v
o
l.
P
A
S
-
9
9
,
n
o
.
3
,
p
p
.
8
6
8
-
8
7
7
,
M
a
y
1
9
8
0
.
[
6
]
K.
Y
.
L
e
e
,
Y
.
M
.
P
a
rk
a
n
d
J.
L
.
Ortiz
,
"
F
u
e
l
-
c
o
st
m
in
im
isa
ti
o
n
f
o
r
b
o
t
h
re
a
l
-
a
n
d
re
a
c
ti
v
e
-
p
o
w
e
r
d
isp
a
tch
e
s,"
in
IEE
Evaluation Warning : The document was created with Spire.PDF for Python.
I
n
t J
I
n
f
&
C
o
m
m
u
n
T
ec
h
n
o
l
I
SS
N:
2252
-
8776
Ta
ilo
r
ed
flo
w
er p
o
llin
a
tio
n
(
TFP
)
a
lg
o
r
ith
m
fo
r
d
imin
u
tio
n
o
f rea
l p
o
w
er lo
s
s
(
K
.
Len
in
)
101
Pr
o
c
e
e
d
in
g
s C
-
Ge
n
e
ra
ti
o
n
,
T
ra
n
sm
issio
n
a
n
d
Distri
b
u
ti
o
n
,
v
o
l.
1
3
1
,
n
o
.
3
,
p
p
.
8
5
-
9
3
,
M
a
y
1
9
8
4
.
[7
]
M
.
K.
M
a
n
g
o
l
i,
a
n
d
K.Y
.
L
e
e
,
“
Op
ti
m
a
l
re
a
l
a
n
d
re
a
c
ti
v
e
p
o
w
e
r
c
o
n
tr
o
l
u
sin
g
l
in
e
a
r
p
ro
g
ra
m
m
in
g
”,
El
e
c
tr
.
Po
we
r
S
y
st.
Res
.
, v
o
l
.
26
,
n
o
.
1
,
p
p
.
1
-
10
,
Ja
n
u
a
ry
1993
.
[
8
]
C.
A
.
Ca
n
iza
re
s,
A
.
C.
Z.
De
S
o
u
z
a
a
n
d
V
.
H.
Q
u
in
tan
a
,
"
Co
m
p
a
riso
n
o
f
p
e
rf
o
rm
a
n
c
e
in
d
ice
s
f
o
r
d
e
tec
ti
o
n
o
f
p
ro
x
im
it
y
to
v
o
lt
a
g
e
c
o
ll
a
p
se
,
"
in
IEE
E
T
r
a
n
sa
c
ti
o
n
s
o
n
Po
we
r
S
y
ste
ms
,
v
o
l.
1
1
,
n
o
.
3
,
p
p
.
1
4
4
1
-
1
4
5
0
,
A
u
g
.
1
9
9
6
.
[9
]
S
.
R.
P
a
ra
n
j
o
th
i
,
a
n
d
K
.
A
n
b
u
ra
ja,
“
Op
ti
m
a
l
p
o
w
e
r
f
lo
w
u
sin
g
re
f
in
e
d
g
e
n
e
ti
c
a
lg
o
rit
h
m
”
,
El
e
c
tr
.
Po
we
r
Co
mp
o
n
.
S
y
st
,
v
o
l.
3
0
,
n
o
.
1
0
,
p
p
.
1
0
5
5
-
1
0
6
3
,
2
0
0
2
.
[
10
]
D.
De
v
a
ra
j
a
n
d
B.
Y
e
g
n
a
n
a
ra
y
a
n
a
,
"
G
e
n
e
ti
c
-
a
lg
o
rit
h
m
-
b
a
se
d
o
p
ti
m
a
l
p
o
w
e
r
f
lo
w
f
o
r
se
c
u
rit
y
e
n
h
a
n
c
e
m
e
n
t,
"
in
IEE
Pr
o
c
e
e
d
in
g
s
-
Ge
n
e
ra
ti
o
n
,
T
ra
n
s
miss
io
n
a
n
d
Distrib
u
ti
o
n
,
v
o
l.
1
5
2
,
n
o
.
6
,
p
p
.
8
9
9
-
9
0
5
,
4
N
o
v
.
2
0
0
5
.
[
1
1
]
A
.
Be
rizz
i,
C.
Bo
v
o
,
M
.
M
e
rl
o
,
a
n
d
M
.
De
lf
a
n
ti
,
“
A
GA
a
p
p
ro
a
c
h
to
c
o
m
p
a
re
o
rp
f
o
b
jec
ti
v
e
f
u
n
c
ti
o
n
s
i
n
c
lu
d
in
g
se
c
o
n
d
a
ry
v
o
lt
a
g
e
re
g
u
latio
n
,
”
El
e
c
tric P
o
we
r S
y
ste
ms
Res
e
a
r
c
h
,
v
o
l.
8
4
,
n
o
.
1
,
p
p
.
1
8
7
–
1
9
4
,
2
0
1
2
.
[
12
]
C.
-
F
.
Y
a
n
g
,
G
.
G
.
L
a
i,
C.
-
H.
L
e
e
,
C.
-
T
.
S
u
,
a
n
d
G
.
W
.
Ch
a
n
g
,
“
Op
ti
m
a
l
se
tt
in
g
o
f
re
a
c
ti
v
e
c
o
m
p
e
n
sa
ti
o
n
d
e
v
ice
s
w
it
h
a
n
im
p
ro
v
e
d
v
o
lt
a
g
e
sta
b
il
it
y
in
d
e
x
f
o
r
v
o
lt
a
g
e
sta
b
il
it
y
e
n
h
a
n
c
e
m
e
n
t,
”
In
ter
n
a
ti
o
n
a
l
J
o
u
rn
a
l
o
f
E
l
e
c
trica
l
Po
we
r
a
n
d
En
e
r
g
y
S
y
ste
ms
,
v
o
l.
3
7
,
n
o
.
1
,
p
p
.
5
0
–
5
7
,
2
0
1
2
.
[1
3
]
P
.
Ro
y
,
S
.
G
h
o
sh
a
l,
a
n
d
S
.
T
h
a
k
u
r
,
“
Op
ti
m
a
l
v
a
r
c
o
n
tro
l
f
o
r
im
p
ro
v
e
m
e
n
ts
in
v
o
lt
a
g
e
p
ro
f
il
e
s
a
n
d
f
o
r
re
a
l
p
o
w
e
r
lo
ss
m
in
i
m
iza
ti
o
n
u
sin
g
b
i
o
g
e
o
g
ra
p
h
y
b
a
se
d
o
p
ti
m
iza
ti
o
n
,
”
In
ter
n
a
ti
o
n
a
l
J
o
u
rn
a
l
o
f
E
lec
trica
l
P
o
w
e
r
a
n
d
E
n
e
r
g
y
S
y
ste
ms
,
v
o
l.
4
3
,
n
o
.
1
,
p
p
.
8
3
0
–
8
3
8
,
2
0
1
2
.
[
14
]
B.
V
e
n
k
a
tes
h
,
G
.
S
a
d
a
siv
a
m
a
n
d
M
.
A
.
Kh
a
n
,
"
A
n
e
w
o
p
ti
m
a
l
re
a
c
ti
v
e
p
o
we
r
sc
h
e
d
u
li
n
g
m
e
th
o
d
f
o
r
lo
s
s
m
in
i
m
iza
ti
o
n
a
n
d
v
o
lt
a
g
e
sta
b
il
it
y
m
a
r
g
in
m
a
x
i
m
iz
a
ti
o
n
u
sin
g
su
c
c
e
ss
i
v
e
m
u
lt
i
-
o
b
jec
ti
v
e
f
u
z
z
y
LP
tec
h
n
iq
u
e
,
"
in
IEE
E
T
r
a
n
sa
c
ti
o
n
s
o
n
Po
we
r
S
y
ste
ms
,
v
o
l.
1
5
,
n
o
.
2
,
p
p
.
8
4
4
-
8
5
1
,
M
a
y
2
0
0
0
.
[
15
]
W
e
i
Y
a
n
,
S
h
u
a
i
L
u
a
n
d
D
.
C.
Y
u
,
"
A
n
o
v
e
l
o
p
ti
m
a
l
re
a
c
ti
v
e
p
o
we
r
d
isp
a
tch
m
e
th
o
d
b
a
se
d
o
n
a
n
im
p
ro
v
e
d
h
y
b
ri
d
e
v
o
lu
ti
o
n
a
ry
p
ro
g
ra
m
m
in
g
tec
h
n
iq
u
e
,
"
in
IEE
E
T
ra
n
s
a
c
ti
o
n
s
o
n
P
o
we
r
S
y
ste
ms
,
v
o
l.
1
9
,
n
o
.
2
,
p
p
.
9
1
3
-
9
1
8
,
M
a
y
2
0
0
4
.
[
16
]
W
e
i
Y
a
n
,
F
a
n
g
L
iu
,
C.
Y
.
C
h
u
n
g
a
n
d
K.
P
.
W
o
n
g
,
"
A
h
y
b
rid
g
e
n
e
ti
c
a
lg
o
rit
h
m
-
in
terio
r
p
o
i
n
t
m
e
th
o
d
f
o
r
o
p
ti
m
a
l
re
a
c
ti
v
e
p
o
we
r
f
lo
w
,
"
in
IEE
E
T
r
a
n
sa
c
ti
o
n
s
o
n
Po
we
r
S
y
ste
ms
,
v
o
l.
2
1
,
n
o
.
3
,
p
p
.
1
1
6
3
-
1
1
6
9
,
A
u
g
.
2
0
0
6
.
[
17
]
J.
Y
u
,
W
.
Y
a
n
,
W
.
L
i,
C.
Y
.
Ch
u
n
g
a
n
d
K.
P
.
W
o
n
g
,
"
A
n
Un
f
ix
e
d
P
ie
c
e
w
ise
-
Op
ti
m
a
l
R
e
a
c
ti
v
e
P
o
w
e
r
-
F
lo
w
M
o
d
e
l
a
n
d
it
s
A
lg
o
rit
h
m
f
o
r
A
C
-
DC S
y
ste
m
s
,
"
in
IEE
E
T
ra
n
sa
c
ti
o
n
s
o
n
P
o
we
r S
y
ste
ms
,
v
o
l.
2
3
,
n
o
.
1
,
p
p
.
1
7
0
-
1
7
6
,
F
e
b
.
2
0
0
8
.
[
18
]
F
.
Ca
p
it
a
n
e
sc
u
,
"
A
ss
e
s
sin
g
Re
a
c
t
iv
e
P
o
w
e
r
Re
se
r
v
e
s
W
it
h
Re
sp
e
c
t
to
Op
e
ra
ti
n
g
Co
n
stra
in
ts
a
n
d
V
o
l
tag
e
S
tab
il
it
y
,
"
in
IEE
E
T
ra
n
sa
c
ti
o
n
s
o
n
Po
we
r
S
y
ste
ms
,
v
o
l.
2
6
,
n
o
.
4
,
p
p
.
2
2
2
4
-
2
2
3
4
,
No
v
.
2
0
1
1
.
[
19
]
Z.
Hu
,
X.
W
a
n
g
,
a
n
d
G
.
T
a
y
lo
r
,
“
S
to
c
h
a
stic
o
p
t
im
a
l
re
a
c
ti
v
e
p
o
we
r
d
isp
a
tch
:
F
o
rm
u
latio
n
a
n
d
so
l
u
ti
o
n
m
e
th
o
d
,
”
In
ter
n
a
t
io
n
a
l
J
o
u
rn
a
l
o
f
E
lec
trica
l
Po
we
r a
n
d
En
e
r
g
y
S
y
s
tem
s
,
v
o
l.
3
2
,
n
o
.
6
,
p
p
.
6
1
5
–
6
2
1
,
Ju
ly
2010.
[
20
]
A
.
Ka
r
g
a
rian
,
M
.
Ra
o
o
f
a
t,
a
n
d
M
.
M
o
h
a
m
m
a
d
i,
“
P
r
o
b
a
b
i
li
stic
r
e
a
c
ti
v
e
p
o
we
r
p
ro
c
u
re
m
e
n
t
in
h
y
b
rid
e
lec
tri
c
it
y
m
a
rk
e
ts
w
it
h
u
n
c
e
rtain
l
o
a
d
s,
”
El
e
c
tric P
o
we
r S
y
ste
ms
Res
e
a
r
c
h
,
v
o
l.
8
2
,
n
o
.
1
,
p
p
.
6
8
–
8
0
,
Ja
n
u
a
ry
2
0
1
2
.
[
21
]
Y
a
n
g
,
X
.
S
.
“
F
lo
we
r
p
o
ll
i
n
a
t
io
n
a
lg
o
rith
m
fo
r
g
lo
b
a
l
o
p
ti
miza
ti
o
n
,
i
n
:
Un
c
o
n
v
e
n
ti
o
n
a
l
C
o
mp
u
ta
t
i
o
n
a
n
d
N
a
tu
r
a
l
Co
mp
u
t
a
ti
o
n
,
”
I
n
:
D
u
ra
n
d
-
L
o
se
J
.
,
Jo
n
o
sk
a
N.
(e
d
s)
Un
c
o
n
v
e
n
ti
o
n
a
l
Co
m
p
u
tatio
n
a
n
d
Na
tu
ra
l
Co
m
p
u
tati
o
n
.
UCN
C
2
0
1
2
.
L
e
c
tu
re
No
tes
i
n
Co
m
p
u
ter
S
c
ien
c
e
,
V
o
l.
7
4
4
5
,
p
p
.
2
4
0
–
2
4
9
,
2
0
1
2
.
[2
2
]
IE
EE
,
“
T
h
e
I
EE
E
3
0
-
b
u
s tes
t
sy
ste
m a
n
d
th
e
IEE
E
1
1
8
-
tes
t
sy
ste
m
”
,
(1
9
9
3
)
,
h
tt
p
:
//
ww
w
.
e
e
.
w
a
sh
in
g
to
n
.
e
d
u
/
trse
a
rc
h
/p
stc
a
/.
[2
3
]
Jia
n
g
tao
Ca
o
,
F
u
li
W
a
n
g
a
n
d
P
i
n
g
L
i,
“
A
n
Im
p
ro
v
e
d
Bio
g
e
o
g
ra
p
h
y
-
b
a
se
d
Op
ti
m
iza
ti
o
n
A
lg
o
rit
h
m
f
o
r
Op
ti
m
a
l
Re
a
c
ti
v
e
P
o
w
e
r
F
lo
w
”
,
In
ter
n
a
ti
o
n
a
l
J
o
u
rn
a
l
o
f
C
o
n
tr
o
l
a
n
d
A
u
to
m
a
ti
o
n
, v
o
l
.
7
,
n
o
.
3
,
p
p
.
1
6
1
-
1
7
6
,
M
a
rc
h
2
0
1
4
.
Evaluation Warning : The document was created with Spire.PDF for Python.