<< Miscellaneous`PhysicalConstants`

Tin = 1000;
cm = 0.01;
a = 50  cm;
b = 50  cm;
c = 70  cm;
Tpovrch = 55;
To = 20;
α = 5;
σ = StefanConstant[[1]];
S1 = 
        2*((a + d1)*(b + d1) + (a + d1)*(c + d1) + (b + 
                    d1)*(c + d1));
d1 = 3  cm;
λ1 = 0.51;
Rt1 = (d1/(λ1*S1));
Print["Tep. odpor vnitri vyzdivky= ", Rt1, " K/W"]

S2 = 
        2*((a + 2  d1 + d2)*(b + 2  d1 + d2) + (a + 2  d1 + 
                    d2)*(c + 2  d1 + d2) + (b + 2  d1 + d2)*(c + 
                    2  d1 + d2));
dataSibral = {{400, 0.105}, {600, 
            0.163}};
λ2 = 
        Fit[dataSibral, {1, T}, T] /. 
          T -> 500;
\(Print["λSibral= ", λ2, \*"" W.(m^(-1)).(K^(-1))""]\)
Rt2 = (d2/(λ2*S2))

S3 = 
        2*((a + 2  d1 + 2  d2)*(b + 2  d1 + 2  d2) + (a + 2  d1 + 
                    2  d2)*(c + 2  d1 + 2  d2) + (b + 2  d1 + 
                    2  d2)*(c + 2  d1 + 
                    2  d2));
Qk = α*S3*(Tpovrch - To);
ϵ = 0.85;
Qr = ϵ*σ*
          S3*((((Tpovrch + 273))^4) - (((To + 273))^4));
rce = 
        Qk + Qr == ((Tin - Tpovrch)/(Rt1 + Rt2));
dosd2 = Solve[rce] // Last;
Print["d2= ", ((d2 /. dosd2)/cm), " cm"]
Pel = Qk + Qr /. dosd2;
Print["Pel= ", Pel, " W"]