(2000./(400^2))*(1-I*Tan[ArcCos[0.7]])

(3000./(400^2))*(1+I*Tan[ArcCos[0.8]])

(10000./(400^2))*(1-I*Tan[ArcCos[0.95]])

kompenzace={-I*0.0141-I*(0.0125/(Sqrt[3]))+I*(0.0625/(Sqrt[3])),+I*0.0128+I*(0.0188/(Sqrt[3]))-I*(0.0625/(Sqrt[3])),I*0.0205-I*(0.0188/(Sqrt[3]))+I*(0.0125/(Sqrt[3]))}

C12=C12/.Solve[kompenzace[[1]]==I*314*C12,C12]

C13=C13/.Solve[kompenzace[[3]]==I*314*C13,C13]
L23=L23/.Solve[kompenzace[[2]]==(1/(I*L23*314)),L23]

PVykonOdSlunceNaPlochu=1000;
fi=70°;
PelNaPlochu=150;
epsilon=0.8;
plocha=3;
TeplotaNebe=-30+273;
TeplotaZem=5+273;
TeplotaVzduch=15+273;
sigma=5.66*(10^(-8));
alfa=18;

PodSlunce=PVykonOdSlunceNaPlochu*plocha*Cos[fi];
Pel=PelNaPlochu*plocha*Cos[fi];
PsalaniDoNebe=sigma*epsilon*plocha*((T^4)-(TeplotaNebe^4));
PsalaniDoZeme=sigma*epsilon*plocha*((T^4)-(TeplotaZem^4));
Pkonvekce=alfa*plocha*(T-TeplotaVzduch);

res=T/.Solve[PodSlunce==Pel+PsalaniDoNebe+PsalaniDoZeme+Pkonvekce,T][[4]];
res2=res-273