Potential and Field 
in a Plane Condensor with a shifted Voltage Strip 
in comparison with
a shifted Semi-Infinite Plane Condensor

B. Schnizer, Th. Heubrandtner, St. Rossegger

Institut fuer Theoretische Physik
Technische Universität Graz
A-8010 Graz, Austria, Europe
Email: schnizer@itp.tu-graz.ac.at

1.  Description

Plane Condensor with a shifted Voltage Strip

The potential and field distributions for the configuration shown below, a twodimensional infinite plane condensor 
of width  D, whose upper electrode has a strip of width  w,  are computed by conformal maps. The strip is at  
a Voltage  V= cvV and can be shifted along the x-axis with the value c;  the other electrodes are grounded.

The main parameters are the Voltage  cvV  applied to the Voltage strip, the width of the condensor   D = cvgd,
the width of the strip  w  =  cvw  and the length of shifting c=cvc. Default values are assigned to these as to all other variables. 

More Information about the Conformal maps and the formulas to derive the field ar given in the notebook
CondShCoStr.nb.

shifted Semi Infinite Plane Condensor

The potential and field distributions for the configuration shown below, a twodimensional semi-infinite plane condensor 
of width  D=huD, whose upper electrode extends from minus infinity to b=hub,  are computed by conformal maps. The upper electrode is at  Voltage V=huV ;  the lower electrode, which is at the x-axis, is grounded.

The main parameters are the Voltage ϕ=huV  applied to the upper electrode, the width of the condensor   D = scgd  and  
the shifting length  b = hub. Default values are assigned to these as to all other variables.

More Information about the Conformal maps and the formulas to derive the field ar given in the notebook
ShSemInf.nb.

Comparison of the two configurations described above

The plots of the potential and the field are overlapped for configurations.

The potential
and the field are plotted,as is the field on the electrodes.The values of the potential and field are computed at a point
which may be selected by the user.
The relative current  I/qv  induced on the upper left electrode by a charge q moving with constant velocity  v  
on various trajectories parallel to the y-axis is shown as function of y.  The user may also chose the x-coordinate 
of such a trajectory.

The cells where values are assigned to the important parameters have a 
background in red (they will be allocated to the specific variables of the two configurations - Background cyan) so that the user may locate them easily; he may assign other values to these variables and evaluate the notebook 
again; by tests we are quite sure that the program will run without difficulty. 
The user may play with some other parameters so that he gets more satisfactory pictures. Such recommendations are also
given in cells with cyan background.
Please, report errors and difficulties per email to the adress given above in the affiliation.

2.  Input

Needs["Graphics`ComplexMap`"]

Off[General::spell1]

In this notebook many variables have similar names. The ensuing warnings are switched off above in the second line.
When working on the program these  warning may be useful; switch them on by the corresponding command On[...].
This is also implemented in the last line of this notebook.

Global input Data

\(\(\(Clear["cv*"]

    \(\(Clear["hu*"]\)\(
\)
    \)

    cvhuV = \ 1; \ \ \ \ \ \  (*\ 
      Potential\ on\ the\ upper\ elecrode\ in\ Volt, \ \ 
\ \ \ \ \ \ \ \ \ \ \ elsewhere\ on\ all\ electrodes\ Potential\  = \ 
        0\ *) 
cvhugd = 
      10; \ \ \ \ \  (*\ \(=\)\(\ \)\(D\  = \ 
          width\ of\ condensor\ in\ mm\)\ *) 
cvhub = \(-10\)\ ; \ \ \  (*\ \(=\)\(\ \)\(b\  = \ 
          position\ of\ the\ right\ end\ of\ the\ upper\ electrode\)\ *) ;

    \(\(cvhuw = \ 
      30\)\(\ \)\); \ \ \ \  (*\ \(=\)\(\ \)\(w\  = \ \ width\ of\ strip\ on\ upper\ electrode\)\ *) 
cvhutcb =  .011; 	 (*\ 
      Relative\ thickness\ of\ lines\ representing\ electrodes\ *) 
HuCol1 = Hue[0.3]; \ \ \ \ \ \  (*\ 
      Color\ of\ the\ ElectricField - Plots\  - \ 
        StripCondensor*) 
\(\(HuCol2 = 
      Hue[0.7]\)\(\ \)\);\)\(\ \ \ \ \ \)\( (*\ 
      Color\ of\ the\ ElectricField - Plots\  - \ 
        SplitCondensor*) \)\)\)

