Remove["Global`*"]//Quiet;

(*m1, c1 = primár fototermiky (T1->T2)
m2, c2 = primár zemního tepelnýho čerpadla (T5->T6)
m3, c3 = sekundár fototermiky (T8 -> T3)
m4, c3 = vytápění (T8 ->T7)
m5, c3 = sekundár tepelnýho čerpadla (T8->T4)
Tin - teplota v domu
Tout - teplota okolí
Q1 teplo z primáru do sekundáru fototermiky
Q2 teplo do sekundáru tepelného čerpadla
Q3 teplo do primáru tepelného čerpadla*)

(*FOTOTERMIKA*)
fototermika1=Psout==α*Skol*(((T1+T2)/2)-Tout)+σ*ϵ*Skol*((((((T1+T2)/2)+273))^4)-(((Tout+273))^4));
fototermika2=Pin-Pout==m1*c1*(T2-T1);

(*Výměník fototermiky*)
fototermika3=m1*c1*T2+m3*c3*T8==m1*c1*T1+m3*c3*T8;
fototermika4=Q1==m1*c1*(T2-T1);
fototermika5=Q1==kS1*(((T2-T3)-(T1-T8))/(Log[((T2-T3)/(T1-T8))]));

(*ZÁSOBNÍK*)
zasobnik=Pel1+T3*c3*m3+T4*c3*m5+T7*c3*m4==c3*T8*(m3*m5*m4);

(*Tepelné čerpadlo*)
cerpadlo1=m2*c2*T6+m5*c3*T8*Pel2==m5*c3*T4+m2*c2*T5;
cerpadlo2=ϵ*Pel2==Pel2+m2*c2*(T6-T5);

(*Zemní výměník tepelného čerpadla*)
cerpadlo3=m2*c2*(T6-T5)==Q2;
cerpadlo4=Q2==kS*(((Tg-T6)-(Tg-T5))/(Log[((Tg-T6)/(Tg-T5))]));

(*Dům*)
dum1=m4*c3*T8==Prad+m4*c3*T7;
dum2=Psz+Prad==Pout;
dum3=Prad==krad*(((T8-T2)/2)-Tin);
dum4=Pout==kdum*(Tin-Tout);

rce={fototermika1,fototermika2,fototermika3,fototermika4,fototermika5,zasobnik,cerpadlo1,cerpadlo2,cerpadlo3,cerpadlo4,dum1,dum2,dum3,dum4};