lambda=0.18;
ϵ=0.65;
sigma=5.66*(10^(-8));
TSalani=25;
TVzduch=20;
TPovrch=36.6;
TPec=1250;
α=8;
rozmerPece=0.5;
tloustaPenosamot=0.02;

rce=(TPec-TPovrch)*(lambda/diz)*6*(rozmerPece^2)==α*6*(((rozmerPece+2*tloustaPenosamot+2*diz))^2)*(TPovrch-TVzduch)+sigma*ϵ*(((rozmerPece+2*tloustaPenosamot+2*diz))^2)*((TPovrch^4)-(TSalani^4));
k=FindRoot[rce,{diz,2}]

(TPec-TPovrch)*(lambda/diz)*6*(rozmerPece^2)/.k