#
ββThe symbols which do not occur in Algorithm 1 , are explained below.
#
π€
π
βvariable symbol in integral for allowing certain transformations.
#
π’
βabscissa of point
πΆ
π€
in the
π€
-plane.
#
ββassume
(
0
<
π
<
1
<
π
<
β
)
;
#
π»
π
π₯
βtangential component of the field strength.
#
π»
π
π¦
βnormal component of the field strength.
#
ββWarning. In order to keep the usual notation for the components of the field strength we have used the indexed names
π»
π₯
and
π»
π¦
on the figures, but in the program we are obliged to modify the indices by putting, for example,
π»
π
π₯
and
π»
π
π¦
as components along the axes, otherwise, if in the program the ordinate
π¦
also exists, errors appear, and the obtained values have no more meaning. For avoiding this circumstance, it is also possible to use only symbols.
π
βΆ
=
0
.
2
5
;
π
βΆ
=
4
.;
π
0
β
βΆ
=
(
π
β
1
)
/
π
;
π»
π¦
0
βΆ
=
1
;
πΉ
βΆ
=
π©
π«
π¨
π
(
π
)
g
e
n
e
r
a
l
π§
,
π
;
#
ββThe transformation function is:
ξ
π
βΆ
=
(
1
/
π
)
β
π€
π
2
β
(
π
2
0
+
2
.
)
π€
π
+
1
.
/
(
π€
π
β
(
π€
π
β
1
.
)
)
;
π§
βΆ
=
i
n
t
(
π
,
π€
π
)
;
π§
βΆ
=
e
v
a
l
(
π§
,
π€
π
=
π
)
+
0
.
5
πΌ
;
#
ββThe last imaginary constant has been introduced for to be in accordance with the configuration of Figure 1 .end proc;
π’
βΆ
=
πΉ
(
β
1
)
;
π’
βΆ
=
e
v
a
l
(
π’
)
;
#
ββCo-ordinate introduced for to have the value zero of the
π₯
-abscissa at the point of the minimum of
π»
π
π¦
.forββ
π€
ββfromββ β0.01ββbyββ β0.25ββtoββ β1.25ββdo
p
r
i
n
(
π€
)
;
π§
βΆ
=
πΉ
(
π€
)
;
π₯
βΆ
=
β
(
π§
β
π’
)
;
π¦
βΆ
=
π
(
π§
)
;
π»
π
π₯
ξ
βΆ
=
β
(
c
o
n
j
u
g
a
t
e
(
(
(
π€
β
1
.
)
/
π€
2
β
(
π
2
0
+
2
.
)
π€
+
1
.
)
πΌ
)
)
;
π»
π
π¦
ξ
βΆ
=
β
β
(
c
o
n
j
u
g
a
t
e
(
(
(
π€
β
1
.
)
/
π€
2
β
(
π
2
0
+
2
.
)
π€
+
1
.
)
πΌ
)
)
;
#
ββThe sign minus has been introduced because, according to the assumed conditions in
π‘
-plane (zeta plane), the field strength is up-down.end do ;
#
p
l
o
t
(
[
π₯
,
π»
π
π¦
,
π€
=
β
1
.
.
β
1
0
0
]
)
;
#
p
l
o
t
(
[
π₯
,
π»
π
π¦
,
π€
=
β
1
0
0
.
.
β
1
0
0
0
]
)
;
#
p
l
o
t
(
[
π₯
,
π»
π
π¦
,
π€
=
β
1
0
0
0
.
.
β
1
0
0
0
0
0
]
)
;
#
p
l
o
t
(
[
π₯
,
π»
π
π¦
,
π€
=
β
1
.
.
β
0
.
0
1
]
)
;
#
p
l
o
t
(
[
π₯
,
π»
π
π¦
,
π€
=
β
0
.
0
1
.
.
β
0
.
0
0
1
]
)
;
#
p
l
o
t
(
[
π₯
,
π»
π
π¦
,
π€
=
β
0
.
0
0
1
.
.
β
0
.
0
0
0
0
1
]
)
;
#
ββFor plotting by using the same program, it suffices to disable the lines for w from β¦tillββdo by putting
#
and enabling the plot statements by removing
#
.
#
ββResults after running the program:ββββββββββββββ
π
βΆ
=
0
.
2
5
ββββββββββββββ
π
βΆ
=
4
. ββββββββββββββ
π
0
βΆ
=
1
.
5
0
0
0
0
0
0
0
0
ββββββββββββββ
π»
π¦
0
βΆ
=
1
Warning, βββfβ ββis ββimplicitly ββdeclared ββlocal ββto ββprocedure βββFβ Warning, βββzβ ββis ββimplicitly ββdeclared ββlocal ββto ββprocedure βββFβ
πΉ
βΆ
=
π©
π«
π¨
π
(
π
)
βlocal
π
,
π§
;
βgeneral*
π§
,
π
; β
π
βΆ
=
s
q
r
t
(
π€
π
2
β
(
π
2
0
+
2
.
)
β
π€
π
+
1
.
)
/
(
π
β
π€
π
β
(
π€
π
β
1
.
)
)
; β
π§
βΆ
=
i
n
t
(
π
,
π€
π
)
; β
π§
βΆ
=
e
v
a
l
(
π§
,
π€
π
=
π
)
+
0
.
5
β
πΌ
end proc ;ββββββββ
π’
βΆ
=
0
.
4
2
0
7
2
7
9
5
1
1
+
1
.
0
0
0
0
0
0
0
0
0
πΌ
ββββββββ
π’
βΆ
=
0
.
4
2
0
7
2
7
9
5
1
1
+
1
.
0
0
0
0
0
0
0
0
0
πΌ
βββββββββββββ
p
r
i
n
(
β
0
.
0
1
)
ββββββββ
π§
βΆ
=
2
.
0
4
8
2
5
5
5
8
9
+
1
.
0
0
0
0
0
0
0
0
0
πΌ
ββββββββββββ
π₯
βΆ
=
1
.
6
2
7
5
2
7
6
3
8
βββββββββββ
π¦
βΆ
=
1
.
0
0
0
0
0
0
0
0
0
βββββββββββββ
π»
π
π₯
βΆ
=
0
. ββββββββββ
π»
π
π¦
βΆ
=
β
0
.
9
8
9
1
5
0
8
1
5
6
ββββββββββββ
p
r
i
n
(
β
0
.
2
6
)
ββββββββ
π§
βΆ
=
0
.
9
4
4
1
0
6
3
1
0
5
+
1
.
0
0
0
0
0
0
0
0
0
πΌ
βββββββββββ
π₯
βΆ
=
0
.
5
2
3
3
7
8
3
5
9
4
βββββββββββ
π¦
βΆ
=
1
.
0
0
0
0
0
0
0
0
0
ββββββββββββββ
π»
π
π₯
βΆ
=
0
. βββββββββββ
π»
π
π¦
βΆ
=
β
0
.
8
5
4
8
3
1
7
7
5
8
βββββββββββββ
p
r
i
n
(
β
0
.
5
1
)
ββββββββ
π§
βΆ
=
0
.
6
8
6
8
8
9
9
8
2
9
+
1
.
0
0
0
0
0
0
0
0
0
πΌ
βββββββββββ
π₯
βΆ
=
0
.
2
6
6
1
6
2
0
3
1
8
βββββββββββ
π¦
βΆ
=
1
.
0
0
0
0
0
0
0
0
0
βββββββββββββ
π»
π
π₯
βΆ
=
0
. ββββββββββ
π»
π
π¦
βΆ
=
β
0
.
8
1
5
6
0
8
7
5
7
4
βββββββββββββ
p
r
i
n
(
β
0
.
7
6
)
ββββββββ
π§
βΆ
=
0
.
5
2
9
8
0
0
3
3
9
7
+
1
.
0
0
0
0
0
0
0
0
0
πΌ
βββββββββββ
π₯
βΆ
=
0
.
1
0
9
0
7
2
3
8
8
6
βββββββββββ
π¦
βΆ
=
1
.
0
0
0
0
0
0
0
0
0
βββββββββββββ
π»
π
π₯
βΆ
=
0
. ββββββββββ
π»
π
π¦
βΆ
=
β
0
.
8
0
2
6
9
1
2
0
5
1
ββββββββββββ
p
r
i
n
(
β
1
.
0
1
)
ββββββββ
π§
βΆ
=
0
.
4
1
6
7
6
8
8
4
5
6
+
1
.
0
0
0
0
0
0
0
0
0
πΌ
ββββββββββ
π₯
βΆ
=
β
0
.
0
0
3
9
5
9
1
0
5
5
βββββββββββ
π¦
βΆ
=
1
.
0
0
0
0
0
0
0
0
0
βββββββββββββ
π»
π
π₯
βΆ
=
0
. ββββββββββ
π»
π
π¦
βΆ
=
β
0
.
8
0
0
0
0
3
5
6
4
3
ββββββββββ( 1)