<< Calculus`VectorAnalysis`
SetCoordinates[Cylindrical[r, fi, z]];
Off[General::"spell1"];

evect = {0, 0, Ez[r]};
bvect = {0, Bfi[r], 0};
hvect = (1/μ)*bvect;

rce1 = Curl[evect][[2]] == -I*ω*
            bvect[[2]];
rce2 = Curl[hvect][[3]] == γ*evect[[3]] // 
          Simplify;

dosB = Solve[rce1, Bfi[r]][[1]];
dos = Union[dosB, D[dosB, r]];
rce = rce2 //. dos // Simplify

dosEzpom = DSolve[rce, Ez[r], r][[1]]

Plot[BesselJ[0, x], {x, 0, 1}, 
        AxesLabel -> {"x", "BJ[0,x]"}, 
        PlotStyle -> {Hue[0], 
            Thickness[0.008]}];
Plot[BesselY[0, x], {x, 0, 1}, 
        AxesLabel -> {"x", "BJ[0,x]"}, 
        PlotStyle -> {Hue[0], Thickness[0.008]}];

dosEz = 
        dosEzpom /. {C[2] -> 0, 
            C[1] -> 1};
try = {γ -> 57*(10^6), μ -> 
            4*Pi*(10^(-7))};
cm = 0.01;
rmax = 1  cm;
absEz[f_] = 
        Abs[Ez[r] //. 
            Union[try, {ω -> 2*Pi*f}, 
              dosEz]];
pl[f_] := 
        Plot[{absEz[f], 0}, {r, -rmax, rmax}, PlotRange -> All, 
          PlotStyle -> {{Hue[0], Thickness[0.008]}, Thickness[0.001]}, 
          AxesLabel -> {"r[m]", "Abs[Ez[r]]"}];
pl[150];

ClearAll[proud];
proud[f_] := γ*NIntegrate[2*Pi*r*absEz[f], {r, 0, rmax}] /. 
          try;
ztraty[f_] := γ*NIntegrate[2*Pi*r*absEz[f]^2, {r, 0, rmax}] /. 
        try
odpor[f_] := ztraty[f]/proud[f]^2;
odporSS = odpor[0];
koefodporu[f_] := odpor[f]/odporSS;
\(\(Plot[{koefodporu[f], 0}, {f, 0, 5000}, PlotRange -> All, 
        PlotStyle -> {{Hue[0], Thickness[0.008]}, Thickness[0.001]}, 
        AxesLabel -> {"f[Hz]", \*""\!\(R[f]\/Rss\)""}];\)\)
\(\(Plot[{koefodporu[f], 0}, {f, 0, 500}, PlotRange -> All, 
        PlotStyle -> {{Hue[0], Thickness[0.008]}, Thickness[0.001]}, 
        AxesLabel -> {"f[Hz]", \*""\!\(R[f]\/Rss\)""}];\)\)

\(\(\(
\)\(\(<< 
        Miscellaneous`ChemicalElements`;\)

    \(k = 15;\)

    \(f0 = 500;\)

    \(tmax = 2*10\^3\/\(k\^2*f0\);\)

    \(Qv = \(k*γ*absEz[f]\^2 /. try\) /. 
          f -> f0;\)

    \(prvek = Copper;\)

    \(λ = \(ThermalConductivity[
            prvek]\)[\([1]\)];\)

    \(ρ = \(Density[prvek]\)[\([1]\)];\)

    \(cp = 10\^3*\(SpecificHeat[prvek]\)[\([1]\)];\)

    \(rce = ρ*cp*D[T[r, t], t] == 
          Div[λ*Grad[T[r, t]]] + Qv;\)

    \(To = 20;\)

    \(α = 10;\)

    \(ϵ = 0.8;\)

    \(σ = 5.67*10\^\(-8\);\)

    \(rmin = 0.01  cm;\)

    \(poc = T[r, 0] == To;\)

    \(okr1 = 0 == D[T[r, t], r] /. r -> rmin;\)

    \(qα = α*\((T[rmax, t] - To)\);\)

    \(qϵ = ϵ*σ*\((\((T[rmax, t] + 273)\)\^4 - \((To + 273)\)\^4)\);\)

    \(okr2 = 
        qα + qϵ == D[T[r, t], r] /. 
          r -> rmax;\)

    \(resT = \(NDSolve[{rce, poc, okr1, okr2}, 
            T[r, t], {r, rmin, rmax}, {t, 0, 
              tmax}]\)[\([1]\)];\)

    \(teplota[r_, t_] = T[r, t] /. resT;\)

    \(Plot3D[teplota[r, t], {r, rmin, rmax}, {t, 0, tmax}, 
        AxesLabel -> {"r[m]  ", "  t[s]", "T[r,t] "}, 
        PlotRange -> All];\)


    \)\)\)