ClearAll["Global`*"];
$HistoryLength=2;
SetDirectory[NotebookDirectory[]];

(*************************************TEPELNÉ ZTRÁTY*************************************)
(*------------------------PARAMETRY PRO CELÝ OBJEKT------------------------*)
(*-----Teploty-----*)
(Subscript[θ, e])={-15,-7,2,7}; (*venkovní výpočtová teplota*)
(Subscript[θ, em])=3; (*průměrná venkovní výpočtová teplota*)
(Subscript[θ, o])=21; (*teplota obytných místností*)
(Subscript[θ, koup])=24;(*teplota koupelen*)
(Subscript[θ, o1])=18; (*teplota obytných místností vytápěných na nižší teplotu*)
(Subscript[θ, t])= 15; (*teplota technické místnosti*)
(Subscript[θ, u])=-8; (*teplota nevytápěné místnosti*)

(*-----Součinitele prostupu tepla-----*)
(Subscript[ΔU, tb])=0.1; (*přírustek (Subscript[U, k]) daný tepelnými mosty*)
Uokna=1.3; (*souč. prostupu oken*)
Udveree=1.5; (*souč. prostupu venkovních dvěří*)
Udverei = 2.6; (*souč. prostupu vnitřních dvěří*)
Udvereg=1.7; (*souč. prostupu garážních dveří*)
UzedB1=0.19; (*souč. prostupu konstrukcí B1*)
UzedB2=0.3; (*souč. prostupu konstrukcí B2*)
UzedC1=0.42; (*souč. prostupu konstrukcí C1*)
UzedC2=0.29; (*souč. prostupu konstrukcí C1*)
Ustrecha=0.2; (*souč. prostupu konstrukcí střechy*)
Upodlaha=0.28; (*souč. prostupu konstrukcí podlahy*)

(*-----Další konstanty-----*)
k=0.34; (*konstanta pro větrání k=ρ*c, při uvažování objemového toku v (m^3)/h*)
(Subscript[n, o]) = 0.3;(*min. výměna vzduchu pro obytné místnosti*)
(Subscript[n, k])=0.6;(*min. výměna vzduchu pro kuchyň a koupelnu s okny*)
(Subscript[e, k])=1; (*korekční činitel*)
(Subscript[f, g1])=1.45; (*obvyklá hodnota*)
(Subscript[G, w])=1; (*méně jak 1m pod podlahou není spodní voda*)
(Subscript[A, g])=17.51+1.88+15.56+48.15+27.51+15.39+10.4+9.74+16.6+7.9;(*Plocha podlhové konstrukce*)
P=14.6+17.9+14.6+17.9-4*0.285;(*Obvod podlahové konstrukce*)
B= ((Subscript[A, g])/(0.5*P));(*Charakteristický parametr podlahové konstrukce*)
ekvi=Interpolation[{{0.5,0.28},{0.25,0.17}},InterpolationOrder->1];
(Subscript[U, ekvi]) = ekvi[0.33];(*Ekvivalentní souč. prostupu tepla pro B stanoven interpolací*)

(*------------------------PARAMETRY PRO MÍSTNOSTI------------------------*)
(*-----Místnost 1.01 + 1.02-----*)
(*plocha podlahy*)
(Subscript[Ap, 1])=17.51+1.88;
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 1])=1.3*2.2;
(Subscript[Ae2, 1])=(2.377+0.288)*2.6-(Subscript[Ae1, 1]);
(Subscript[Ae3, 1])=(Subscript[Ap, 1]);
(Subscript[Uke1, 1])=Udveree;
(Subscript[Uke2, 1])=UzedB1;
(Subscript[Uke3, 1])=Ustrecha;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 1])=6.36*2.6;
(Subscript[Aeu2, 1])=0;
(Subscript[Ukeu1, 1])=UzedC2;
(Subscript[Ukeu2, 1])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 1])=(Subscript[θr2, 1])=(Subscript[θr3, 1])=(Subscript[θ, o]);
(Subscript[Aj1, 1])=2*0.8*1.97;
(Subscript[Aj2, 1])=(4.375+2.887)*2.6-(Subscript[Aj1, 1]);
(Subscript[Aj3, 1])=(1.561+0.275)*2.6;
(Subscript[Ukj1, 1])=Udverei;
(Subscript[Ukj2, 1])=UzedB2;
(Subscript[Ukj3, 1])=UzedC1;
(Subscript[θr4, 1])=(Subscript[θr5, 1])=(Subscript[θr6, 1])=(Subscript[θr7, 1])=(Subscript[θr8, 1])=0;
(Subscript[Aj4, 1])=(Subscript[Aj5, 1])=(Subscript[Aj6, 1])=(Subscript[Aj7, 1])=(Subscript[Aj8, 1])=0;
(Subscript[Ukj4, 1])=(Subscript[Ukj5, 1])=(Subscript[Ukj6, 1])=(Subscript[Ukj7, 1])=(Subscript[Ukj58, 1])=0;

