<<PhysicalConstants`

mm=(10^(-3));
dingot=200mm;
μ=N[VacuumPermeability][[1]];
γ=50*(10^6);
f=50;
ω=2.*Pi*f;
a=(Sqrt[(2/(ω*γ*μ))]);
(*Print["a= ",(10^3)a,"  mm"];*)
d2=dingot-a;
el2=1000mm;
S2=el2*a;
R2=((Pi*d2)/(S2*γ));
(*Print["R2= ",(10^6)R2,"  μΩ"]*)

mezera=10mm;
dinduktor=dingot+2.*mezera;
d1=dinduktor+a;

el1=el2;
S1=el1*a;
R1=1.25((Pi*d1)/(S1*γ))*(n^2);

SL1=Pi*((d1^2)/4);
L1nekonecna=μ*(n^2)*(SL1/el1);
α1=0.93;
L1=α1*L1nekonecna;
SL2=Pi*((d2^2)/4);
L2nekonecna=μ*(SL2/el2);
α2=0.95;
L2=α2*L2nekonecna;
F=0.95;
M1=F*n*L2;

hustota=8500;
c=383;
T1=20;
T2=800;
V=Pi*((dingot^2)/4.)*el2;
m=hustota*V;
energie=m*c*(T2-T1);
t=5*60;
Print["vykon= ",Pideal=(energie/t),"  W"];
P=1.2*Pideal;