SetDirectory["FileName"/.NotebookInformation[EvaluationNotebook[]]/.FrontEnd`FileName[d_List,nam_,___]:>ToFileName[d]];

cisilka={M->1,Icl->1,W->0,RH->60};
Off[FindRoot::"nlnum",Union::"heads",ReplaceAll::"reps"];
PMV1=(0.303*Exp[-2.1*M]+0.028)*(58.15*(M-W)-3.05*(10^(-3))*(5733-406.7*(M-W)-pa)-24.42*((M-W)-1)-(10^(-3))*M*(5867-pa)-0.0814*M*(34-1.0*ta)-3.96*(10^(-8))*fcl*((((tcl+273))^4)-(((tr+273))^4))-fcl*hc*(tcl-1.0*ta));
dosazeni1={hc->2.38*(((tcl-1.0*ta))^0.25),fcl->1.05+0.1*Icl,pa->((RH/100))*Exp[((1500.3+23.5*ta)/(234+ta))]}/.cisilka;
PMV2=PMV1/.dosazeni1/.cisilka;
pravastrana=35.7-1.6282*(M-W)-0.155*Icl*(3.96*(10^(-8))*fcl*((((tcl+273))^4)-(((tr+273))^4))+fcl*hc*(tcl-ta))-tcl;
prav[ta_,tr_]=pravastrana/.dosazeni1/.cisilka;
dosadzateploty[taa_,trr_]:=Union[{ta->taa,tr->trr},FindRoot[prav[taa,trr]==0,{tcl,taa,35}]];
PeEmVe[taa_,trr_]:=PMV2/.dosadzateploty[taa,trr];
najditrr[ta_,pmvvv_]:=trr/.FindRoot[Evaluate[PeEmVe[ta,trr]]==pmvvv,{trr,15,30}]

Plot[najditrr[ta,0],{ta,16,27},AxesLabel->{"ta[°C]","tr[°C]"}]

datapohoda=Table[{ta,najditrr[ta,0]},{ta,16,26,0.5}];
datapohoda>>datpohoda

nasaepmv=PeEmVe[16.2,15.5]

f[pmv_]:=100-95*Exp[-0.03353*(pmv^4)-0.2179*(pmv^2)];
Plot[f[pmv],{pmv,-3,3}]

f[nasaepmv]

Plot[Exp[((1500.3+23.5*ta)/(234+ta))],{ta,0,100}]

mvzduch=(349/282.)*5*7*4

mvzduch*(14-7.5)

Pchlad=((mvzduch*(72-27.5))/(3*3600.))
Pohrev=((mvzduch*(37-27.5))/(3*3600.))

750/3600.

(750000/3600.)/(11.6-4.1)

h=4000.;
g=9.81;
m=70;
energy=m*g*h
energyInOneBeer=0.5*1850*1000.;
(energy/energyInOneBeer)