(*-----Místnost 1.03-----*)
(Subscript[Ap, 2])=15.56;(*plocha podlahy*)
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 2])=1.68*0.96+1*2.2;
(Subscript[Ae2, 2])=(1.775+4.524+3.611)*2.6-(Subscript[Ae1, 2]);
(Subscript[Ae3, 2])=(Subscript[Ap, 2]);
(Subscript[Uke1, 2])=Uokna;
(Subscript[Uke2, 2])=UzedB1;
(Subscript[Uke3, 2])=Ustrecha;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 2])=(Subscript[Aeu2, 2])=0;
(Subscript[Ukeu1, 2])=(Subscript[Ukeu2, 2])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 2])=(Subscript[θ, o1]);
(Subscript[Aj1, 2])=(1.561+0.275)*2.6;
(Subscript[Ukj1, 2])=UzedC1;
(Subscript[θr2, 2])=(Subscript[θr3, 2])=(Subscript[θr4, 2])=(Subscript[θr5, 2])=(Subscript[θr6, 2])=(Subscript[θr7, 2])=(Subscript[θr8, 2])=0;
(Subscript[Aj2, 2])=(Subscript[Aj3, 2])=(Subscript[Aj4, 2])=(Subscript[Aj5, 2])=(Subscript[Aj6, 2])=(Subscript[Aj7, 2])=(Subscript[Aj8, 2])=0;
(Subscript[Ukj2, 2])=(Subscript[Ukj3, 2])=(Subscript[Ukj4, 2])=(Subscript[Ukj5, 2])=(Subscript[Ukj6, 2])=(Subscript[Ukj7, 2])=(Subscript[Ukj8, 2])=0;

(*-----Místnost 1.04-----*)
(Subscript[Ap, 3])=48.15;(*plocha podlahy*)
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 3])=2*(2.180*2.360+2.280*0.4+1.8*2.38);
(Subscript[Ae2, 3])=(2.5+7.174+1.712)*3.65-(Subscript[Ae1, 3]);
(Subscript[Ae3, 3])=(Subscript[Ap, 3]);
(Subscript[Uke1, 3])=Uokna;
(Subscript[Uke2, 3])=UzedB1;
(Subscript[Uke3, 3])=Ustrecha;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 3])=(Subscript[Aeu2, 3])=0;
(Subscript[Ukeu1, 3])=(Subscript[Ukeu2, 3])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 3])=(Subscript[θr2, 3])=(Subscript[θ, o1]);
(Subscript[Aj1, 3])=0.8*1.97;
(Subscript[Aj2, 3])=(4.375+1.3+0.8)*2.6-(Subscript[Aj1, 3]);
(Subscript[Ukj1, 3])=Udverei;
(Subscript[Ukj2, 3])=UzedC1;
(Subscript[θr3, 3])=(Subscript[θr4, 3])=(Subscript[θr5, 3])=(Subscript[θr6, 3])=(Subscript[θr7, 3])=(Subscript[θr8, 3])=0;
(Subscript[Aj3, 3])=(Subscript[Aj4, 3])=(Subscript[Aj5, 3])=(Subscript[Aj6, 3])=(Subscript[Aj7, 3])=(Subscript[Aj8, 3])=0;
(Subscript[Ukj3, 3])=(Subscript[Ukj4, 3])=(Subscript[Ukj5, 3])=(Subscript[Ukj6, 3])=(Subscript[Ukj7, 3])=(Subscript[Ukj8, 3])=0;

(*-----Místnost 1.05-----*)
(*plocha podlahy*)
(Subscript[Ap, 4])=27.51;
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 4])=2.18*2.38;
(Subscript[Ae2, 4])=(6.1+4.08)*2.6-(Subscript[Ae1, 4]);
(Subscript[Ae3, 4])=(Subscript[Ap, 4])-2.612*0.6;
(Subscript[Uke1, 4])=Uokna;
(Subscript[Uke2, 4])=UzedB1;
(Subscript[Uke3, 4])=Ustrecha;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 4])=(Subscript[Aeu2, 4])=0;
(Subscript[Ukeu1, 4])=(Subscript[Ukeu2, 4])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 4])=(Subscript[θ, o1]);
(Subscript[Aj1, 4])=1.2*2.6;
(Subscript[Ukj1, 4])=UzedC1;
(Subscript[θr2, 4])=(Subscript[θr3, 4])=(Subscript[θr4, 4])=(Subscript[θr5, 4])=(Subscript[θr6, 4])=(Subscript[θr7, 4])=(Subscript[θr8, 4])=0;
(Subscript[Aj2, 4])=(Subscript[Aj3, 4])=(Subscript[Aj4, 4])=(Subscript[Aj5, 4])=(Subscript[Aj6, 4])=(Subscript[Aj7, 4])=(Subscript[Aj8, 4])=0;
(Subscript[Ukj2, 4])=(Subscript[Ukj3, 4])=(Subscript[Ukj4, 4])=(Subscript[Ukj5, 4])=(Subscript[Ukj6, 4])=(Subscript[Ukj7, 4])=(Subscript[Ukj8, 4])=0;

