<< Miscellaneous`PhysicalConstants`
mH2015 = 770*0.15
mdreva15 = 770*(1 - 0.15)

voda1 = 
        mvoda18 /. Solve[
              (mvoda18/(mdreva15 + mvoda18)) == 
                0.18][[1]];
voda2 = 
        mvoda08 /. Solve[
              (mvoda08/(mdreva15 + mvoda08)) == 
                0.08][[1]];
Δvoda1kubik = %% - %;

V = 3*1.5*0.01;

Δm = Δvoda1kubik*V

r = 2.5*(10^6);
cpdt = 100*4186;
%/%%

ΔW = (r + cpdt)*Δm

VacuumPermittivity[[1]]

Pstr = Δ(W/(2*3600));
dos = {ω -> 2*Pi*f, ϵ0 -> 
            VacuumPermittivity[[1]], ϵrtgδ -> 
            3*(10^(-2)), U -> 1000, 
          d -> 0.015};
rce = 
        Pstr == 
            V*ω*ϵ0*ϵrtgδ*((((U/d)))^2) /. 
          dos;
resf = Solve[rce, f][[1]]
Solve[2  tgδ == 0.03, 
        tgδ];
Q = (Pstr/0.015)
Cdesk = ((VacuumPermittivity[[1]]*2*1.5*3)/0.015)
L = (1/(Cdesk*(((2*Pi*f))^2))) /. resf
dosr = (Solve[Pstr == ((U^2)/R), R] /. 
            dos)[[1]]

doscelk = 
        Union[dosr, resf, 
          dos, {cap -> 5.3*(10^(-9)), 
            L -> 212*(10^(-9))}];

par[z1_, z2_] := ((z1*z2)/(z1 + z2));
1000/
      Abs[par[(1/(I*ω*cap)), R]/(par[(1/(I*ω*cap)), R] + I*ω*L) //. doscelk]