These Input Data are given as inputs to the two diffrent arangements

2.1  Input Data - Strip

cvV   =  cvhuV;        
cvgd = 
        cvhugd;      
cvw     =  
        cvhuw ;    
\(\(\(cvc\  = \ 
      cvhub\)\(\ \)\(;\)\( (*\ \(=\)\(\ \)\(c\  = \ 
          position\ of\ the\ right\ end\ of\ the\ strip\ on\ the\ x - 
            axis\)\ *) \)\(
\)\(cvtcb = 
      cvhutcb\)\(;\)\(	\)\( (*\ 
      Relative\ thickness\ of\ lines\ representing\ electrodes\ *) \)\(;\)\(
\)\)\)

\(Print[
      SequenceForm[
"  Width of condensor:                        D = ", \ cvgd, \ " mm."]]\)
\(Print[
      SequenceForm[
"  Width of strip at upper electrode :        w = ", \ cvw, \ " mm."]]\)
Print[
      SequenceForm[
"  Position of the right end of the \
strip:    c = ", cvc ,  " mm."]]
cvmc = cvc - cvw/2;
Print[
      SequenceForm["  Middle Postition of the strip:             mc = ", 
        cvmc ,  " mm."]]

2.2 Functions for the conformal map

2.2.1   Finding the value of the parameter p belonging to the dimensions selected in 2.1

The conformal map needs one parameter  α = cva = (e^(((π w)/(4 D)))),whose value depends on the ratio of the two geometric
parameters D and w.

cva  =  Exp[(π cvw)/(4 cvgd)];

2.2.2  The function  z  =  f(W)

The three conformal maps (Subscript[f, 1]), (Subscript[f, 2]), (Subscript[f, 3]) and the resulting map f  as described above are implemented below. 
Notations are:  W = cvWw,  t = cvtt, t' = cvtp,  α  = cva ; z = f(W) = cvf[cvWw].

cvf1[cvWw_]  =  Exp[cvWw] ;
cvf2[
          cvtt_]  =  - (cvtt cva^2   -  
                1)/(cvtt   -   cva^2) ;
cvf3[cvtp_]  =  
        cvgd/π  Log[cvtp]  + cvmc;
cvf[cvWw_]      =   cvf3[cvf2[cvf1[cvWw]]];

2.2.3 The inverse function  W = (f^(-1))(z)

The inverse functions to  (Subscript[f, 1])to (Subscript[f, 3]) may be expressed as analytic functions;  W = (f^(-1))(z)= cvi[sczz] :

cvi1[cvtt_] =  Log[cvtt] ;
cvi2[
          cvtp_] =  (cvtp  cva^2   +  
              1)/(cvtp  +   cva^2)  ;
cvi3[
          cvzz_] =   Exp[(cvzz  - cvmc) π/
              cvgd];
cvi[cvzz_]    =  cvi1[cvi2[cvi3[cvzz]]];

2.2.4  The complex electric field

The complex electric field   (Subscript["E", "x"])"\:f39f+ i "(Subscript["E", "y"])  at a point   z = x + i y   is given by the function  cvfe, which depends on 
W =cvWw, or by the function cvfez, which depends on cvzz. The value of W belonging to a given value  z  is 
found by  cvi. 
Note: Here MKSA units are used; but D = cvgd is given in mm, so the electric field computed from that formula is in  
kV/m.  (Subscript[E, z])(z) = cvfez[vszz] , (Subscript[E, z])(W) = cvfe[vsWw] .