(*-----Místnost 1.06-----*)
(*plocha podlahy*)
(Subscript[Ap, 5])=15.39;
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 5])=0.98*2.38;
(Subscript[Ae2, 5])=4.735*2.6-(Subscript[Ae1, 5]);
(Subscript[Ae3, 5])=0;
(Subscript[Uke1, 5])=Uokna;
(Subscript[Uke2, 5])=UzedB1;
(Subscript[Uke3, 5])=0;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 5])=0;
(Subscript[Aeu2, 5])=0;
(Subscript[Ukeu1, 5])=0;
(Subscript[Ukeu2, 5])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 5])=(Subscript[θ, koup]);
(Subscript[Aj1, 5])=(2.208+1.255+0.6-0.125)*2.6;
(Subscript[Ukj1, 5])=UzedC1;
(Subscript[θr2, 5])=(Subscript[θr3, 5])=(Subscript[θr4, 5])=(Subscript[θr5, 5])=(Subscript[θr6, 5])=(Subscript[θr7, 5])=(Subscript[θr8, 5])=0;
(Subscript[Aj2, 5])=(Subscript[Aj3, 5])=(Subscript[Aj4, 5])=(Subscript[Aj5, 5])=(Subscript[Aj6, 5])=(Subscript[Aj7, 5])=(Subscript[Aj8, 5])=0;
(Subscript[Ukj2, 5])=(Subscript[Ukj3, 5])=(Subscript[Ukj4, 5])=(Subscript[Ukj5, 5])=(Subscript[Ukj6, 5])=(Subscript[Ukj7, 5])=(Subscript[Ukj8, 5])=0;

(*-----Místnost 1.07-----*)
(*plocha podlahy*)
(Subscript[Ap, 6])=10.4;
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 6])=(Subscript[Ae2, 6])=(Subscript[Ae3, 6])=0;
(Subscript[Uke1, 6])=(Subscript[Uke2, 6])=(Subscript[Uke3, 6])=0;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 6])=(Subscript[Aeu2, 6])=0;
(Subscript[Ukeu1, 6])=(Subscript[Ukeu2, 6])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 6])=(Subscript[θr2, 6])=(Subscript[θr3, 6])=(Subscript[θr4, 6])=(Subscript[θ, o]);
(Subscript[Aj1, 6])=0.7*1.97;
(Subscript[Aj2, 6])=(3.7+2.81*2)*2.6-(Subscript[Aj1, 6]);
(Subscript[Aj3, 6])=3.7*2.6;
(Subscript[Aj4, 6])=(Subscript[Ap, 6]);
(Subscript[Ukj1, 6])=Udverei;
(Subscript[Ukj2, 6])=UzedC1;
(Subscript[Ukj3, 6])=UzedB2;
(Subscript[Ukj4, 6])=Upodlaha;
(Subscript[θr5, 6])=(Subscript[θr6, 6])=(Subscript[θr7, 6])=(Subscript[θr8, 6])=0;
(Subscript[Aj5, 6])=(Subscript[Aj6, 6])=(Subscript[Aj7, 6])=(Subscript[Aj8, 6])=0;
(Subscript[Ukj5, 6])=(Subscript[Ukj6, 6])=(Subscript[Ukj7, 6])=(Subscript[Ukj8, 6])=0;

(*-----Místnost 1.08-----*)
(*plocha podlahy*)
(Subscript[Ap, 7])=9.74;
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 7])=(Subscript[Ae2, 7])=(Subscript[Ae3, 7])=0;
(Subscript[Uke1, 7])=(Subscript[Uke2, 7])=(Subscript[Uke3, 7])=0;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 7])=1.2*2.6;
(Subscript[Aeu2, 7])=0;
(Subscript[Ukeu1, 7])=UzedB2;
(Subscript[Ukeu2, 7])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 7])=(Subscript[θr2, 7])=(Subscript[θr3, 7])=(Subscript[θr4, 7])=(Subscript[θ, o1]);
(Subscript[θr5, 7])=(Subscript[θr6, 7])=(Subscript[θ, koup]);
(Subscript[θr7, 7])=(Subscript[θr8, 7])=(Subscript[θ, t]);
(Subscript[Aj1, 7])=(Subscript[Aj7, 7])=0.8*1.97;
(Subscript[Aj2, 7])=2.887*2.6-(Subscript[Aj1, 7]);
(Subscript[Aj3, 7])=(Subscript[Aj5, 7])=0.7*1.97;
(Subscript[Aj4, 7])=(3.7+2.81)*2.6-(Subscript[Aj3, 7]);
(Subscript[Aj6, 7])=2.292*2.6-(Subscript[Aj5, 7]);
(Subscript[Aj8, 7])=2.6*2.6-(Subscript[Aj7, 7]);
(Subscript[Ukj1, 7])=(Subscript[Ukj3, 7])=(Subscript[Ukj5, 7])=(Subscript[Ukj7, 7])=Udverei;
(Subscript[Ukj2, 7])=UzedB2;
(Subscript[Ukj4, 7])=(Subscript[Ukj6, 7])=(Subscript[Ukj8, 7])=UzedC1;

