I
n
d
o
n
e
s
i
a
n
J
o
u
r
n
a
l
o
f
E
l
e
c
t
r
i
c
a
l
E
n
g
i
n
e
e
r
i
n
g
a
n
d
C
o
m
p
u
t
e
r
S
c
i
e
n
c
e
V
o
l
.
11
,
N
o
.
3
,
S
e
p
t
e
m
b
e
r
201
8
,
p
p
.
1109
~
1113
I
S
S
N
:
2
5
0
2
-
4
7
5
2
,
D
O
I
:
1
0
.
1
1
5
9
1
/
i
j
e
e
c
s
.
v
1
1
.
i
3
.
p
p
1
1
0
9
-
1
1
1
3
1109
J
o
u
r
n
a
l
h
o
m
e
p
a
g
e
:
h
t
t
p
:
/
/
i
a
e
s
c
o
r
e
.
c
o
m
/
j
o
u
r
n
a
l
s
/
i
n
d
e
x
.
p
h
p
/
i
j
e
e
c
s
C
u
b
e
A
r
i
th
m
e
ti
c
:
I
m
p
r
o
v
i
n
g
Eu
l
e
r
M
e
th
o
d
fo
r
O
r
d
i
n
a
r
y
D
i
ffe
r
e
n
ti
a
l
Eq
u
a
ti
o
n
U
s
i
n
g
C
u
b
e
M
e
a
n
N
o
o
r
a
i
d
a
S
a
m
s
u
d
i
n
,
N
u
r
h
a
f
i
z
a
h
M
o
z
i
y
a
n
a
M
o
h
d
Y
u
s
o
p
,
S
y
a
h
r
u
l
F
a
h
m
y
,
A
n
i
s
S
h
a
h
i
d
a
N
i
z
a
b
i
n
t
i
M
o
k
h
t
a
r
T
A
T
I U
ni
ve
rs
i
t
y
Col
l
e
g
e
, T
e
l
uk K
a
l
ong
,
K
e
m
a
m
a
n,
T
e
re
ng
g
a
nu
,
U
ni
ve
rs
i
t
i
P
e
rt
a
ha
na
n N
a
s
i
ona
l
M
a
l
a
y
s
i
a
, K
e
m
S
ung
a
i
Be
s
i
, K
ua
l
a
L
um
p
ur
,
M
a
l
a
y
s
i
a
A
r
t
i
c
l
e
I
n
f
o
A
B
S
T
R
A
C
T
A
r
t
i
c
l
e
h
i
s
t
o
r
y
:
R
e
c
e
i
v
e
d
M
a
y
1
,
2
0
1
8
R
e
v
i
s
e
d
J
u
n
2
,
2
0
1
8
A
c
c
e
p
t
e
d
J
u
n
11
,
2
0
1
8
T
he
E
ul
e
r
m
e
t
hod
i
s
a
fi
rs
t
-
orde
r
num
e
ri
c
a
l
p
roc
e
dure
for
s
ol
vi
ng
O
rdi
na
ry
D
i
ffe
re
nt
i
a
l
E
qua
t
i
on
(O
D
E
s
)
p
robl
e
m
s
. It
i
s
a
n e
ffe
c
t
i
ve
a
nd e
a
s
y
m
e
t
hod t
o
s
ol
ve
i
ni
t
i
a
l
va
l
ue
p
robl
e
m
s
.
A
l
t
houg
h
E
ul
e
r
p
rovi
de
s
s
i
m
p
l
e
p
roc
e
dure
for
s
ol
vi
ng
O
D
E
s
,
t
he
re
ha
ve
be
e
n
i
s
s
ue
s
s
uc
h a
s
c
om
p
l
e
x
i
t
y
, t
i
m
e
of p
roc
e
s
s
i
ng
a
nd
a
c
c
ura
c
y
t
ha
t
c
om
p
e
l
l
e
d
t
he
us
e
of
ot
he
r,
m
ore
c
om
p
l
e
x
,
m
e
t
hods
.
Im
p
rove
m
e
nt
s
t
o
t
he
E
ul
e
r
m
e
t
hod ha
ve
a
t
t
ra
c
t
e
d m
uc
h a
t
t
e
nt
i
on re
s
ul
t
i
ng
i
n
num
e
rous
m
odi
fi
e
d
E
ul
e
r
m
e
t
hods
.
T
hi
s
p
a
p
e
r
p
rop
os
e
s
Cube
A
ri
t
hm
e
t
i
c
, a
m
odi
fi
e
d
E
ul
e
r
m
e
t
hod
w
i
t
h
i
m
p
rove
d
a
c
c
ura
c
y
.
T
he
e
ffi
c
i
e
nc
y
of
Cube
A
ri
t
hm
e
t
i
c
w
a
s
c
om
p
a
re
d
w
i
t
h
E
ul
e
r
A
ri
t
hm
e
t
i
c
a
nd
t
e
s
t
e
d
us
i
ng
S
CIL
A
B
a
g
a
i
ns
t
e
x
a
c
t
s
ol
ut
i
ons
.
Re
s
ul
t
s
i
ndi
c
a
t
e
t
ha
t
not
onl
y
Cube
A
ri
t
hm
e
t
i
c
p
rovi
de
d
s
ol
ut
i
ons
t
ha
t
a
re
s
i
m
i
l
a
r
t
o
e
x
a
c
t
s
ol
ut
i
ons
a
t
s
m
a
l
l
s
t
e
p
s
i
z
e
,
but
a
l
s
o a
t
hi
g
he
r
s
t
e
p
s
i
z
e
, he
nc
e
p
roduc
i
ng
m
ore
a
c
c
ura
t
e
re
s
ul
t
s
.
K
e
y
w
o
r
d
s
:
C
u
b
e
a
r
i
t
h
m
e
t
i
c
E
u
l
e
r
S
C
I
L
A
B
Copy
r
i
ght
©
201
8
Ins
t
i
t
ut
e
of
A
dv
anc
e
d E
ngi
ne
e
r
i
ng
and Sc
i
e
nc
e
.
A
l
l
r
i
ght
s
r
e
s
e
r
v
e
d
.
C
o
r
r
e
s
p
o
n
d
i
n
g
A
u
t
h
o
r
:
N
o
o
r
a
i
d
a
S
a
m
s
u
d
i
n
,
T
A
T
I
U
n
i
v
e
r
s
i
t
y
C
o
l
l
e
g
e
,
T
e
l
u
k
K
a
l
o
n
g
,
K
e
m
a
m
a
n
,
T
e
r
e
n
g
g
a
n
u
,
U
n
i
v
e
r
s
i
t
i
P
e
r
t
a
h
a
n
a
n
N
a
s
i
o
n
a
l
M
a
l
a
y
s
i
a
,
K
e
m
S
u
n
g
a
i
B
e
s
i
,
K
u
a
l
a
L
u
m
p
u
r
,
M
a
l
a
y
s
i
a
.
E
m
a
i
l
:
n
o
o
r
a
i
d
a
@t
a
t
i
u
c
.
e
d
u
.
m
y
1.
I
N
T
R
O
D
U
C
T
I
O
N
N
u
m
e
r
i
c
a
l
m
e
t
h
o
d
i
s
w
i
d
e
l
y
u
s
e
d
i
n
c
o
m
p
u
t
e
r
s
c
i
e
n
c
e
.
I
t
h
e
l
p
s
r
e
s
e
a
r
c
h
e
r
s
,
e
n
g
i
n
e
e
r
s
,
a
n
d
s
t
u
d
e
n
t
s
t
o
a
p
p
r
o
x
i
m
a
t
e
t
h
e
s
o
l
u
t
i
o
n
o
f
c
o
m
p
l
i
c
a
t
e
d
p
r
o
b
l
e
m
s
s
o
t
h
a
t
t
h
e
f
i
n
a
l
s
o
l
u
t
i
o
n
c
o
n
s
i
s
t
s
o
f
o
n
l
y
a
d
d
i
t
i
o
n
,
s
u
b
t
r
a
c
t
i
o
n
a
n
d
m
u
l
t
i
p
l
i
c
a
t
i
o
n
o
p
e
r
a
t
i
o
n
s
.
N
u
m
e
r
i
c
a
l
m
e
t
h
o
d
f
o
r
m
s
a
n
i
m
p
o
r
t
a
n
t
p
a
r
t
o
f
s
o
l
v
i
n
g
I
n
i
t
i
a
l
V
a
l
u
e
P
r
o
b
l
e
m
s
(
I
V
P
s
)
i
n
O
r
d
i
n
a
r
y
D
i
f
f
e
r
e
n
t
i
a
l
E
q
u
a
t
i
o
n
s
(
O
D
E
s
)
,
e
s
p
e
c
i
a
l
l
y
i
n
c
a
s
e
s
w
h
e
r
e
t
h
e
e
x
a
c
t
s
o
l
u
t
i
o
n
i
s
d
i
f
f
i
c
u
l
t
t
o
o
b
t
a
i
n
o
r
i
n
t
h
e
a
b
s
e
n
c
e
o
f
c
l
o
s
e
d
f
o
r
m
a
n
a
l
y
t
i
c
f
o
r
m
u
l
a
.
T
h
e
E
u
l
e
r
m
e
t
h
o
d
i
s
a
f
i
r
s
t
-
o
r
d
e
r
n
u
m
e
r
i
c
a
l
p
r
o
c
e
d
u
r
e
f
o
r
s
o
l
v
i
n
g
O
D
E
s
w
i
t
h
a
g
i
v
e
n
i
n
i
t
i
a
l
v
a
l
u
e
.
I
t
i
s
a
n
e
f
f
e
c
t
i
v
e
[
1
]
a
n
d
t
h
e
e
a
s
i
e
s
t
m
e
t
h
o
d
[
2
]
t
o
s
o
l
v
e
I
V
P
s
i
n
O
D
E
s
.
E
u
l
e
r
m
e
t
h
o
d
h
a
s
b
e
e
n
u
s
e
i
n
v
a
r
i
e
t
y
f
i
e
l
d
s
u
c
h
a
s
i
n
c
o
m
m
u
n
i
c
a
t
i
o
n
a
n
d
s
e
c
u
r
i
t
y
t
o
g
e
t
b
e
t
t
e
r
s
o
l
u
t
i
o
n
.
R
e
g
a
r
d
i
n
g
[
3
]
,
E
u
l
e
r
m
e
t
h
o
d
h
a
s
b
e
t
t
e
r
p
e
r
f
o
r
m
a
n
c
e
i
n
m
o
d
e
l
i
n
g
u
n
c
e
r
t
a
i
n
t
i
e
s
.
I
t
i
s
s
u
p
p
e
o
t
e
d
b
y
[
4
]
w
h
e
r
e
a
s
e
x
p
e
r
i
m
e
n
t
a
l
r
e
s
u
l
t
s
i
n
d
i
c
a
t
e
t
h
a
t
i
t
i
s
m
u
c
h
s
i
m
p
l
e
r
a
n
d
f
a
s
t
e
r
,
m
o
r
e
s
e
n
s
i
t
i
v
e
a
n
d
r
o
b
u
s
t
.
A
l
t
h
o
u
g
h
E
u
l
e
r
p
r
o
v
i
d
e
s
s
i
m
p
l
e
p
r
o
c
e
d
u
r
e
f
o
r
s
o
l
v
i
n
g
O
D
E
s
,
t
h
e
r
e
a
r
e
i
s
s
u
e
s
t
h
a
t
h
a
v
e
c
o
m
p
e
l
l
e
d
t
h
e
u
s
e
o
f
o
t
h
e
r
,
m
o
r
e
c
o
m
p
l
e
x
,
m
e
t
h
o
d
s
s
u
c
h
a
s
R
u
n
g
e
K
u
t
t
e
.
O
n
e
o
f
t
h
e
m
i
s
“
a
c
c
u
r
a
c
y
”
a
s
r
e
p
o
r
t
e
d
b
y
[
5
-
6
]
.
I
m
p
r
o
v
e
m
e
n
t
s
t
o
t
h
e
E
u
l
e
r
m
e
t
h
o
d
h
a
v
e
a
t
t
r
a
c
t
e
d
m
u
c
h
a
t
t
e
n
t
i
o
n
i
n
c
l
u
d
i
n
g
[
1
-
2
]
,
[
5
-
7
]
,
r
e
s
u
l
t
i
n
g
i
n
m
o
d
i
f
i
e
d
E
u
l
e
r
m
e
t
h
o
d
s
.
T
h
i
s
p
a
p
e
r
i
s
a
i
m
t
o
d
i
s
c
o
v
e
r
a
n
e
w
a
l
g
o
r
i
t
h
m
a
s
a
c
c
u
r
a
t
e
a
s
p
o
s
s
i
b
l
e
w
i
t
h
e
x
a
c
t
s
o
l
u
t
i
o
n
b
y
u
s
i
n
g
a
v
e
r
a
g
e
c
o
n
c
e
p
t
.
T
h
e
n
e
w
a
l
g
o
r
i
t
h
m
p
r
o
p
o
s
e
d
i
n
t
h
i
s
p
a
p
e
r
n
a
m
e
i
t
a
s
C
u
b
e
A
r
i
t
h
m
e
t
i
c
.
T
h
i
s
i
d
e
a
t
o
d
e
v
e
l
o
p
t
h
i
s
n
e
w
a
l
g
o
r
i
t
h
m
i
s
b
y
e
n
h
a
n
c
e
o
t
h
e
r
a
l
g
o
r
i
t
h
m
p
r
o
p
o
s
e
d
b
y
Q
u
r
e
s
h
i
[
7
]
a
n
d
N
u
r
h
a
f
i
z
a
h
[
8
]
.
Q
u
r
e
s
h
i
u
s
e
d
a
v
e
r
a
g
e
o
f
a
r
i
t
h
m
e
t
i
c
m
e
a
n
f
o
r
t
w
o
p
o
i
n
t
s
a
t
c
o
o
r
d
i
n
a
t
e
y
=
[
f
(
x
n
,
y
n
)
+
f
(
x
n
+
1
,
y
n
+
1
)
]
/
2
.
T
h
e
e
q
u
a
t
i
o
n
i
s
r
e
f
e
r
r
i
n
g
a
s
E
u
l
e
r
A
r
i
t
h
m
e
t
i
c
s
h
o
w
s
i
m
p
r
o
v
e
m
e
n
t
i
n
a
c
c
u
r
a
c
y
a
n
d
s
p
e
e
d
c
o
m
p
a
r
e
d
t
o
E
u
l
e
r
M
e
t
h
o
d
.
[
8
]
u
s
e
s
c
o
n
c
e
p
t
o
f
Q
u
r
e
s
h
i
[
7
]
b
u
t
c
h
o
o
s
e
H
a
r
m
o
n
i
c
m
e
a
n
a
n
d
C
o
n
t
r
a
h
a
r
m
o
n
i
c
m
e
a
n
i
n
t
h
e
e
q
u
a
t
i
o
n
.
Evaluation Warning : The document was created with Spire.PDF for Python.
I
S
S
N
:
2502
-
4
7
5
2
I
n
d
o
n
e
s
i
a
n
J
E
l
e
c
E
n
g
&
C
o
m
p
S
c
i
,
V
o
l
.
11
,
N
o
.
3
,
S
e
p
t
e
m
b
e
r
2
0
1
8
:
1109
–
1113
1110
T
h
u
s
,
t
h
i
s
p
a
p
e
r
p
r
o
p
o
s
e
s
a
m
o
d
i
f
i
e
d
E
u
l
e
r
m
e
t
h
o
d
u
s
i
n
g
t
h
e
C
u
b
e
m
e
a
n
c
o
n
c
e
p
t
.
T
h
e
m
e
t
h
o
d
w
a
s
c
o
n
s
t
r
u
c
t
e
d
b
y
e
x
t
r
a
c
t
i
n
g
t
h
e
E
u
l
e
r
A
r
i
t
h
m
e
t
i
c
(
E
A
)
m
e
t
h
o
d
i
n
[
9
]
a
n
d
u
s
i
n
g
C
u
b
e
m
e
a
n
i
n
t
h
e
e
q
u
a
t
i
o
n
.
P
r
o
p
o
s
e
a
l
g
o
r
i
t
h
m
i
s
c
o
m
p
a
r
e
d
b
e
t
w
e
e
n
m
o
d
i
f
i
e
d
E
u
l
e
r
k
n
o
w
n
a
s
E
u
l
e
r
A
r
i
t
h
m
e
t
i
c
w
i
t
h
t
h
e
e
x
a
c
t
s
o
l
u
t
i
o
n
.
T
h
e
a
c
c
u
r
a
c
y
o
f
t
h
e
p
r
o
p
o
s
e
d
m
e
t
h
o
d
i
s
d
e
m
o
n
s
t
r
a
t
e
d
b
y
t
h
e
c
o
m
p
a
r
i
n
g
e
r
r
o
r
r
e
s
u
l
t
s
o
f
w
i
t
h
s
t
e
p
s
i
z
e
s
o
f
h
0
.
1
,
0
.
0
1
a
n
d
0
.
0
0
1
.
2.
P
R
O
P
O
S
E
D
A
L
G
O
R
I
T
H
M
T
h
e
a
l
g
o
r
i
t
h
m
i
s
b
a
s
e
d
o
n
t
h
e
C
u
b
e
M
e
a
n
i
n
E
u
l
e
r
A
r
i
t
h
m
e
t
i
c
.
E
u
l
e
r
A
r
i
t
h
m
e
t
i
c
(
E
A
)
i
s
a
w
e
l
l
-
k
n
o
w
n
t
e
c
h
n
i
q
u
e
f
o
r
i
m
p
r
o
v
i
n
g
t
h
e
E
u
l
e
r
m
e
t
h
o
d
.
E
u
l
e
r
m
e
t
h
o
d
s
u
s
e
d
i
n
[
1
]
a
n
d
[
2
]
w
e
r
e
e
x
t
r
a
c
t
e
d
a
s
b
a
s
i
s
f
o
r
t
h
e
a
l
g
o
r
i
t
h
m
.
C
u
b
e
m
e
a
n
w
i
t
h
i
n
t
w
o
c
o
o
r
d
i
n
a
t
e
p
o
i
n
t
s
o
f
f
u
n
c
t
i
o
n
w
e
r
e
u
s
e
d
t
o
i
m
p
r
o
v
e
t
h
e
m
e
t
h
o
d
s
.
E
q
u
a
t
i
o
n
(
1
)
i
s
a
b
a
s
i
c
f
o
r
m
u
l
a
o
f
t
h
e
E
u
l
e
r
m
e
t
h
o
d
.
)
,
(
0
0
1
n
y
t
tf
y
y
n
(
1
)
E
q
u
a
t
i
o
n
(
2
)
b
y
[
1
]
w
a
s
m
o
d
i
f
i
e
d
u
s
i
n
g
t
h
e
a
v
e
r
a
g
e
c
o
n
c
e
p
t
.
)
2
)
,
(
(
,
(
2
,
2
(
,
1
n
n
n
n
n
n
n
n
n
n
y
x
hf
y
h
x
f
y
x
f
h
y
h
x
hf
y
y
(
2
)
T
h
e
p
r
o
p
o
s
e
d
a
v
e
r
a
g
e
i
n
c
u
b
e
m
e
a
n
f
o
r
t
h
e
t
w
o
m
i
d
p
o
i
n
t
s
i
s
w
r
i
t
t
e
n
a
s
E
q
u
a
t
i
o
n
(
3
)
3
3
3
,
1
)
2
))
,
(
(
,
(
(
2
,
2
(
n
n
n
n
n
n
n
n
n
n
y
x
hf
y
h
x
f
y
x
f
h
y
h
x
hf
y
y
(
4
)
M
o
d
i
f
i
c
a
t
i
o
n
s
t
o
t
h
e
s
l
o
p
e
o
f
t
h
e
f
u
n
c
t
i
o
n
a
t
e
s
t
i
m
a
t
e
d
m
i
d
p
o
i
n
t
s
o
f
0
0
,
y
x
a
n
d
1
1
,
y
x
w
o
u
l
d
i
m
p
r
o
v
e
t
h
e
s
t
a
b
i
l
i
t
y
a
n
d
a
c
c
u
r
a
c
y
o
f
E
u
l
e
r
.
3.
R
E
S
E
A
R
C
H
M
E
T
H
O
D
T
h
e
r
e
a
r
e
a
f
e
w
p
r
o
c
e
s
s
e
s
t
o
p
r
e
s
e
n
t
a
m
a
t
h
e
m
a
t
i
c
a
l
m
o
d
e
l
.
F
i
r
s
t
,
i
t
c
a
n
b
e
p
r
e
s
e
n
t
e
d
b
y
m
a
t
h
e
m
a
t
i
c
a
l
s
o
f
t
w
a
r
e
a
n
d
a
l
g
o
r
i
t
h
m
d
e
v
e
l
o
p
m
e
n
t
.
B
e
f
o
r
e
u
s
i
n
g
t
h
e
m
a
t
h
e
m
a
t
i
c
a
l
s
o
f
t
w
a
r
e
,
a
l
g
o
r
i
t
h
m
d
e
v
e
l
o
p
m
e
n
t
m
u
s
t
b
e
d
o
n
e
f
i
r
s
t
.
A
l
g
o
r
i
t
h
m
r
e
f
e
r
s
a
s
i
n
s
t
r
u
c
t
i
o
n
s
f
o
r
s
o
l
v
i
n
g
p
r
o
b
l
e
m
u
s
i
n
g
s
i
m
p
l
e
l
a
n
g
u
a
g
e
[
1
0
-
1
1
]
.
A
c
c
o
r
d
i
n
g
t
o
[
1
2
]
,
c
o
n
s
t
r
u
c
t
i
o
n
p
r
o
c
e
s
s
i
n
m
a
t
h
e
m
a
t
i
c
a
l
s
o
f
t
w
a
r
e
i
n
c
l
u
d
e
s
a
s
f
o
l
l
o
w
s
.
a)
T
h
e
d
e
s
i
g
n
a
n
d
a
n
a
l
y
s
i
s
a
l
g
o
r
i
t
h
m
s
b)
A
l
g
o
r
i
t
h
m
c
o
d
i
n
g
c)
T
e
s
t
i
n
g
d)
D
e
t
a
i
l
s
d
o
c
u
m
e
n
t
a
t
i
o
n
e)
D
i
s
t
r
i
b
u
t
i
o
n
a
n
d
m
a
i
n
t
e
n
a
n
c
e
o
f
t
h
e
s
o
f
t
w
a
r
e
A
l
g
o
r
i
t
h
m
t
h
a
t
h
a
s
b
e
e
n
d
e
v
e
l
o
p
e
d
w
i
l
l
b
e
c
o
n
v
e
r
t
i
n
g
t
o
t
h
e
p
r
o
g
r
a
m
m
i
n
g
l
a
n
g
u
a
g
e
.
I
n
t
h
i
s
s
t
u
d
y
,
t
h
e
c
o
d
e
p
r
o
g
r
a
m
i
s
w
r
i
t
t
e
n
i
n
S
C
I
L
A
B
6
.
0
t
o
t
e
s
t
t
h
e
e
f
f
e
c
t
i
v
e
n
e
s
s
o
f
t
h
e
a
l
g
o
r
i
t
h
m
.
T
h
e
e
f
f
e
c
t
i
v
e
n
e
s
s
o
f
t
h
e
a
l
g
o
r
i
t
h
m
w
i
l
l
b
e
i
d
e
n
t
i
f
y
i
n
g
b
y
t
a
k
i
n
g
i
t
a
c
c
u
r
a
c
y
a
n
d
t
i
m
e
o
f
t
a
k
e
n
t
o
p
r
o
c
e
s
s
t
h
e
a
l
g
o
r
i
t
h
m
.
T
h
i
s
p
a
p
e
r
w
i
l
l
f
o
c
u
s
o
n
a
c
c
u
r
a
c
y
b
y
c
o
m
p
a
r
i
n
g
t
h
e
e
x
a
c
t
s
o
l
u
t
i
o
n
a
n
d
s
o
l
u
t
i
o
n
b
e
t
w
e
e
n
p
r
o
p
o
s
e
d
m
e
t
h
o
d
k
n
o
w
n
a
s
C
u
b
e
E
u
l
e
r
A
r
i
t
h
m
e
t
i
c
a
n
d
m
o
d
i
f
i
e
d
E
u
l
e
r
n
a
m
i
n
g
a
s
A
r
i
t
h
m
e
t
i
c
E
u
l
e
r
.
C
u
b
e
A
r
i
t
h
m
e
t
i
c
(
C
A
)
,
t
h
e
p
r
o
p
o
s
e
d
a
l
g
o
r
i
t
h
m
i
n
E
q
u
a
t
i
o
n
(
3
)
,
w
a
s
r
e
p
r
e
s
e
n
t
e
d
u
s
i
n
g
S
C
I
L
A
B
6
.
0
.
S
C
I
L
A
B
i
s
a
p
o
w
e
r
f
u
l
,
l
o
w
-
c
o
s
t
[
1
3
]
-
[
1
4
]
a
n
d
f
l
e
x
i
b
l
e
t
o
o
l
f
o
r
r
e
s
e
a
r
c
h
e
r
s
,
e
n
g
i
n
e
e
r
s
,
a
n
d
s
t
u
d
e
n
t
s
t
o
u
s
e
i
n
m
a
t
h
e
m
a
t
i
c
a
l
c
o
m
p
u
t
e
r
a
p
p
l
i
c
a
t
i
o
n
.
S
c
i
l
a
b
w
a
s
d
e
v
e
l
o
p
e
d
b
y
t
h
e
F
r
e
n
c
h
N
a
t
i
o
n
a
l
I
n
s
t
i
t
u
t
e
f
o
r
I
n
f
o
r
m
a
t
i
c
s
a
n
d
A
u
t
o
m
a
t
i
o
n
R
e
s
e
a
r
c
h
(
I
N
R
I
A
)
[
1
5
]
.
T
h
e
e
f
f
e
c
t
i
v
e
n
e
s
s
o
f
C
A
w
a
s
d
e
t
e
r
m
i
n
e
d
t
h
r
o
u
g
h
t
h
e
c
o
m
p
a
r
i
s
o
n
o
f
b
o
t
h
C
A
a
n
d
m
o
d
i
f
i
e
d
E
A
s
o
l
u
t
i
o
n
s
,
f
o
c
u
s
i
n
g
o
n
a
c
c
u
r
a
c
y
a
n
d
p
r
o
c
e
s
s
i
n
g
t
i
m
e
.
4.
R
E
S
U
L
T
S
A
N
D
A
N
A
L
Y
S
I
S
R
e
s
u
l
t
s
o
f
t
h
r
e
e
f
i
r
s
t
o
r
d
e
r
O
D
E
s
w
e
r
e
r
e
c
o
r
d
e
d
.
T
h
e
s
e
r
e
s
u
l
t
s
w
e
r
e
o
b
t
a
i
n
e
d
b
y
t
e
s
t
i
n
g
C
A
w
i
t
h
t
h
r
e
e
s
e
t
s
o
f
l
i
n
e
a
r
O
D
E
s
u
s
i
n
g
t
h
r
e
e
d
i
f
f
e
r
e
n
t
s
t
e
p
s
i
z
e
s
:
0
.
1
,
0
.
0
1
a
n
d
0
.
0
0
1
.
T
h
e
s
e
s
e
t
s
w
e
r
e
t
e
s
t
e
d
f
o
r
e
f
f
i
c
i
e
n
c
y
u
s
i
n
g
t
h
e
m
i
n
i
m
u
m
a
n
d
m
a
x
i
m
u
m
e
r
r
o
r
f
o
r
t
h
e
w
h
o
l
e
c
y
c
l
e
.
T
a
b
l
e
1
i
l
l
u
s
t
r
a
t
e
s
t
h
e
e
x
a
c
t
s
o
l
u
t
i
o
n
s
Evaluation Warning : The document was created with Spire.PDF for Python.
I
n
d
o
n
e
s
i
a
n
J
E
l
e
c
E
n
g
&
C
o
m
p
S
c
i
I
S
S
N
:
2502
-
4
7
5
2
T
C
u
b
e
A
r
i
t
h
m
e
t
i
c
:
I
m
p
r
o
v
i
n
g
E
u
l
e
r
M
e
t
h
o
d
f
o
r
O
r
d
i
n
a
r
y
D
i
f
f
e
r
e
n
t
i
a
l
E
q
u
a
t
i
o
n
…
(
N
o
o
r
a
i
d
a
S
a
m
s
u
d
i
n
)
1111
f
o
r
t
h
e
O
D
E
s
.
T
a
b
l
e
2
i
l
l
u
s
t
r
a
t
e
s
t
h
e
c
o
m
p
a
r
i
s
o
n
r
e
s
u
l
t
s
b
e
t
w
e
e
n
E
A
a
n
d
C
A
m
e
t
h
o
d
s
.
R
e
l
a
t
i
v
e
e
r
r
o
r
b
a
s
e
d
o
n
[
1
7
]
w
a
s
c
a
l
c
u
l
a
t
e
d
a
s
:
E
r
r
o
r
=
]
[
Ev
Ex
,
Ex
=
E
x
a
c
t
_
v
a
l
u
e
a
n
d
Ev
=
E
u
l
e
r
’
s
_
m
o
d
i
f
i
e
d
_
v
a
l
u
e
T
a
b
l
e
1
.
S
e
t
o
f
P
r
o
b
l
e
m
F
i
r
s
t
O
r
d
e
r
O
D
E
s
E
q
u
a
t
i
o
n
E
x
a
c
t
S
o
l
u
t
i
o
n
I
n
i
t
i
a
l
V
a
l
u
e
s
I
n
t
e
r
v
a
l
o
f
I
n
t
e
g
r
a
t
i
o
n
S
o
u
r
c
e
y
’=
-
0
.
5
y
e
(
-
0
.
5
x
)
y
(
0
)
=
0
0
≤
x
≤
2
0
[
1
6
]
y
’
=
-
y
e
(
-
x
)
y
(
0
)
=
1
0
≤
x
≤
2
0
[
1
6
]
y
’
=
-
10
y
e
(
-
10
x
)
y
(
0
)
=
1
0
≤
x
≤
2
0
[
1
6
]
T
a
b
l
e
2
.
R
e
s
u
l
t
s
o
f
E
r
r
o
r
s
U
s
i
n
g
V
a
r
i
o
u
s
S
t
e
p
S
i
z
e
h
M
e
t
h
o
d
E
u
l
e
r
A
r
i
t
h
m
e
t
i
c
C
u
b
e
A
r
i
t
h
m
e
t
i
c
S
t
e
p
S
i
z
e
h
=
0
.
0
0
1
h
=
0
.
0
1
h
=
0
.
1
h
=
0
.
0
0
1
h
=
0
.
0
1
h
=
0
.
1
P
r
o
b
l
e
m
1
0
0
0
.
0
0
0
0
8
3
0
0
.
0
0
0
0
0
1
0
.
0
0
0
0
7
6
P
r
o
b
l
e
m
2
0
0
.
0
0
0
0
0
3
0
.
0
0
0
3
5
5
0
0
.
0
0
0
0
0
3
0
.
0
0
3
0
2
P
r
o
b
l
e
m
3
0
0
.
0
0
0
4
0
1
0
.
1
1
7
9
0
0
0
0
.
0
0
0
3
0
0
0
.
0
2
9
0
2
0
T
h
e
r
e
s
u
l
t
s
f
o
r
e
a
c
h
m
e
t
h
o
d
w
e
r
e
c
o
m
p
a
r
e
d
t
o
t
h
e
e
x
a
c
t
s
o
l
u
t
i
o
n
u
s
i
n
g
m
a
x
i
m
u
m
e
r
r
o
r
.
I
t
i
s
e
v
i
d
e
n
t
f
r
o
m
T
a
b
l
e
2
t
h
a
t
t
h
e
p
r
o
p
o
s
e
d
C
A
m
e
t
h
o
d
p
r
o
v
i
d
e
d
a
m
o
r
e
a
c
c
u
r
a
t
e
r
e
s
u
l
t
t
h
a
n
E
A
a
t
h
i
g
h
e
r
h
s
i
z
e
s
.
S
m
a
l
l
e
r
s
t
e
p
s
i
z
e
p
r
o
d
u
c
e
d
s
i
m
i
l
a
r
r
e
s
u
l
t
s
f
o
r
b
o
t
h
m
e
t
h
o
d
s
.
P
r
o
b
l
e
m
1
d
e
m
o
n
s
t
r
a
t
e
d
t
h
a
t
t
h
e
r
e
s
u
l
t
s
o
f
t
h
e
C
A
m
e
t
h
o
d
i
s
b
e
t
t
e
r
t
h
a
n
E
A
a
t
h
=
0
.
1
w
i
t
h
a
m
a
x
i
m
u
m
e
r
r
o
r
o
f
0
.
0
0
0
0
7
6
a
n
d
0
.
0
0
0
0
8
3
r
e
s
p
e
c
t
i
v
e
l
y
.
F
i
g
u
r
e
1
i
l
l
u
s
t
r
a
t
e
s
t
h
e
d
i
f
f
e
r
e
n
c
e
s
b
e
t
w
e
e
n
t
h
e
t
w
o
m
e
t
h
o
d
s
.
R
e
s
u
l
t
s
a
t
h
=
0
.
0
1
s
h
o
w
e
d
t
h
a
t
C
A
s
c
o
r
e
d
0
.
0
0
0
0
0
1
c
o
m
p
a
r
e
d
t
o
E
A
a
t
0
.
B
o
t
h
m
e
t
h
o
d
s
s
c
o
r
e
d
t
h
e
s
a
m
e
a
t
h
=
0
.
0
0
1
.
P
r
o
b
l
e
m
2
d
e
m
o
n
s
t
r
a
t
e
d
t
h
a
t
t
h
e
r
e
s
u
l
t
s
f
o
r
b
o
t
h
C
A
a
n
d
E
A
m
e
t
h
o
d
s
a
r
e
a
l
m
o
s
t
t
h
e
s
a
m
e
a
t
h
=
0
.
1
,
0
.
0
1
a
n
d
0
.
0
0
1
,
a
s
s
h
o
w
n
i
n
F
i
g
u
r
e
2
.
F
i
g
u
r
e
1
.
D
i
f
f
e
r
e
n
c
e
s
b
e
t
w
e
e
n
t
h
e
t
w
o
m
e
t
h
o
d
s
F
i
g
u
r
e
1
.
T
h
e
r
e
s
u
l
t
s
f
o
r
b
o
t
h
C
A
a
n
d
E
A
m
e
t
h
o
d
s
Evaluation Warning : The document was created with Spire.PDF for Python.
I
S
S
N
:
2502
-
4
7
5
2
I
n
d
o
n
e
s
i
a
n
J
E
l
e
c
E
n
g
&
C
o
m
p
S
c
i
,
V
o
l
.
11
,
N
o
.
3
,
S
e
p
t
e
m
b
e
r
2
0
1
8
:
1109
–
1113
1112
P
r
o
b
l
e
m
3
d
e
m
o
n
s
t
r
a
t
e
d
t
h
a
t
t
h
e
C
A
m
e
t
h
o
d
s
c
o
r
e
d
b
e
t
t
e
r
r
e
s
u
l
t
s
c
o
m
p
a
r
e
d
t
o
E
A
a
t
e
a
c
h
s
t
e
p
s
i
z
e
.
A
t
h
=
0
.
1
,
C
A
s
c
o
r
e
d
0
.
0
2
9
0
2
0
w
h
i
l
s
t
E
A
s
c
o
r
e
d
0
.
1
1
7
9
0
0
,
a
n
d
a
t
h
=
0
.
0
1
,
C
A
s
c
o
r
e
d
0
.
0
0
0
3
0
0
w
h
i
l
s
t
E
A
s
c
o
r
e
d
0
.
0
0
0
4
0
,
a
s
s
h
o
w
n
i
n
F
i
g
u
r
e
3
.
F
i
g
u
r
e
2
.
T
h
e
r
e
s
u
l
t
s
f
o
r
b
o
t
h
C
A
a
n
d
E
A
m
e
t
h
o
d
s
I
t
c
a
n
b
e
s
u
m
m
a
r
i
z
e
d
t
h
a
t
t
h
e
p
r
o
p
o
s
e
d
C
A
m
e
t
h
o
d
p
r
o
v
i
d
e
s
a
m
o
r
e
a
c
c
u
r
a
t
e
r
e
s
u
l
t
t
h
a
n
E
A
i
n
s
o
l
v
i
n
g
O
D
E
s
a
t
h
i
g
h
e
r
s
t
e
p
s
i
z
e
.
5.
C
O
N
C
L
U
S
I
O
N
T
h
i
s
p
a
p
e
r
h
a
s
p
r
o
p
o
s
e
d
C
u
b
e
A
r
i
t
h
m
e
t
i
c
,
a
n
e
w
m
e
t
h
o
d
u
s
i
n
g
m
o
d
i
f
i
e
d
E
u
l
e
r
a
s
f
i
n
d
i
n
g
o
f
t
h
i
s
s
t
u
d
y
.
T
h
e
f
e
a
s
i
b
i
l
i
t
y
o
f
t
h
e
p
r
o
p
o
s
e
d
a
l
g
o
r
i
t
h
m
w
a
s
c
o
m
p
a
r
e
d
w
i
t
h
E
u
l
e
r
A
r
i
t
h
m
e
t
i
c
a
n
d
t
e
s
t
e
d
u
s
i
n
g
S
C
I
L
A
B
.
S
u
b
s
e
q
u
e
n
t
l
y
,
t
h
e
C
u
b
e
A
r
i
t
h
m
e
t
i
c
c
o
m
p
a
r
e
d
w
i
t
h
E
u
l
e
r
A
r
i
t
h
m
e
t
i
c
s
c
h
e
m
e
w
i
t
h
t
h
e
e
x
a
c
t
s
o
l
u
t
i
o
n
.
U
s
u
a
l
l
y
,
t
h
e
o
r
d
i
n
a
r
y
E
u
l
e
r
m
e
t
h
o
d
u
s
i
n
g
a
s
m
a
l
l
s
t
e
p
s
i
z
e
g
i
v
e
s
t
h
e
s
o
l
u
t
i
o
n
a
l
m
o
s
t
t
o
t
h
e
e
x
a
c
t
s
o
l
u
t
i
o
n
.
H
o
w
e
v
e
r
,
i
n
t
h
i
s
s
t
u
d
y
h
a
s
b
e
e
n
d
e
m
o
n
s
t
r
a
t
e
d
t
h
a
t
C
u
b
e
A
r
i
t
h
m
e
t
i
c
p
r
o
v
i
d
e
d
s
o
l
u
t
i
o
n
s
t
h
a
t
a
r
e
s
i
m
i
l
a
r
t
o
e
x
a
c
t
s
o
l
u
t
i
o
n
s
a
t
s
m
a
l
l
s
t
e
p
s
i
z
e
a
n
d
a
l
s
o
a
t
h
i
g
h
e
r
s
t
e
p
s
i
z
e
.
T
h
e
b
e
n
e
f
i
t
s
o
f
u
s
i
n
g
s
m
a
l
l
s
t
e
p
s
i
z
e
wi
l
l
g
i
v
e
a
h
i
g
h
e
r
a
c
c
u
r
a
c
y
.
T
h
e
r
e
f
o
r
e
a
h
i
g
h
e
r
s
t
e
p
s
i
z
e
w
o
u
l
d
r
e
d
u
c
e
c
o
m
p
l
e
x
i
t
y
a
n
d
p
r
o
c
e
s
s
i
n
g
t
i
m
e
.
A
s
a
c
o
n
c
l
u
s
i
o
n
,
C
u
b
e
A
r
i
t
h
m
e
t
i
c
c
a
n
b
e
u
s
e
d
a
s
a
n
a
l
t
e
r
n
a
t
i
v
e
a
l
g
o
r
i
t
h
m
t
o
s
o
l
v
e
O
D
E
p
r
o
b
l
e
m
s
.
A
C
K
N
O
WL
E
D
G
M
E
N
T
S
T
h
e
a
u
t
h
o
r
s
g
r
a
t
e
f
u
l
l
y
a
c
k
n
o
w
l
e
d
g
e
t
h
e
F
a
c
u
l
t
y
o
f
S
c
i
e
n
c
e
a
n
d
T
e
c
h
n
o
l
o
g
y
D
e
f
e
n
s
e
,
U
n
i
v
e
r
s
i
t
i
P
e
r
t
a
h
a
n
a
n
N
a
s
i
o
n
a
l
M
a
l
a
y
s
i
a
.
T
h
e
a
u
t
h
o
r
N
o
o
r
a
i
d
a
S
a
m
s
u
d
i
n
a
l
s
o
w
o
u
l
d
l
i
k
e
t
o
t
h
a
n
k
T
A
T
I
U
n
i
v
e
r
s
i
t
y
C
o
l
l
e
g
e
f
o
r
t
h
e
S
h
o
r
t
T
e
r
m
G
r
a
n
t
(
S
T
G
1
/
2
0
1
8
)
.
R
E
F
E
R
E
N
C
E
S
[1]
A
M
odi
fi
e
d
O
D
E
S
ol
ve
r
for
A
ut
onom
ous
Ini
t
i
a
l
V
a
l
ue
P
robl
e
m
s
Z
a
i
b
-
un
-
N
i
s
a
M
e
m
on,
S
a
ni
a
Q
ure
s
hi
,
A
s
i
f
A
l
i
S
ha
i
kh,
M
uha
m
m
a
d
S
a
l
e
e
m
Cha
ndi
o,
M
a
t
he
m
a
t
i
c
a
l
T
he
ory
a
nd
M
ode
l
i
ng
w
w
w
.i
i
s
t
e
.org
IS
S
N
2224
-
5804
(P
a
p
e
r)
IS
S
N
2225
-
0522 (O
nl
i
ne
)
V
ol
.4,
N
o.3, 2014
[2]
Z
.
S
a
l
l
e
h,
O
rdi
na
ry
D
i
ffe
re
nt
i
a
l
E
qua
t
i
ons
(O
D
E
)
U
s
i
ng
E
ul
e
r'
s
T
e
c
hni
que
a
nd
S
c
i
l
a
b
P
rog
ra
m
m
i
ng
,
M
a
t
he
m
a
t
i
c
a
l
M
ode
l
s
a
nd M
e
t
hods
i
n M
ode
rn S
c
i
e
nc
e
20(4) (2012), p
p
. 264
-
269
[3]
A.
J
a
bba
r
,
F
.M
.
M
a
l
i
k
,
S
a
m
p
l
e
d
-
D
a
t
a
Ba
c
k
s
t
e
p
p
i
ng
Cont
rol
of
a
Q
ua
drot
or
U
nm
a
nne
d
A
e
ri
a
l
V
e
hi
c
l
e
.
IA
E
S
Int
e
rna
t
i
ona
l
J
ourna
l
of
Robot
i
c
s
a
nd
A
ut
om
a
t
i
on
(IJ
RA
)
,
V
ol
4 N
o 2, 2015;
p
p
.
124
-
134.
[4]
Y
.L
i
,
Bre
a
ki
ng
t
he
D
i
g
i
t
a
l
V
i
de
o
S
t
e
g
a
nog
ra
p
hy
.
T
E
L
K
O
M
N
IK
A
Indone
s
i
a
n
J
ourna
l
of
E
l
e
c
t
ri
c
a
l
E
ng
i
ne
e
ri
ng
,
V
ol
11
No
3,
2013,
PP
1691
-
1696
[5]
P
. H
e
nri
c
i
,
D
i
s
c
re
t
e
V
a
ri
a
bl
e
M
e
t
hods
i
n O
rdi
na
ry
D
i
ffe
re
nt
i
a
l
E
qua
t
i
on, J
ohn W
i
l
e
y
&
S
ons
, 1962, p
p
. 187
-
190.
[6]
E
. K
re
y
s
z
i
g
, A
dva
nc
e
d E
ng
i
ne
e
ri
ng
M
a
t
he
m
a
t
i
c
s
, 9t
h e
d., J
ohn W
i
l
e
y
a
nd S
ons
, 2006, p
p
. 19
-
20
[7]
S
.
Q
ure
s
hi
,
M
.S
.
Cha
ndi
o,
I.A
.
J
une
j
o,
A
.A
.
S
ha
i
kh
a
nd
Z
.
-
U
.N
.
M
e
m
on,
O
n
E
rror
Bound
for
L
oc
a
l
T
runc
a
t
i
on
E
rror
of
a
n
E
x
p
l
i
c
i
t
It
e
ra
t
i
ve
A
l
g
ori
t
hm
i
n
O
rdi
na
ry
D
i
ffe
re
nt
i
a
l
E
qua
t
i
ons
,S
c
i
.Int
.(L
a
hore
)26(3)
(2014),
p
p
.
1413
-
1415.
[8]
N
.M
.M
.
Y
us
op
a
nd
M
.K
.
H
a
s
a
n
a
nd
M
.
Ra
hm
a
t
,
C
om
p
a
ri
s
on
N
e
w
A
l
g
ori
t
hm
M
odi
fi
e
d
E
ul
e
r
i
n
O
rdi
na
ry
D
i
ffe
re
nt
i
a
l
E
qua
t
i
on U
s
i
ng
S
c
i
l
a
b
P
rog
ra
m
m
i
ng
.
L
e
c
t
ure
N
ot
e
s
on S
oft
w
a
re
E
ng
i
ne
e
ri
ng
,
2015. 3(3):
p
p
. 199
-
202
Evaluation Warning : The document was created with Spire.PDF for Python.
I
n
d
o
n
e
s
i
a
n
J
E
l
e
c
E
n
g
&
C
o
m
p
S
c
i
I
S
S
N
:
2502
-
4
7
5
2
T
C
u
b
e
A
r
i
t
h
m
e
t
i
c
:
I
m
p
r
o
v
i
n
g
E
u
l
e
r
M
e
t
h
o
d
f
o
r
O
r
d
i
n
a
r
y
D
i
f
f
e
r
e
n
t
i
a
l
E
q
u
a
t
i
o
n
…
(
N
o
o
r
a
i
d
a
S
a
m
s
u
d
i
n
)
1113
[9]
S
.
F
a
dug
ba
,
B.
O
g
unri
nde
a
nd
T
.
O
kunl
ol
a
,
E
ul
e
r'
s
M
e
t
hod
for
S
ol
vi
ng
Ini
t
i
a
l
V
a
l
ue
P
robl
e
m
s
I
nO
rdi
na
ry
D
i
ffe
re
nt
i
a
l
E
qua
t
i
ons
, T
he
P
a
c
i
fi
c
J
ourna
l
of S
c
i
e
nc
e
T
e
c
hnol
og
y
13(2)(2012), p
p
.152
-
158.
[10]
N.M.M.
Y
us
op,
M.K.
H
as
an,
M.
Wook
,
M.F
.M.
A
m
r
an,
S.R
.
A
hm
ad,
Com
p
a
ri
s
on
N
e
w
A
l
g
ori
t
hm
M
odi
fi
e
d
E
ul
e
r
Ba
s
e
d
O
n
H
a
rm
oni
c
-
P
ol
y
g
on
A
p
p
roa
c
h
F
or
S
ol
vi
ng
O
rdi
na
ry
D
i
ffe
re
nt
i
a
l
E
qua
t
i
on,
J
ourna
l
of T
e
l
e
c
om
m
uni
c
a
t
i
on,
E
l
e
c
t
roni
c
a
nd
Com
p
ut
e
r E
ng
i
ne
e
ri
ng
,
V
ol
um
e
9, Is
s
ue
2
-
11, 2017, p
p
. 29
-
32
[11]
N.M.M.
Y
us
op,
M.K.
H
as
an,
M.
Wook
,
M.F
.M.
A
m
r
an,
S.R
.
A
hm
ad
,
A
ne
w
E
ul
e
r
s
c
he
m
e
ba
s
e
d
on
ha
rm
oni
c
-
p
ol
y
g
on
a
p
p
roa
c
h
for
s
ol
vi
ng
fi
rs
t
orde
r
ordi
na
ry
di
ffe
re
nt
i
a
l
e
qua
t
i
on.
T
he
2nd
Int
e
rna
t
i
ona
l
Confe
re
nc
e
O
n
A
p
p
l
i
e
d S
c
i
e
nc
e
A
nd T
e
c
hnol
og
y
2017 (ICA
S
T
’17)
[12]
F
.
A
g
us
,
H
a
vi
l
uddi
n
,
S
c
i
l
a
b
s
oft
w
a
re
a
s
a
n
a
l
t
e
rna
t
i
ve
l
ow
-
c
os
t
c
om
p
ut
i
ng
i
n
s
ol
vi
ng
t
he
l
i
ne
a
r
e
qua
t
i
ons
p
robl
e
m
.
A
IP
Confe
re
nc
e
P
roc
e
e
di
ng
s
(2017)
[13]
N
.M
.M
.
Y
us
op
a
nd
M
.K
.
H
a
s
a
n,
D
e
ve
l
op
m
e
nt
of
N
e
w
H
a
rm
oni
c
E
ul
e
r
U
s
i
ng
N
ons
t
a
nda
rd
F
i
ni
t
e
D
i
ffe
re
nc
e
T
e
c
hni
que
for
S
ol
vi
ng
S
t
i
ff P
robl
e
m
s
,
J
urna
l
T
e
knol
og
i
(S
c
i
e
nc
e
s
a
nd
E
ng
i
ne
e
ri
ng
),
7(22) (2014), p
p
. 19
-
24
[14]
M
.A
.Is
l
a
m
,
A
c
c
ura
c
y
A
na
l
y
s
i
s
of
N
um
e
ri
c
a
l
s
ol
ut
i
ons
of
i
ni
t
i
a
l
va
l
ue
p
robl
e
m
s
(IV
P
)
for
ordi
na
ry
di
ffe
re
nt
i
a
l
e
qua
t
i
ons
(O
D
E
)
.
IO
S
R J
ourna
l
of M
a
t
he
m
a
t
i
c
s
. V
ol
um
e
11, Is
s
ue
3 V
e
r.
III (2015), P
P
18
-
23
[15]
IN
RIA
-
E
N
P
C (2017). S
c
i
l
a
b
(V
e
rs
i
on
6.0). U
RL
ht
t
p
:
/
/
w
w
w
.s
c
i
l
a
b.org
/
[16]
Z
a
ri
na
,
B.I.,
M
oha
m
e
d,
S
.,
K
ha
i
ri
l
,
I.
a
nd
Z
a
na
ri
a
h,
M
.,
Bl
oc
k
M
e
t
hod
for
G
e
ne
ra
l
i
s
e
d
M
ul
t
i
s
t
e
p
A
da
m
s
a
nd
Ba
c
kw
a
rds
D
i
ffe
re
nt
i
a
t
i
on
F
orm
ul
a
e
In S
ol
vi
ng
F
i
rs
t
O
rde
r O
D
E
s
. M
A
T
E
M
A
T
IK
A
, 2005. 21(1):
p
. 25
-
11
[17]
R. J
a
i
s
w
a
l
, A
. A
. P
a
t
ha
n,
S
t
udy
of N
um
e
ri
c
a
l
A
na
l
y
s
i
s
–
D
i
ffe
re
nt
i
a
l
E
qua
t
i
on
Int
e
rna
t
i
ona
l
J
ourna
l
of A
dva
nc
e
d
Re
s
e
a
rc
h
i
n Com
p
ut
e
r
S
c
i
e
nc
e
a
nd S
oft
w
a
re
E
ng
i
ne
e
ri
ng
,
V
ol
um
e
5, Is
s
ue
10, (2015)
Evaluation Warning : The document was created with Spire.PDF for Python.