cvfez[
          cvzz_] =   -I cvV  / cvgd/
              2 * (Tanh[π (Conjugate [cvzz - cvmc] + 
                          cvw/2)/(2 cvgd)]  + 
Tanh[π ( 
                        cvw/2  -  
                          Conjugate [cvzz - cvmc])/(2 cvgd)]) // 
          Chop ;

cvfe[
          cvWw_]  =   I cvV/
            cvgd (2 Cosh[Conjugate[cvWw]]  -  (cva^2  +  
                    1/cva^2))/(cva^2  -  1/cva^2);

2.3  Input Data

huV   = cvhuV;     ;
hugd = 
      cvhugd;       ;
hub = 
      cvhub;           ;
\(hutcb = cvhutcb; 	 (*\ 
      Relative\ thickness\ of\ lines\ representing\ electrodes\ *) ;\)

\(Print[
      SequenceForm[
"  Width of condensor:                                  D = ", \ hugd, \ " mm."]]\)
\(Print[
      SequenceForm[
"  position of the right end of the upper electrode :   b = ", \ \ hub, \ " mm."]]\)

2.4 Functions for the conformal map

2.2.1  The function  z  =  f(W)

The three conformal maps (Subscript[f, 1]), (Subscript[f, 2]), (Subscript[f, 3]), the length of shifting and the resulting map f  as described above are implemented below. 
Notations are:  W = huWw,  t = hutt, t' = hutp.

huf1[huWw_] = Exp[huWw];
huf2[hutt_] = hutt;
huf3[hutp_] = hugd/π (1 + Log[hutp] + hutp);
huf[scWw_]  =  huf3[huf2[huf1[scWw]]] + hub;

2.2.2 The inverse function  W = (f^(-1))(z)

It is not possible to express the inverse function (Subscript[f, 3]^(-1))(z)  as a closed analytic expression. So it must be computed 
numerically; this is done by a Newton iteration implemented in the command  FindRoot. This in turn needs an 
appropriate starting value depending on  z = huzz. Mathematica permits one to implement 
this procedure in such a way that the inverse function  W=(f^(-1))(z) = hui may be specified as a function of z.

hui[huzz_] :=  
        hupqq   /.  
          FindRoot[(huf[hupqq] ) - huzz, {hupqq, 
              1.2 Sign[Im[huzz] - hugd] + I π/2}]

2.2.3  The complex electric field

The complex electric field   (Subscript["E", "x"])"\:f39f+ i "(Subscript["E", "y"])  at a point   z = x + i y   is given by the function  hufe, which depends on 
W = huWw. The value of W belonging to the given value  z  is found by  hui.
Note: Here MKSA units are used; but D = hugd is given in 
mm, so the electric field given by that formula is in  kV/m.

hufe[
            huWw_] = -I (huV /hugd)*1/Conjugate[1 + Exp[huWw]];

3.  Section with Field Graphics

Most evaluations for finding values or graphs are done by using the values of W=U+i V belonging to a given z.For
curves appropriate ranges of U must be found.The range of V is almost always 0≤V≤π.

cvhuneq = 
        20  ;

3.1  Variables needed in the pictures - Strip

cvneq = 
        cvhuneq  ;

cvneq  specifies the number of equipotential and field lines;  this value in turn specifies the step widths used for the
variables  U ( = cvdu) and  V = π/cvneq used in cvcm. The range of the variable U is bounded by cvlu  and 
by cvru. These values have been determined by trial and error. The horizontal plot range in all pictures is 
determined by cvlx and cvrx. In some figures the vertical plot range is determined by cvmaE and multiples of it.

cvmaE  =  cvV/cvgd ;
cvplr  =  (cvgd  +  
            Abs[2 cvc]);
cvims  =  500;
cvlu  =  -cvi[cvplr];
cvru  =   cvi[cvplr];
cvdu  =  cvru/5;