(*-----Místnost 1.09-----*)
(*plocha podlahy*)
(Subscript[Ap, 8])=16.6;
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 8])=2*0.6*2.2;
(Subscript[Ae2, 8])=(2.4+6.155)*2.6-(Subscript[Ae1, 8]);
(Subscript[Ae3, 8])=0;
(Subscript[Uke1, 8])=Uokna;
(Subscript[Uke2, 8])=UzedB1;
(Subscript[Uke3, 8])=0;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 8])=(Subscript[Aeu2, 8])=0;
(Subscript[Ukeu1, 8])=(Subscript[Ukeu2, 8])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 8])=(Subscript[θr2, 8])=(Subscript[θr3, 8])=(Subscript[θ, o]);
(Subscript[θr4, 8])=(Subscript[θ, t]);
(Subscript[Aj1, 8])=0.7*1.97;
(Subscript[Aj2, 8])=(2.229+2.887)*2.6-(Subscript[Aj1, 8]);
(Subscript[Aj3, 8])=3*3.7;
(Subscript[Aj4, 8])=3*2.6;
(Subscript[Ukj1, 8])=Udverei;
(Subscript[Ukj2, 8])=(Subscript[Ukj4, 8])=UzedC1;
(Subscript[Ukj3, 8])=Upodlaha;
(Subscript[θr5, 8])=(Subscript[θr6, 8])=(Subscript[θr7, 8])=(Subscript[θr8, 8])=0;
(Subscript[Aj5, 8])=(Subscript[Aj6, 8])=(Subscript[Aj7, 8])=(Subscript[Aj8, 8])=0;
(Subscript[Ukj5, 8])=(Subscript[Ukj56, 8])=(Subscript[Ukj7, 8])=(Subscript[Ukj9, 8])=0;

(*-----Místnost 1.10-----*)
(*plocha podlahy*)
(Subscript[Ap, 9])=7.9;
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 9])=1*0.5;
(Subscript[Ae2, 9])=2.6*2.6-(Subscript[Ae1, 9]);
(Subscript[Ae3, 9])=0;
(Subscript[Uke1, 9])=Uokna;
(Subscript[Uke2, 9])=UzedB1;
(Subscript[Uke3, 9])=0;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 9])=0.8*1.97;
(Subscript[Aeu2, 9])=3*2.6-(Subscript[Aeu1, 9]);
(Subscript[Ukeu1, 9])=Udverei;
(Subscript[Ukeu2, 9])=UzedB2;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 9])=(Subscript[θr2, 9])=(Subscript[θ, o]);
(Subscript[θr3, 9])=(Subscript[θ, koup]);
(Subscript[θr4, 9])=(Subscript[θ, o1]);
(Subscript[Aj1, 9])=0.8*1.97;
(Subscript[Aj2, 9])=(2.229+2.887)*2.6-(Subscript[Aj1, 9]);
(Subscript[Aj3, 9])=3*2.6;
(Subscript[Aj4, 9])=3*2.6;
(Subscript[Ukj1, 9])=Udverei;
(Subscript[Ukj2, 9])=(Subscript[Ukj3, 9])=UzedC1;
(Subscript[Ukj4, 9])=Upodlaha;
(Subscript[θr5, 9])=(Subscript[θr6, 9])=(Subscript[θr7, 9])=(Subscript[θr8, 9])=0;
(Subscript[Aj5, 9])=(Subscript[Aj6, 9])=(Subscript[Aj7, 9])=(Subscript[Aj8, 9])=0;
(Subscript[Ukj5, 9])=(Subscript[Ukj6, 9])=(Subscript[Ukj7, 9])=(Subscript[Ukj8, 9])=0;

(*-----Místnost 1.11 (schodiště 1p. namísto garáže)-----*)
(*plocha podlahy*)
(Subscript[Ap, 10])=4.34;
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 10])=(Subscript[Ae2, 10])=(Subscript[Ae3, 10])=0;
(Subscript[Uke1, 10])=(Subscript[Uke2, 10])=(Subscript[Uke3, 10])=0;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 10])=(Subscript[Aeu2, 10])=0;
(Subscript[Ukeu1, 10])=(Subscript[Ukeu2, 9])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 10])=(Subscript[θr2, 10])=(Subscript[θr3, 10])=(Subscript[θr4, 10])=(Subscript[θr5, 10])=(Subscript[θr6, 10])=(Subscript[θr7, 10])=(Subscript[θr8, 10])=0;
(Subscript[Aj1, 10])=(Subscript[Aj2, 10])=(Subscript[Aj3, 10])=(Subscript[Aj4, 10])=(Subscript[Aj5, 10])=(Subscript[Aj6, 10])=(Subscript[Aj7, 10])=(Subscript[Aj8, 10])=0;
(Subscript[Ukj1, 10])=(Subscript[Ukj2, 10])=(Subscript[Ukj3, 9])=(Subscript[Ukj4, 10])=(Subscript[Ukj5, 10])=(Subscript[Ukj6, 10])=(Subscript[Ukj7, 10])=(Subscript[Ukj8, 10])=0;

