cp = 470;
cm = 0.01;
el = 3;
ro = 7850;
γFe = (0.000000245^(-1));
d = 0.3;
m = Pi*((d^2)/4)*el*ro;
minuta = 60;
cas = 20  minuta;
ΔT = 1000;
energie = 1.  ro*cp*ΔT;
vykon = (energie/cas);
Pt = 1.5*vykon;

<< Miscellaneous`PhysicalConstants`

μ = (VacuumPermeability // 
            N)[[1]]

ω = 100.  Pi;
γCu = 56*(10^6);
aFe = (Sqrt[(2/(ω*γFe*μ))]);
aCu = (Sqrt[(2/(ω*γCu*μ))]);

dingot = 0.3;
d2 = dingot - aFe;
R2 =  Pi*d2/(γFe*aFe*el)

mezera = 3  cm;
dinduktor = dingot + 2*mezera;
d1 = dinduktor + aCu;
R1 = (n1^2)* Pi*d1/(γCu*aCu*el)

\(\(Print[\*""\!\(d2\/el\)= "", 
        d2\/el];\)\)
α2 = 0.97;
L2 =  Pi*(((d2^2)*(1^2))/(4  el))*μ*α2;

\(\(Print[\*""\!\(d1\/el\)= "", 
        d1\/el];\)\)
α1 = 0.96;
L1 =  Pi*(((d1^2)*(n1^2))/(4  el))*μ*α1;

datF = {{0.2, 0.2, 1, 0.966}, {0.2, 
            0.2, 1.1, 0.880}, {0.3, 0.2, 1, 0.930}, {0.3, 0.3, 1, 
            0.931}};

    nahF = Fit[datF, {1, p1, p2, pl}, {p1, p2, pl}]

    naseF = 
      nahF /. {p1 -> 0.12, p2 -> 0.09, pl -> 1}

F = naseF
M =  Pi*(((d2^2)*1*n1)/(4  el))*μ*F

data = 
        Table[{x, Sin[x]}, {x, -Pi, Pi, 
            Pi/10. }];
intOrd = 2;
nahSin[
          x_] = Interpolation[data, InterpolationOrder -> intOrd][
          x];
Plot[nahSin''[x], {x, -Pi, Pi}]

datF = {{0.2, 0.2, 1, 0.966}, {0.2, 
            0.2, 1.1, 0.880}, {0.3, 0.2, 1, 0.930}, {0.3, 0.3, 1, 
            0.931}};

    nahF = Fit[datF, {1, p1, p2, pl}, {p1, p2, pl}]

    naseF = 
      nahF /. {p1 -> 0.12, p2 -> 0.09, pl -> 1}

U = 400;
rc1 = 
        U == 
          R1*i1 + L1*I*ω*i1 + 
            M*I*ω*i2;
rc2 = 
        0 == 
          R2*i2 + L2*I*ω*i2 + 
            M*I*ω*i1;
rc3 = R2*Abs[i2]^2 == Pt;
res = Solve[{rc1, rc2}, {i1, 
            i2}][[1]]
rcevykon = rc3 /. res

FindRoot[rcevykon, {n1, 10, 1000}]

(rcevykon /. Abs[tt_] -> (Sqrt[Re[tt]^2 + Im[tt]^2]) // Simplify) // 
      Solve

dosn1 = {n1 -> 35};
dosproudy = res /. dosn1

R2*Abs[i2]^2 /. dosproudy /. dosn1

Pt