SetDirectory["c:\data"];
<<  ./cteni1.m

\( (*\ 
      field\ concetrator\ parameters\ *) 
 (*\ \(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(----\)\)\)\)\)\)\)\)\)\)\)\)\)\)\ *) 
\(DeleniRfc = 1;\)\)
DeleniZfc = 1;
R1fc = 15;
R2fc = 20;
\(\(Zfc = 5;\)
 (*\ 
      charge\ parameters\ *) 
 (*\ \(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(----\)\)\)\)\)\)\)\)\)\)\)\)\)\)\ *) \)
DeleniZcharge = 2;
DeleniRcharge = 2;
R1charge = 10;
R2charge = 14;
\(\(Zcharge = 5;\)
 (*\ 
      coil\ parameters\ *) 
 (*\ \(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(----\)\)\)\)\)\)\)\)\)\)\)\)\)\)\ *) \)
DeleniZcoil = 1;
DeleniRcoil = 1;
R1coil = 20;
R2coil = 30;
Zcoil = 5;

Deleni2Civek[vzdalenostCivek_, 
R11_, R21_, 
        Z1_, 
R12_, R22_, 
        Z2_, 
deleni1R_,  
        deleni1Z_, 
deleni2R_,  
        deleni2Z_
] := 
      Module[{dr1, dr2, dz1, dz2, i, j, R1aR2a, Z1aZ2a, RZ1, R1bR2b, Z1bZ2b, 
          RZ2}, 
dr1 = (R21 - R11)/
            deleni1R; 
dz1 = 
          Z1/deleni1Z; 
dr2 = (R22 - R12)/
            deleni2R; 
dz2 = 
          Z2/deleni2Z; 

R1aR2a = 
          Table[N[{R11 + i*dr1, R11 + (i + 1)*dr1}], {i, 0, 
              deleni1R - 1}]; 
Z1aZ2a = 
          Table[N[{i*dz1, (i + 1)*dz1}], {i, 0, 
              deleni1Z - 1}]; 

RZ1 = 
          Flatten[Table[
Flatten[{R1aR2a[[i]], 
                  Z1aZ2a[[j]]}]
, {j, 1, 
                Length[Z1aZ2a]}, {i, 1, Length[R1aR2a]}], 
            1]; 
R1bR2b = 
          Table[N[{R12 + i*dr2, R12 + (i + 1)*dr2}], {i, 0, 
              deleni2R - 1}]; 
Z1bZ2b = 
          Table[N[{vzdalenostCivek + i*dz2, 
                vzdalenostCivek + (i + 1)*dz2}], {i, 0, 
              deleni2Z - 1}]; 

RZ2 = 
          Flatten[Table[
Flatten[{R1bR2b[[i]], 
                  Z1bZ2b[[j]]}]
, {j, 1, 
                Length[Z1bZ2b]}, {i, 1, Length[R1bR2b]}], 
            1]; 
{RZ1
, RZ2}]

pomcoil = 
        Deleni2Civek[0, 
R1coil, R2coil, 
          Zcoil, 
R1coil, R2coil, 
          Zcoil, 
DeleniRcoil, 
          DeleniZcoil, 
DeleniRcoil, 
          DeleniZcoil];
Lcoil = Lab[pomcoil[[1, 1]], pomcoil[[2, 1]]]

pomfc = 
        Deleni2Civek[0, 
R1fc, R2fc, 
          Zfc, 
R1fc, R2fc, 
          Zfc, 
DeleniRfc, 
          DeleniZfc, 
DeleniRfc, DeleniZfc];
Lfc = Lab[pomfc[[1, 1]], pomfc[[2, 1]]]

(* ------------------------------ *) 
 (* 
      coil - field concentrator induction *) 
pomcoilfc = 
        Deleni2Civek[10, 
R1coil, R2coil, 
          Zcoil, 
R1fc, R2fc, 
          Zfc, 
DeleniRcoil, 
          DeleniZcoil, 
DeleniRfc,  
          DeleniZfc];
Lcoilfc = Lab[pomcoilfc[[1, 1]], pomcoilfc[[2, 1]]]

(* ------------------------------ *) 
 (* 
      charge - field concentrator induction *) 
pomchargefc = 
        Deleni2Civek[5, 
R1charge, R2charge, 
          Zcharge, 
R1fc, R2fc, 
          Zfc, 
DeleniRcharge, 
          DeleniZcharge, 
DeleniRfc,  
          DeleniZfc];
Lchargejfc = 
      Table[Lab[pomchargefc[[1, i]], pomchargefc[[2, 1]]], {i, 1, 
          Length[pomchargefc[[1]]]}]

pomchargecoil = 
        Deleni2Civek[5, 
R1charge, R2charge, 
          Zcharge, 
R1coil, R2coil, 
          Zcoil, 
DeleniRcharge, 
          DeleniZcharge, 
DeleniRcoil, 
          DeleniZcoil];
Lchargejcoil = 
      Table[Lab[pomchargecoil[[1, i]], pomchargecoil[[2, 1]]], {i, 1, 
          Length[pomchargecoil[[1]]]}]

pomchargeji = 
        Deleni2Civek[0, 
R1charge, R2charge, 
          Zcharge, 
R1charge, R2charge, 
          Zcharge, 
DeleniRcharge, 
          DeleniZcharge, 
DeleniRcharge, 
          DeleniZcharge];
Lchargeji = 
      Table[Lab[pomchargeji[[1, i]], pomchargeji[[2, j]]], {i, 1, 
          Length[pomchargeji[[1]]]}, {j, 1, 
          Length[pomchargeji[[2]]]}]