(*-----Místnost 2.01-----*)
(*plocha podlahy*)
(Subscript[Ap, 11])=11.24+(Subscript[Ap, 10]);
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 11])=1*0.6;
(Subscript[Ae2, 11])=3.7*2.1+1.1*2.73-(Subscript[Ae1, 11]);
(Subscript[Ae3, 11])=(Subscript[Ap, 11]);
(Subscript[Uke1, 11])=Uokna;
(Subscript[Uke2, 11])=UzedB1;
(Subscript[Uke3, 11])=Ustrecha;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 11])=(Subscript[Aeu2, 11])=0;
(Subscript[Ukeu1, 11])=(Subscript[Ukeu2, 11])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 11])=(Subscript[θr2, 11])=(Subscript[θ, o1]);
(Subscript[θr3, 11])=(Subscript[θr4, 11])=(Subscript[θr5, 11])=(Subscript[θ, koup]);
(Subscript[Aj1, 11])=(Subscript[Aj3, 11])=0.7*1.97;
(Subscript[Aj2, 11])=(Subscript[Aj4, 11])=2.4*(((2.1+2.373)/2))-(Subscript[Aj1, 11]);
(Subscript[Aj5, 11])=3*3.7;
(Subscript[Ukj1, 11])=(Subscript[Ukj3, 11])=Udverei;
(Subscript[Ukj2, 11])=(Subscript[Ukj4, 11])=UzedC1;
(Subscript[Ukj5, 11])=Upodlaha;
(Subscript[θr6, 11])=(Subscript[θr7, 11])=(Subscript[θr8, 11])=0;
(Subscript[Aj6, 11])=(Subscript[Aj7, 11])=(Subscript[Aj8, 11])=0;
(Subscript[Ukj6, 11])=(Subscript[Ukj7, 11])=(Subscript[Ukj8, 11])=0;

(*-----Místnost 2.02-----*)
(*plocha podlahy*)
(Subscript[Ap, 12])=6.18;
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 12])=1.2*0.5;
(Subscript[Ae2, 12])=2.574*2.1+2.4*(((2.1+2.373)/2))-(Subscript[Ae1, 11]);
(Subscript[Ae3, 12])=(Subscript[Ap, 12]);
(Subscript[Uke1, 12])=Uokna;
(Subscript[Uke2, 12])=UzedB1;
(Subscript[Uke3, 12])=Ustrecha;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 12])=(Subscript[Aeu2, 12])=0;
(Subscript[Ukeu1, 12])=(Subscript[Ukeu2, 12])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 12])=(Subscript[θr2, 12])=(Subscript[θ, o]);
(Subscript[θr3, 12])=(Subscript[θ, t]);
(Subscript[Aj1, 12])=0.7*1.97;
(Subscript[Aj2, 12])=2.4*(((2.1+2.373)/2))+2.574*2.373-(Subscript[Aj1, 12]);
(Subscript[Aj3, 12])=2.4*2.574;
(Subscript[Ukj1, 12])=Udverei;
(Subscript[Ukj2, 12])=UzedC1;
(Subscript[Ukj3, 12])=Upodlaha;
(Subscript[θr4, 12])=(Subscript[θr5, 12])=(Subscript[θr6, 12])=(Subscript[θr7, 12])=(Subscript[θr8, 12])=0;
(Subscript[Aj4, 12])=(Subscript[Aj5, 12])=(Subscript[Aj6, 12])=(Subscript[Aj7, 12])=(Subscript[Aj8, 12])=0;
(Subscript[Ukj4, 12])=(Subscript[Ukj5, 12])=(Subscript[Ukj6, 12])=(Subscript[Ukj7, 12])=(Subscript[Ukj8, 12])=0;

(*-----Místnost 2.03-----*)
(*plocha podlahy*)
(Subscript[Ap, 13])=16.79;
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 13])=1.5*1.2+2*1.8*0.6;
(Subscript[Ae2, 13])=3.737*2.73+4.679*(((2.373+2.73)/2))-(Subscript[Ae1, 13]);
(Subscript[Ae3, 13])=(Subscript[Ap, 13]);
(Subscript[Uke1, 13])=Uokna;
(Subscript[Uke2, 13])=UzedB1;
(Subscript[Uke3, 13])=Ustrecha;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 13])=(Subscript[Aeu2, 13])=0;
(Subscript[Ukeu1, 13])=(Subscript[Ukeu2, 13])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 13])=(Subscript[θr2, 13])=(Subscript[θ, o1]);
(Subscript[Aj1, 13])=2.574*2.373;
(Subscript[Aj2, 13])=2.5*2.81;
(Subscript[Ukj1, 13])=UzedC1;
(Subscript[Ukj2, 13])=Upodlaha;
(Subscript[θr3, 13])=(Subscript[θr4, 13])=(Subscript[θr5, 13])=(Subscript[θr6, 13])=(Subscript[θr7, 13])=(Subscript[θr8, 13])=0;
(Subscript[Aj3, 13])=(Subscript[Aj4, 13])=(Subscript[Aj5, 13])=(Subscript[Aj6, 13])=(Subscript[Aj7, 13])=(Subscript[Aj8, 13])=0;
(Subscript[Ukj3, 13])=(Subscript[Ukj4, 13])=(Subscript[Ukj5, 13])=(Subscript[Ukj6, 13])=(Subscript[Ukj7, 13])=(Subscript[Ukj8, 13])=0;

