<< Miscellaneous`PhysicalConstants`
Off[General::"spell1", Solve::"ifun"];

ω = 314;
Dinduktor = 0.4;
Dingot = 0.3;
Linduktor = 4;
Lingot = 4;
ρFe = 1.2*(10^(-6));
ρCu = 2*(10^(-8));
μ0 = (VacuumPermeability // 
              N)[[1]];
aFe = (Sqrt[((2*ρFe)/(ω*μ0))]);
aCu = (Sqrt[((2*ρCu)/(ω*μ0))]);
Print["aFe= ", 100  aFe, " cm"]
Print["aCu= ", 100  aCu, " cm"]

R2 = Ringot = ((ρFe*Pi*(Dingot - aFe))/(aFe*Lingot));
Print["R2= ", (10^6) 
        Ringot, " μΩ"]

R1 = Rinduktor = 
          1.2  (n^2)*((ρCu*Pi*(Dinduktor + aCu))/(aCu*Linduktor));
Print["R1= ", (10^6) 
        Rinduktor, " μΩ"]

\(\(Print[\*""\!\(D\/el\)= "", \(Dingot - aFe\)\/Lingot];\)\)
α2 = 0.98;
L2 = α2*μ0*(1^2)*((Pi*(((Dingot - aFe))^2))/(4*Lingot));

\(\(Print[\*""(D/el)= "", \(Dinduktor + aCu\)\/Linduktor];\)\)
α1 = 0.96;
L1 = α1*μ0*
        (n^2)*((Pi*(((Dinduktor + aCu))^2))/(4*Linduktor))

F =  .98;
M = 
      F*μ0*1*
        n*((Pi*(((Dinduktor + aCu))^2))/(4*Linduktor))

j = I;
U1 = 400;
P2zadane = (10^6);
P2 = Abs[I2]^2*R2;
U2 == j*ω*(L2*I2 + M*I1);
I2 == U2/R2;
U1 == 
        R1*I1 + j*ω (L1*I1 + M*I2);
res = 
        Solve[{U2 == j*ω*(L2*I2 + M*I1), I2 == U2/R2, 
              U1 == R1*I1 + j*ω (L1*I1 + M*I2)}, {I1, I2, 
              U2}] // Chop;
rceVykonu = 
      P2zadane == P2 /. res[[1]]
Solve[rceVykonu][[2]]

dosN = {n -> 70};
resN = 
      Solve[{U2 == j*ω*(L2*I2 + M*I1), I2 == U2/R2, 
              U1 == R1*I1 + j*ω (L1*I1 + M*I2)} /. 
            dosN, {I1, I2, U2}] // Chop
P2Skutecny = P2 /. resN[[1]]

S1 = U1*Conjugate[I1] /. resN[[1]]
P1 = Re[S1]
Q1 = Im[S1]