I*{4-(7./(Sqrt[3]))+(3./(Sqrt[3])),-8+(5./(Sqrt[3]))-(3./(Sqrt[3])),(7./(Sqrt[3]))-(5./(Sqrt[3]))}

epsilonPanelVen=1;
epsilonPanelZem=1;
alfa=6;
qOdSlunce=1000.;
PmaxPanel=100;
plochaPanel=2;
Tnebe=20+273;
Tzem=30+273;
Tvzduch=20+273;
uhelDopadu=0*Degree;
sigma=5.66*(10^(-8));
ucinnost=(PmaxPanel/(plochaPanel*qOdSlunce));

PSlunce=qOdSlunce*plochaPanel*Cos[uhelDopadu]*epsilonPanelVen;
PKonvekce=alfa*plochaPanel*2*(T-Tvzduch);
PSalaniDoZem=epsilonPanelZem*sigma*plochaPanel*((T^4)-(Tzem^4));
PSalaniDoNebe=epsilonPanelVen*sigma*plochaPanel*((T^4)-(Tnebe^4));
Pel=qOdSlunce*Cos[uhelDopadu]*plochaPanel*ucinnost;


res=Solve[PSlunce==PKonvekce+PSalaniDoZem+PSalaniDoNebe+Pel,T];
Tvysl=Select[T/.res,Re[#]>0&&Im[#]==0&][[1]]-273;
Print["Tvysl= ",Tvysl," °C"]

ClearAll[alfa]
Tpec=1300+273;
dIzolace
Tokoli=20;
alfaVnejsi=6;
epsilonPec=0.9;