(*-----Místnost 2.04-----*)
(*plocha podlahy*)
(Subscript[Ap, 14])=16.72;
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 14])=0.8*2.2+2*1.8*0.6;
(Subscript[Ae2, 14])=3.738*2.73+4.679*(((2.373+2.73)/2))-(Subscript[Ae1, 14]);
(Subscript[Ae3, 14])=(Subscript[Ap, 14]);
(Subscript[Uke1, 14])=Uokna;
(Subscript[Uke2, 14])=UzedB1;
(Subscript[Uke3, 14])=Ustrecha;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 14])=(Subscript[Aeu2, 14])=0;
(Subscript[Ukeu1, 14])=(Subscript[Ukeu2, 14])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 14])=(Subscript[θ, koup]);
(Subscript[Aj1, 14])=2.3*2.373;
(Subscript[Ukj1, 14])=UzedC1;
(Subscript[θr2, 14])=(Subscript[θr3, 14])=(Subscript[θr4, 14])=(Subscript[θr5, 14])=(Subscript[θr6, 14])=(Subscript[θr7, 14])=(Subscript[θr8, 14])=0;
(Subscript[Aj2, 14])=(Subscript[Aj3, 14])=(Subscript[Aj4, 14])=(Subscript[Aj5, 14])=(Subscript[Aj6, 14])=(Subscript[Aj7, 14])=(Subscript[Aj8, 14])=0;
(Subscript[Ukj2, 14])=(Subscript[Ukj3, 14])=(Subscript[Ukj4, 14])=(Subscript[Ukj5, 14])=(Subscript[Ukj6, 14])=(Subscript[Ukj7, 14])=(Subscript[Ukj8, 14])=0;

(*-----Místnost 2.05-----*)
(*plocha podlahy*)
(Subscript[Ap, 15])=5.52;
(*konstrukce s sousedící s exterierém*)
(Subscript[Ae1, 15])=1.2*0.5;
(Subscript[Ae2, 15])=2.3*2.1+2.4*(((2.1+2.373)/2))-(Subscript[Ae1, 15]);
(Subscript[Ae3, 15])=(Subscript[Ap, 15]);
(Subscript[Uke1, 15])=Uokna;
(Subscript[Uke2, 15])=UzedB1;
(Subscript[Uke3, 15])=Ustrecha;
(*konstrukce s sousedící s nevytápěným prostorem*)
(Subscript[Aeu1, 15])=(Subscript[Aeu2, 15])=0;
(Subscript[Ukeu1, 15])=(Subscript[Ukeu2, 15])=0;
(*konstrukce s sousedící s prostorem s výrazně odlišnou teplotou*)
(Subscript[θr1, 15])=(Subscript[θr2, 15])=(Subscript[θ, o]);
(Subscript[Aj1, 15])=0.7*1.97;
(Subscript[Aj2, 15])=2.3*2.373+2.4*(((2.1+2.373)/2))-(Subscript[Aj1, 15]);
(Subscript[Ukj1, 15])=Udverei;
(Subscript[Ukj2, 15])=UzedC1;
(Subscript[θr3, 15])=(Subscript[θr4, 15])=(Subscript[θr5, 15])=(Subscript[θr6, 15])=(Subscript[θr7, 15])=(Subscript[θr8, 15])=0;
(Subscript[Aj3, 15])=(Subscript[Aj4, 15])=(Subscript[Aj5, 15])=(Subscript[Aj6, 15])=(Subscript[Aj7, 15])=(Subscript[Aj8, 15])=0;
(Subscript[Ukj3, 15])=(Subscript[Ukj4, 15])=(Subscript[Ukj5, 15])=(Subscript[Ukj6, 15])=(Subscript[Ukj7, 15])=(Subscript[Ukj8, 15])=0;

