ClearAll["Global`*"];
$HistoryLength=2;
Off[InterpolatingFunction::dmval];
SetDirectory[NotebookDirectory[]];

(*známé energetické parametry z tabulky*)
T={15,7,2,-7,-15};
n=Length[T];
Ptop35={8.12,8.12,5.68,5.17,3.92};
Pel35={1.93,1.93,1.7,1.84,1.69};
Ptop45={7.72,7.72,6.64,5.12,3.83};
Pel45={2.3,2.3,2.24,2.35,2.08};
Ptop55={7.1,7.1,7.1*(6.64/7.72),4.25,4.2*(3.83/5.12)};
Pel55={2.73,2.73,2.73*(2.3/2.54),2.44,2.4*(2.35/2.08)};
k=1000; (*kW->W*)
(*vytvoření listu dat s hodnotami*)
radekEl35[i_]:={T[[i]],35.,Pel35[[i]]};
radekEl45[i_]:={T[[i]],45.,Pel45[[i]]};
radekEl55[i_]:={T[[i]],55.,Pel55[[i]]};
radekTop35[i_]:={T[[i]],35.,Ptop35[[i]]};
radekTop45[i_]:={T[[i]],45.,Ptop45[[i]]};
radekTop55[i_]:={T[[i]],55.,Ptop55[[i]]};
dataEl=Flatten[Table[{radekEl35[i],radekEl45[i],radekEl55[i]},{i,n}],1]; 
dataTop=Flatten[Table[{radekTop35[i],radekTop45[i],radekTop55[i]},{i,n}],1];
(*vyjádření průběhů zavislosti energetických hodnot na Tvzduch, Tvoda*)
(*Elektriký příkon*)
Pel[Tvzduch_,Tvoda_]=k*Interpolation[dataEl,InterpolationOrder->{1,1}][Tvzduch,Tvoda];
(*Topný výkon*)
Ptop[Tvzduch_,Tvoda_]=k*Interpolation[dataTop,InterpolationOrder->{1,1}][Tvzduch,Tvoda];
(*Topný faktor*)
eps[Tvzduch_,Tvoda_]=((Ptop[Tvzduch,Tvoda])/(Pel[Tvzduch,Tvoda]));
(*uložení průběhů*)
Put[{Pel[Tvzduch,Tvoda],Ptop[Tvzduch,Tvoda],eps[Tvzduch,Tvoda]},"NIBE AMS10-8"]

(*grafy*)
ts=0.008;
ps={{Thickness[ts],RGBColor[0,0,1]},{Thickness[ts],RGBColor[0.5,0.5,0]},{Thickness[ts],RGBColor[1,0,0]},{GrayLevel[0]}};
Plot[{Ptop[Tvzduch,35],Ptop[Tvzduch,45],Ptop[Tvzduch,50],0},{Tvzduch,-15,15},AxesLabel->{"Tv [°C]","Ptop [W]"},PlotStyle->ps,PlotLegends->{"35 °C","45 °C","55 °C"}]
Plot[{Pel[Tvzduch,35],Pel[Tvzduch,40],Pel[Tvzduch,50],0},{Tvzduch,-15,15},AxesLabel->{"Tv [°C]","Pel [W]"},PlotStyle->ps,PlotLegends->{"35 °C","45 °C","55 °C"}]
Plot[{eps[Tvzduch,35],eps[Tvzduch,40],eps[Tvzduch,50],0},{Tvzduch,-15,15},AxesLabel->{"Tv [°C]","ϵ [-]"},PlotStyle->ps,PlotLegends->{"35 °C","45 °C","55 °C"}]