The user may want to improve the pictures drawn below. For this, he may try to obtain better values for 
cvlu, cvru, cvdu, cvmaE. If he wants to do so it is recommended to remove the semicolons of the last block
at the end of the preceeding subsection, so that the values are printed. Then he may try a slight change of 
these values and see what comes out. It is commendable to insert these new values below in the
commands for  cvp2, cvp1, cvp23, cvcm. It is extremely commendable to save this  trial notebook
under a new name beforehand so that the notebook with the orginal specifications remains intact!

cvp2  computes the electric field at the inside of the upper electrode (cvp1 of the lower) .The corresponding pictureis shown in figure cvp23.  cvcm plots the equipotential and electric field lines.  cvp1 shows  the electric field at the bottom electrode. 
In the two cvp pictures the unit of the electric field strength is  kV/m  since in the definition of the functions  cvfe 
and cvfez in 2.2.4  D = cvgd is used, whose value is defined in mm in  1.1, while we use  MKSA units.

cvtru  =  { Hue[0], AbsoluteThickness[2], 
          Line[{{-cvw/2 + cvmc, -cvmaE}, {-cvw/2 + cvmc, 
                cvmaE}}], 
Line[{{cvw/2  + cvmc, - 
                  cvmaE}, {cvw/2 + cvmc, 
                cvmaE}}]};
cvtrl  =  { Hue[0], 
          Line[{{-cvw/2 + cvmc, 
                0}, {-cvw/2 + 
                  cvmc, -cvmaE}}], 
Line[{{cvw/2 + 
                  cvmc, 0}, {cvw/2 + cvmc, -cvmaE}}]};

