Off[General::"newpkg",General::"obspkg"]

<<Calculus`VectorAnalysis`;
Off[General::"spell1"];
<<Miscellaneous`PhysicalConstants`
<<Miscellaneous`ChemicalElements`

SetCoordinates[Cartesian[x,y,z]];

VacuumPermittivity

EVektor={Ex[y],0,0}
Curl[EVektor]
HVektor={0,0,Hz[y]};
Curl[HVektor]

Thread[Curl[EVektor]==-I*ω*μ*HVektor]
rce1=%[[3]]

Thread[Curl[HVektor]==γ*EVektor]
rce2=%[[1]]

FullSimplify[DSolve[{rce1,rce2,Hz[0]==H0,Ex[0]==E0},{Ex[y],Hz[y]},y][[1]],{μ>0,γ>0,ω>0}]

{rce1,rce2}

rceH=rce1/.D[Solve[rce2,Ex[y]][[1]],y]

dosres=Simplify[DSolve[rceH,Hz[y],y][[1]],{μ>0,γ>0,ω>0}]

(((-1))^(1/4))//N
((1+I)/(Sqrt[2]))//N

dosres/.C[1]->0
%/.y->0

hacko=Hz[y]/.(defres=dosres/.{C[1]->0,C[2]->H0})
ecko=(1/γ)Curl[{0,0,hacko}][[1]]

cm=0.01;
D1=15cm;
T1=20;
T2=500;
delka=0.5;
minuta=60;
tmax=5minuta;
c=400;
m=ρ*V;
ρ=8400;
V=((D1^2)/4)*Pi*delka;
P=Pt/.Solve[Pt*tmax==m*c*(T2-T1)][[1]];
Print["cisty tepelny vykon = ",P," W"]

γCu=57*(10^6);
γBrass=0.28γCu;
ω=2*Pi*f0;
μ0=N[VacuumPermeability][[1]];
aBrass=(Sqrt[(2/(ω*γBrass*μ0))]);
Solve[aBrass==0.12*D1]
f=50;