(*------------------------VÝPOČET------------------------*)
For[x=1; pom1=1,x<16,x++, pom1++;
(*nastavení vypisování textu*)
If[x==1,Print[Column[{"************************VÝSLEDKY************************","--------------------ZTRÁTY MÍSTNOSTÍ--------------------",{"pro teplotu (Subscript[θ, e]) = ", (Subscript[θ, e])[[1]], "°C"},"--------První patro--------"}]]];
If[x==11,Print["--------Druhé patro--------"]];
If[pom1> 11,Print["místnost ", pom1-11],
If[pom1==2,Print["místnost 1+2"],Print["místnost ", pom1 ]]];
(*teploty*)
If[x==8||x==15,(Subscript[θ, x])=(Subscript[θ, koup]),(Subscript[θ, x])=(Subscript[θ, o])];
If[x==1||x==6||x==12,(Subscript[θ, x])=(Subscript[θ, o1])];
If[x==9,(Subscript[θ, x])=(Subscript[θ, t])];
(*měrná ztráta do exteriéru*)
(Subscript[A, 1])=(Subscript[Ae1, x]);(Subscript[A, 2])=(Subscript[Ae2, x]);(Subscript[A, 3])=(Subscript[Ae3, x]);
(Subscript[Uk, 1])=(Subscript[Uke1, x]);(Subscript[Uk, 2])=(Subscript[Uke2, x]);(Subscript[Uk, 3])=(Subscript[Uke3, x]);
(Subscript[H, Te])=(Subscript[e, k])*((∑_(i=1)^3)(Subscript[A, i])*((Subscript[Uk, i])+(Subscript[ΔU, tb])));
(*měrná ztráta přes nevytápěnou místnost do exteriéru*)
bu= (((Subscript[θ, x])-(Subscript[θ, u]))/((Subscript[θ, x])-(Subscript[θ, e])));
If[(Subscript[θ, u])≤(Subscript[θ, e]),bu=1];
(Subscript[A, 1])=(Subscript[Aeu1, x]);(Subscript[A, 2])=(Subscript[Aeu2, x]);
(Subscript[Uk, 1])=(Subscript[Ukeu1, x]);(Subscript[Uk, 2])=(Subscript[Ukeu2, x]);
(Subscript[H, Tue])=bu*((∑_(i=1)^2)(Subscript[A, i])*(Subscript[Uk, i]));
(*měrná ztráta přes podlahu do zeminy*)
If[x>10, (Subscript[f, g1])=0];
(Subscript[f, g2])= (((Subscript[θ, x])-(Subscript[θ, em]))/((Subscript[θ, x])-(Subscript[θ, e])));
(Subscript[H, Tg])=(Subscript[f, g1])*(Subscript[f, g2])*(Subscript[U, ekvi])*(Subscript[Ap, x])*(Subscript[G, w]);
(*měrná ztráta větráním*)
If[x==4||x==8||x==15,(Subscript[n, x])=(Subscript[n, k]),(Subscript[n, x])=(Subscript[n, o])];
If[x==3, (Subscript[V, x])=(Subscript[Ap, x])*3.65,(Subscript[V, x])=(Subscript[Ap, x])*2.6];
If[x==10, (Subscript[V, x])=(Subscript[Ap, x])*(2.6/2)];
If[x>10, (Subscript[V, x])=(Subscript[Ap, x])*((2.73+2.1)/2)];
(Subscript[V, i])=(Subscript[V, x])*(Subscript[n, x]);
(Subscript[H, V])=(Subscript[V, i])*k;
(*měrná ztráta do místnosti s výrazně odlišnou teplotou*)
(Subscript[bj, 1])= (((Subscript[θ, x])-(Subscript[θr1, x]))/((Subscript[θ, x])-(Subscript[θ, e]))); (Subscript[bj, 2])= (((Subscript[θ, x])-(Subscript[θr2, x]))/((Subscript[θ, x])-(Subscript[θ, e]))); (Subscript[bj, 3])= (((Subscript[θ, x])-(Subscript[θr3, x]))/((Subscript[θ, x])-(Subscript[θ, e]))); (Subscript[bj, 4])= (((Subscript[θ, x])-(Subscript[θr4, x]))/((Subscript[θ, x])-(Subscript[θ, e]))); (Subscript[bj, 5])= (((Subscript[θ, x])-(Subscript[θr5, x]))/((Subscript[θ, x])-(Subscript[θ, e]))); (Subscript[bj, 6])= (((Subscript[θ, x])-(Subscript[θr6, x]))/((Subscript[θ, x])-(Subscript[θ, e]))); (Subscript[bj, 7])= (((Subscript[θ, x])-(Subscript[θr7, x]))/((Subscript[θ, x])-(Subscript[θ, e]))); (Subscript[bj, 8])= (((Subscript[θ, x])-(Subscript[θr8, x]))/((Subscript[θ, x])-(Subscript[θ, e])));
(Subscript[A, 1])=(Subscript[Aj1, x]);(Subscript[A, 2])=(Subscript[Aj2, x]); (Subscript[A, 3])=(Subscript[Aj3, x]); (Subscript[A, 4])=(Subscript[Aj4, x]); (Subscript[A, 5])=(Subscript[Aj5, x]); (Subscript[A, 6])=(Subscript[Aj6, x]); (Subscript[A, 7])=(Subscript[Aj7, x]); (Subscript[A, 8])=(Subscript[Aj8, x]);
(Subscript[Uk, 1])=(Subscript[Ukj1, x]);(Subscript[Uk, 2])=(Subscript[Ukj2, x]); (Subscript[Uk, 3])=(Subscript[Ukj3, x]); (Subscript[Uk, 4])=(Subscript[Ukj4, x]); (Subscript[Uk, 5])=(Subscript[Ukj5, x]);
(Subscript[Uk, 6])=(Subscript[Ukj6, x]); (Subscript[Uk, 7])=(Subscript[Ukj7, x]); (Subscript[Uk, 8])=(Subscript[Ukj8, x]);

(Subscript[H, Tj])=((∑_(i=1)^8)(Subscript[A, i])*(Subscript[Uk, i])*(Subscript[bj, i]));
(*Ztráty místnosti*)
(Subscript[Φt, x])=((Subscript[H, Te])+(Subscript[H, Tue])+(Subscript[H, Tg])+(Subscript[H, Tj]))*((Subscript[θ, x])-(Subscript[θ, e]));
(Subscript[Φv, x])=(Subscript[H, V])*((Subscript[θ, x])-(Subscript[θ, e])); 
(Subscript[Φ, x])=(Subscript[Φt, x])+(Subscript[Φv, x]); 
Print["θ = ",(Subscript[θ, x])," °C, Φt = ",(Subscript[Φt, x])[[1]]," W, Φv = ",(Subscript[Φv, x])[[1]]," W, Φ = ",(Subscript[Φ, x])[[1]]," W"];
]
(*celkový ztráty pater*)
Φ1t=(∑_(i=1)^10)(Subscript[Φt, i]); Φ1v=(∑_(i=1)^10)(Subscript[Φv, i]);Φ1p=(∑_(i=1)^10)(Subscript[Φ, i]);
Φ2t=(∑_(i=11)^15)(Subscript[Φt, i]); Φ2v=(∑_(i=11)^15)(Subscript[Φv, i]); Φ2p=(∑_(i=11)^15)(Subscript[Φ, i]);
Print["--------Celkem první patro--------"];Print["Φt = ",Φ1t[[1]]," W, Φv = ",Φ1v[[1]]," W, Φ = ",Φ1p[[1]]," W"];Print["--------Celkem druhé patro--------"];
Print["Φt = ",Φ2t[[1]]," W, Φv = ",Φ2v[[1]]," W, Φ = ",Φ2p[[1]]," W"];