cvp2 = Plot[{cvfez[xx + I cvgd]//Im},{xx, -cvplr, cvplr},DisplayFunction ->  Identity,PlotStyle -> {HuCol1}];
cvp1 = Plot[{cvfez[xx ] // Im},{xx, -cvplr, cvplr},DisplayFunction ->  Identity,PlotStyle -> {HuCol1}];

cvp23  =  
            Show[cvp2,  ImageSize  ->  cvims, 
              Epilog  ->  
                cvtru, 
AxesLabel -> {"x [mm]", 
                  None},  
              PlotLabel  ->  
SequenceForm[ 
                  Subscript["E", "y"],  "  [kV/m]"]];

  
cvcm = CartesianMap[cvf,{5  cvlu, 5  cvru, cvdu},{$MachineEpsilon, π (1 - $MachineEpsilon), π (1 - 2  $MachineEpsilon)/cvneq},
 ImageSize  ->  cvims ,PlotStyle -> {HuCol1},
 PlotRange  -> {cvplr {-1, 1}, {-0.1, 
                      2  cvgd }},
\(AxesLabel -> {"x [mm]", "y [mm]"}\),
Epilog->{Thickness[cvtcb],Line[2 {{-cvplr, 0}, {cvplr, 0}}],
Line[{{-cvplr, 
                            cvgd}, {cvc - cvw, cvgd}}],Line[{{ cvplr, cvgd}, {cvc, cvgd}}],
Hue[0],Line[{{cvc, cvgd}, {cvc - cvw, cvgd}}],
Text[SequenceForm["Φ = ", 
                            cvV, " Volt"],{cvmc,1.15 cvgd}]}];

  
Show[cvp1, 
            Graphics[cvtrl], 
 
            ImageSize  -> 
              cvims, 
AxesLabel -> {"x [mm]", 
                None},  
            PlotLabel  ->  
SequenceForm[ 
                Subscript["E", "y"],  "  [kV/m]"], 
            Epilog  ->  cvtrl];

3.2  Variables needed in all pictures - SemInf

hufur[hurgd_, hurb_]  =  
        2.5  +  
          UnitStep[
              hurgd/Abs[hurb]  -  2]  (0.5  hurgd/Abs[hurb]  -  
                1);

hufyr[hurgd_, hurb_]  =  
        3  +  UnitStep[
              3  -  hurgd/Abs[hurb] ]  (1.5  - 
                0.5 hurgd/Abs[hurb] );

huneq = 
        cvhuneq  ;

huneq  specifies the number of equipotential and field lines;  The range of the variable U(=hudu) is bounded by hulu  and 
by huru. This value and the function hufur have been determined by trial and error. This applies also to the function
hufyr rendering a value hufy needed for the plot range in picture hup23. The horizontal plot range in all pictures is 
determined by huplr.

huplr  =   (hugd + Abs[2 hub]);   huims  =  
      500;
hulu  =  Re[hui[-huplr]];
huru  =  hufur[hugd, Abs[hub]] ;
hudu  =  huru/5;
hufy  =  hufyr[hugd, Abs[hub]];

The user may want to improve the pictures drawn below. For this, he may try to obtain better values for 
hulu, huru, hudu, hufy. If he wants to do so it is recommended to remove the semicolons of the last block
at the end of the preceeding subsection, so that the values are printed. Then he may try a slight change of
these values and see what comes out. It is commendable to insert these new values below in the
commands for  hup2, hup3, hup23, hucm, hup1. It is extremely commendable to save this  trial notebook
under a new name beforehand so that the notebook with the orginal specifications remains intact!

hup2 (hup3) gives the electric field at the inside and outside of the upper electrode.
These are then combined to gives figure hup23.  hucm plots the equipotential and electric field lines. 
hup1 shows  the electric field at the bottom electrode. In all four hup pictures the unit of the electric
field strength is  kV/m  since in the definition of the function  hufe in 2.2.3  D = hugd is used, whose value
is defined in mm in  1.1, while we use  MKSA units.

hup2 = ParametricPlot[{huf[ww  +  I π] // Re, 
              hufe[ww +  I π] // Im},{ww,-hulu,0},  DisplayFunction -> Identity,PlotStyle -> {HuCol2}];

hup3 = ParametricPlot[{huf[ww  +  I π] // Re, 
              hufe[ww +  I π] // Im},
{ww,0,hulu},  DisplayFunction -> Identity,PlotStyle -> {HuCol2}];

\(Print["Electric field distribution on upper electrode"]\)
hup23 = Show[hup2,hup3,
DisplayFunction ->  $DisplayFunction,
 PlotRange  ->  {huplr {-1, 1}, 
                        5 huV/hugd {-1, 1}}, ImageSize  ->  huims,Epilog -> { Hue[0], 
                        AbsoluteThickness[
                          2], 
Line[{{hub, -5 huV/
                                  hugd}, {hub, 5  huV/hugd}}]},
AxesLabel -> {"x [mm]", None}, \(PlotLabel\  -> \ 
SequenceForm[Subscript["         ", "y"], "    
          ", Subscript["E", "y"], \ "  [kV/m]"]\)];

Print["Equipotential and electric field lines"]
hucm = CartesianMap[huf,{hulu, huru, hudu},{$MachineEpsilon, 
                  Pi (1 - $MachineEpsilon), π/huneq},
 ImageSize  ->  huims ,
\(AxesLabel -> {"x [mm]", "y [mm]"}\),PlotRange  -> {huplr {-1, 1}, {-0.2, 
                      2  hugd }},PlotStyle -> {HuCol2},Epilog->{Thickness[hutcb],Hue[0],Line[{{-huplr, hugd},{ hub, hugd}}],
Text[
                        SequenceForm["Φ = ", 
                          huV, " Volt"], {- 0.6 huplr, 
                          1.2  hugd }],GrayLevel[0],Line[{{huplr, 0}, {-huplr, 0}}]}]

Print["Electric field distribution on lower electrode"]
hupl = ParametricPlot[{huf[ww] // Re, hufe[ww] // Im},{ww, -hulu, hulu},PlotStyle -> {HuCol2},
ImageSize  ->  huims,Epilog -> { Hue[0], 
                    AbsoluteThickness[
                      2], 
Line[{{hub, -huV/
                            hugd}, {hub, 0}}]}, PlotRange  ->  {huplr {-1, 1}, 
                    huV/hugd {-1, 0}},
AxesLabel -> {"x [mm]", 
                    SequenceForm[
                      Subscript["E", "y"], "  [kV/m]"]}];

3.3  Comparison of the described plane condensors

Print["Electric field distribution on upper electrodes"]
Show[cvp23, hup23, 
        DisplayFunction ->  $DisplayFunction];

    
Print["Equipotential and electric field lines"]
Show[hucm, cvcm, 
        DisplayFunction ->  $DisplayFunction];

    
Print["Electric field distribution on lower electrode"]
Show[cvp1, hupl, Graphics[cvtrl],  
        DisplayFunction ->  $DisplayFunction, 
 
        ImageSize  -> cvims, 
        PlotRange  ->  {cvplr {-1, 1}, 
            cvhuV {-0.1, 
                0}}, 
AxesLabel -> {"x [mm]", 
            None},  PlotLabel  ->  
SequenceForm[ 
            Subscript["E", "y"],  "  [kV/m]"], 
        Epilog  ->  cvtrl];

4.  Field and Potential at a Selected Point

Below the coordinates (xp,yp ) of a point zp may be chosen by the user.The values of the potential and
of the electric field at this point are computed and printed.The equipotential line and the field line passing through that point are shown.

cvhuxp  =  -10 ; cvhuyp  =  
        5 ;

4.1  Strip

cvWp is the point  Wp  corresponding to  zp   computed by the inverse mapping function  W = (f^(-1))(z) = cvi.  
This value is used to compute the complex potential in the W-condensor, which is equal to the potential at zp 
in the real condensor. Wp is also used to compute the complex electric field by the formula given in  2.2.4.  
cvpc (cvec) displays the equipotential (electric field ) line passing through  zp.  The factor  1000  in the second 
Print command (giving the electric field strength)  is  needed since in the definition of the function  cvfe in 2.2.4  
D = cvgd is used, whose value is defined  in mm in 1.1, while we use  MKSA units.

cvxp  =  cvhuxp ; cvyp  =  
        cvhuyp ;

cvWp = cvi[cvxp   + I cvyp];
Print[
        SequenceForm["At point (cvxp,cvyp) = (", cvxp, ", ", 
          cvyp, ") mm : ", "Potential = ", 
          N[Re[-I cvV/π cvWp]] // Chop, " V;"]];
Print[
          SequenceForm["             (", 
            Subscript["E", "x"], ",", 
            Subscript["E", "y"], ") = (", 
            N[Re[1000 cvfe[cvWp]]], ", ", 
            N[Im[1000 cvfe[cvWp]]] ") V/m "]];

cvpc = 
        ParametricPlot[
          Evaluate[{Re[cvf[I Im[cvWp] + s]], 
              Im[cvf[I Im[cvWp] + s]]}], {s, -20, 40}, 
          DisplayFunction -> Identity, 
          PlotStyle -> {HuCol1, Dashing[{ .01}]}];
cvec=ParametricPlot[Evaluate[{Re[cvf[Re[cvWp] + I s]], 
                Im[cvf[Re[cvWp] + I s]]}],{s, 0, Pi},DisplayFunction -> Identity,PlotStyle -> {HuCol1}];
cvpoi=Show[cvpc,cvec,AspectRatio -> Automatic,DisplayFunction -> $DisplayFunction,ImageSize -> 600,TextStyle -> {FontSize -> 16},AxesLabel -> {"x [mm]", "y [mm]"},
 PlotRange  -> {cvplr {-1, 1}, {-0.1, 
                  cvgd 1.3}},Epilog->{Thickness[cvtcb],Line[2 {{-cvplr, 0}, {cvplr, 0}}],
Line[{{-cvplr, 
                        cvgd}, {cvc - cvw, cvgd}}],Line[{{ cvplr, cvgd}, {cvc, cvgd}}],
Hue[0],Line[{{cvc, cvgd}, {cvc - cvw, cvgd}}],
Text[SequenceForm["Φ = ", 
                        cvV, " Volt"],{cvmc,1.15 cvgd}],
AbsolutePointSize[4],Point[{cvxp, cvyp}]}];

4.2  SemInf

huWp is the point  Wp  corresponding to  zp   computed by the inverse mapping function  
W = (f^(-1))(z) = hui.  This value is used to compute the complex potential in the W-condensor,
which is equal to the potential at zp in the real condensor. Wp is also used to compute the
complex electric field by the formula given in  2.2.3.  hupc (huec) displays the equipotential 
(electric field) line passing through  zp.  The factor  1000  in the second Print command 
(giving the electric field strength)  is  needed since in the definition of the function  hufe in 2.2.3  
D = hugd is used, whose value is defined  in mm in  1.1, while we use  MKSA units.

huxp  =  cvhuxp ; huyp  = 
        cvhuyp ;

huWp = hui[huxp  + I huyp];
Print[
        SequenceForm["At point (huxp,huyp) = (", huxp, ", ", 
          huyp, ") mm : ", "Potential = ", 
          N[huV/π Im[huWp]] // Chop, " V;"]];
Print[
        SequenceForm["             (", 
          Subscript["E", "x"], ",", 
          Subscript["E", "y"], ") = (", 
          Re[1000 hufe[huWp]], ", ", 
          Im[1000 hufe[huWp]],  ") V/m "]];

    
hupc = 
        ParametricPlot[
          Evaluate[{Re[huf[I Im[huWp] + s]], Im[huf[I Im[huWp] + s]]}], {s, 
            hulu, huru}, PlotStyle -> {HuCol2, Dashing[{ .01}]}, 
          DisplayFunction -> Identity];
huec = 
        ParametricPlot[
          Evaluate[{Re[huf[Re[huWp] + I s]], Im[huf[Re[huWp] + I s]]}], {s, 
            0, Pi}, DisplayFunction -> Identity, 
          PlotStyle -> {HuCol2}];
hupoi = 
        Show[hupc, huec, AspectRatio -> Automatic, 
          DisplayFunction -> $DisplayFunction, ImageSize -> 600, 
          TextStyle -> {FontSize -> 16}, 
          AxesLabel -> {"x [mm]", "y [mm]"}, 
          PlotRange -> {Max[Abs[huxp], huplr] {-1, 1}, {- .2, 
                2  hugd }}, 
          Epilog -> {Thickness[hutcb], Hue[0], 
              Line[{{-Max[Abs[huxp], huplr], hugd}, { hub, 
                    hugd}}], 
Text[
                SequenceForm["Φ = ", 
                  huV, " Volt"], {- 0.6 huplr, 1.2  hugd }], 
              GrayLevel[0], 
              Line[{{-Max[Abs[huxp], huplr], 0}, {Max[Abs[huxp], huplr], 
                    0}}], AbsolutePointSize[4], 
              Point[{huxp, huyp}]}];

4.3  Comparison

Print["----------- Shifted Voltage Strip ----------"];
Print[
        SequenceForm["At point (cvxp,cvyp) = (", cvxp, ", ", 
          cvyp, ") mm : ", "Potential = ", 
          N[Re[-I cvV/π cvWp]] // Chop, " V;"]];
Print[
        SequenceForm["             (", 
          Subscript["E", "x"], ",", 
          Subscript["E", "y"], ") = (", 
          N[Re[1000 cvfe[cvWp]]], ", ", 
          N[Im[1000 cvfe[
                    cvWp]]] ") V/m "]];
Print["----------- Shifted Semi Infinite ----------"];
Print[
        SequenceForm["At point (huxp,huyp) = (", huxp, ", ", 
          huyp, ") mm : ", "Potential = ", 
          N[huV/π Im[huWp]] // Chop, " V;"]];
Print[
        SequenceForm["             (", 
          Subscript["E", "x"], ",", 
          Subscript["E", "y"], ") = (", 
          Re[1000 hufe[huWp]], ", ", 
          Im[1000 hufe[huWp]],  ") V/m "]];

    
Show[hupoi, cvpoi, 
        DisplayFunction -> $DisplayFunction];

5.   Current Induced on the Electrode for a Selected Trajectory

cvhuxs  =   -20;

The user may select the horizontal coordinate  x = cvhuxs  of a vertical trajectory. The first figure shows the position of this 
trajectory in the real condensor.

5.1 Trajectory - Strip

The next one gives the current ratio induced in the upper left electrode as a function of 
the vertical position  y = cvyt  of the point charge, using W = (f^(-1))(z) =cvi  in the function cvfe for the complex electric field.

cvxs  =   cvhuxs;

\(\(\( (*\(Show[
          Graphics[{Thickness[0.005], \ Line[{{cvxs, 0}, {cvxs, cvgd}}]}], 
          Graphics[{Thickness[cvtcb], 
              Line[cvplr {{\(-1\), 0}, {1, 
                      0}}], 
Line[{{\(-cvplr\), 
                    cvgd}, {cvc - cvw, cvgd}}], 
              Line[{{\ cvplr, cvgd}, {cvc, cvgd}}], 
Hue[
                0], Line[{{cvc, cvgd}, {cvc - cvw, 
                    cvgd}}]}], 
Axes\  -> \ True, \ 
          ImageSize\  -> \ 
            500];\)*) \)\(
\)\(cvpls\  = \ 
        Plot[\(-\ Im[cvfe[cvi[cvxs + \ I\ cvyt]]]\)/cvV, {cvyt, 0, 
            cvgd}, \ 
PlotStyle\  -> \ {HuCol1, 
              Thickness[0.005]}, \ ImageSize\  -> \ 500, 
          DisplayFunction -> \ 
            Identity, 
PlotLabel\  -> \ \ SequenceForm["x = ", 
              cvxs, " mm"], 
AxesLabel\  -> \ \ {"y [mm]", "I/qv [1/mm]"}];\)\)\)

5.2 Trajectory - SemInf

The next one gives the current ratio induced in the upper electrode as a function of 
the vertical position  y = huyt  of the point charge, using W = (f^(-1))(z) = hui  in the function hufe for the complex electric field

huxs  =   cvhuxs;

Show[Graphics[{Thickness[0.005],  Line[{{huxs, 0}, {huxs, hugd}}]}], 
        Graphics[{Thickness[hutcb], 
            Line[Max[Abs[huxs], 
                  huplr] {{-1, 0}, {1, 0}}], 
Hue[0],  
            Line[{{-Max[Abs[huxp], huplr], hugd}, {hub, 
                  hugd}}]}], 
     Axes  ->  
          True,  ImageSize  ->  500];
hupls  =  
        Plot[- Im[hufe[hui[huxs  +  I huyt]]]/huV, {huyt, 0, 
            hugd},  
PlotStyle  ->  {HuCol2, 
              Thickness[0.005]},  
          ImageSize  ->  
            500, 
PlotLabel  ->   SequenceForm["x \
= ", huxs, " mm"], 
AxesLabel  ->   {"y \
[mm]", "ID/qv"}];

5.3 Comparison

Show[Graphics[{Thickness[0.005],  Line[{{huxs, 0}, {huxs, hugd}}]}], 
        Graphics[{Thickness[hutcb], 
            Line[Max[Abs[huxs], 
                  huplr] {{-1, 0}, {1, 0}}], 
Hue[0],  
            Line[{{-Max[Abs[huxp], huplr], hugd}, {hub, 
                  hugd}}]}], 
     Axes  ->  
          True,  ImageSize  ->  500];
Show[hupls, cvpls, 
        DisplayFunction ->  $DisplayFunction];