a=0.4;
b=0.5;
c=0.3;
cm=0.01;
λPlastizol=0.44;
dPlastizol=2cm;
λSibral=0.11;
SStredniPlastizol=0.5(2*(a*b+b*c+a*c)+2*((a+2*dPlastizol)*(b+2*dPlastizol)+(b+2*dPlastizol)*(c+2*dPlastizol)+(a+2*dPlastizol)*(c+2*dPlastizol)));
SStredniSibral=0.5*(2*((a+2*dPlastizol)*(b+2*dPlastizol)+(b+2*dPlastizol)*(c+2*dPlastizol)+(a+2*dPlastizol)*(c+2*dPlastizol))+2*((a+2*dPlastizol+2*dSibral)*(b+2*dPlastizol+2*dSibral)+(b+2*dPlastizol+2*dSibral)*(c+2*dPlastizol+2*dSibral)+(a+2*dPlastizol+2*dSibral)*(c+2*dPlastizol+2*dSibral)));
Svnejsi=2*((a+2*dPlastizol+2*dSibral)*(b+2*dPlastizol+2*dSibral)+(b+2*dPlastizol+2*dSibral)*(c+2*dPlastizol+2*dSibral)+(a+2*dPlastizol+2*dSibral)*(c+2*dPlastizol+2*dSibral));
α=6;
Qkonvekce=Svnejsi*α*(TpovrchPece-Tokoli);
ϵ=0.3;
σ=5.67*(10^(-8));
Qsalani=ϵ*σ*Svnejsi*((((TpovrchPece+273))^4)-(((Tokoli+273))^4));
Tokoli=20;
Tpec=1200;
RthPlastizol=(dPlastizol/(λPlastizol*SStredniPlastizol));
RthSibral=(dSibral/(λSibral*SStredniSibral));
rce=((Tpec-TpovrchPece)/(RthPlastizol+RthSibral))==Qkonvekce+Qsalani;
rce1=rce/.TpovrchPece->50;
dos1=Solve[rce1][[-1]];
Print["minimalni prikon pece= ",Qkonvekce+Qsalani/.dos1/.TpovrchPece->50]

rce2=rce/.dSibral->0.1;
dos2=Solve[rce2][[-1]]
Print["minimalni prikon pece= ",Qkonvekce+Qsalani/.dos2/.dSibral->0.1]

ρPlastizol=500;
cPlastizol=1140;
ρSibral=160;
cSibral=cPlastizol;
prechod[odkud_,kam_]:=Module[{a,b,resab,strmost=4000},
resab=Solve[{a-b==odkud,a+b==kam},{a,b}][[1]];
a+b*Tanh[strmost*(x-dPlastizol)]/.resab
];
pom=prechod[1,2];
Plot[pom,{x,0,2*dPlastizol},AxesOrigin->{0,0}];

rocepe=prechod[ρPlastizol*cPlastizol,ρSibral*cSibral]
λ=prechod[λPlastizol,λSibral]

rceNabeh=rocepe*D[T[x,t],t]==D[λ*D[T[x,t],x],x]