(*CELKOVÉ ZTRÁTY*)
Print[Column[{"","","--------------------ZTRÁTY OBJEKTU--------------------",{"pro teplotu (Subscript[θ, e]) = ", (Subscript[θ, e]), "°C"}}]];
(Subscript[Φ, c])=Φ1p+Φ2p;Print[(Subscript[Φ, c])"W"];
Print[Column[{"","","*********ZTRÁTY POKRYTÉ RADIÁTORY********"}]];
(Subscript[Φ, KR])=Φ2p+((Subscript[Φ, 1])+(Subscript[Φ, 5]));Print[(Subscript[Φ, KR]) "W"];
Print[Column[{"","","*********ZTRÁTY POKRYTÉ PODLAHOVÝM VYTÁPĚNÍM********"}]];
(Subscript[Φ, KP])=Φ1p-((Subscript[Φ, 1])+(Subscript[Φ, 5]));Print[(Subscript[Φ, KP]) "W"];
(*interpolace a uložení výsledku*)
dvojice[i_]:={(Subscript[θ, e])[[i]],(Subscript[Φ, c])[[i]]};
dataC=Table[dvojice[i],{i,Length[(Subscript[θ, e])]}];
dvojice[i_]:={(Subscript[θ, e])[[i]],(Subscript[Φ, KR])[[i]]};
dataR=Table[dvojice[i],{i,Length[(Subscript[θ, e])]}];
dvojice[i_]:={(Subscript[θ, e])[[i]],(Subscript[Φ, KP])[[i]]};
dataP=Table[dvojice[i],{i,Length[(Subscript[θ, e])]}];
PtopnyP[Tvzduch_]=Interpolation[dataP,InterpolationOrder->1][Tvzduch];
PtopnyR[Tvzduch_]=Interpolation[dataR,InterpolationOrder->1][Tvzduch];
PtopnyC[Tvzduch_]=Interpolation[dataC,InterpolationOrder->1][Tvzduch];
Put[{PtopnyP[Tvzduch],PtopnyR[Tvzduch],PtopnyC[Tvzduch]},"Ptop"]

Plot[{PtopnyP[Tvzduch],PtopnyR[Tvzduch],PtopnyC[Tvzduch]},{Tvzduch,-15,26},AxesLabel->{"(Subscript[θ, e]) [°C]","Φ [W]"},PlotStyle->{{Thickness[0.008],RGBColor[0,0,1]},{Thickness[0.008],RGBColor[0.6,0.4,0.0]},{Thickness[0.008],RGBColor[1,0,0]}},PlotLegends->{"(Subscript[Φ, podl.])","(Subscript[Φ, rad.])","(Subscript[Φ, c])"}]