%% r        - polomer vaku v [cm]
%% h        - vyšška vaku v [cm]
%% V_ml     - objem krvi v [ml] vo vaku
%% Ro_krv	- hustota krvi [kg/m3], od 1048-1066

%% Typ plastu:
%%              1) Polyvinylchlorid (PVC)
%%              2) Polypropylén (PP)
%%              3) Polyetylén (PE)
%%          	4) Etylenvinyl acetát (EVA)

%% T_krv_p - pociatocna teplota krvi
%% T_krv_k - teplota na ktorú chcete krv ochladit
%% T_chlad - teplota chladnicky

function [Xaxes, Yaxes] = chladenie(r, h, V_ml, Ro_krv, typ_plast, T_krv_p, T_krv_k, T_chlad )
   clc;
   lam_PVC=0.20;
   lam_PP=0.22;
   lam_PE=0.35;
   lam_EVA=0.34;
   c_krv = 3617; %% súčiniteľ tepelnej vodivosti krvi
   x=0.025; %% hrúbka plastu
   %% volanie funkcie
   lam_plast = get_lam_plastu(typ_plast);
   %% premeny jednotiek na zakladne
   r_vak = r*0.01;
   h_vak = h*0.01;
   T_krv_p_k=T_krv_p+273.15;
   T_krv_k_k=T_krv_k+273.15;
   T_chlad_k=T_chlad+273.15;
   V=V_ml*0.000001;

   C1=Ro_krv*V*c_krv; %% hmotnostna tepelna kapacita krvi
   S=(r_vak*h_vak)+(r_vak+h_vak); %% plocha vaku
   T_roz=(T_krv_p_k-T_chlad_k);
   Q1=C1*(T_krv_p_k-T_chlad_k); %% teplo, ktore ma krv odovzdat
   D1=Q1*x/lam_plast*S;
   D_PVC=Q1*x/lam_PVC*S;
   D_PP=Q1*x/lam_PP*S;
   D_PE=Q1*x/lam_PE*S;
   D_EVA=Q1*x/lam_EVA*S;
   Q_krv_k=C1*T_krv_k_k;


   X1 = [];
   Y1 = [];
   X2 = [];
   Y2 = [];
   X3 = [];
   Y3 = [];

   T_zmena1=T_krv_p_k;
   T_zmena2=T_krv_p_k;
          %%abc = T_zmena-273.15;

          %%X=[X;0];
          %%Y=[Y;abc];
for i =1:T_roz
          t1=D_PVC*(1/(T_zmena1-T_chlad_k)); %% vypocet casu
          t1=t1-(D1*(1/(T_roz)));
          t2=D_PP*(1/(T_zmena1-T_chlad_k)); %% vypocet casu
          t2=t1-(D1*(1/(T_roz)));
          t3=D_PE*(1/(T_zmena1-T_chlad_k)); %% vypocet casu
          t3=t1-(D1*(1/(T_roz)));
          t4=D_EVA*(1/(T_zmena1-T_chlad_k)); %% vypocet casu
          t4=t1-(D1*(1/(T_roz)));
          T_vys1=T_zmena1-273.15; %% ak dosiahne teplotu T_krv_k=T_vys1, bod vyznacit
          T_zmena1=T_zmena1-1;
          X1=[X1;t1]; %% vyskreslovať všetko do jedneho grafu
          Y1=[Y1;T_vys1]; %%
          X2=[X1;t2];%%
          Y2=[Y1;T_vys1];%%
          X3=[X1;t3];%%
          Y3=[Y1;T_vys1]; %%
          X4=[X1;t4];%%
          Y4=[Y1;T_vys1]; %%


   end

   for i =1:T_roz
          t2=D1*(1/(T_zmena2-T_krv_k_k)); %% vypocet casu
          t2=t2-(D1*(1/(T_roz)));
          Q_krv=C1*T_zmena2; %%porovnavat s Q_krv_k ak je Q_krv=Q_krv_k tento bod vyznacit
          T_vys2=T_zmena2-273.15; %% ak dosiahne teplotu T_krv_k=T_vys2, bod vyznacit
          T_zmena2=T_zmena2-1;

          X2=[X2;t2];
          Y2=[Y2;T_vys2];
          X3=[X3;t2];
          Y3=[Y3;Q_krv];

   end
 Xaxes = X1;
   Yaxes = Y1;
   figure;
   plot(X1, Y1);

   Xaxes = X2;
   Yaxes = Y2;
   figure;
   plot(X2, Y2);

   Xaxes = X3;
   Yaxes = Y3;
   figure;
   plot(X3, Y3);


end
%%volana funkcia, sucinitel prestupu tepla
function lam_plastu = get_lam_plastu(typ_plast)
   switch logical(true)
       case (typ_plast==1 ),  lam_plastu = 0.20;
       case (typ_plast==2 ),  lam_plastu = 0.22;
       case (typ_plast==3 ),  lam_plastu = 0.35;
       case (typ_plast==4 ),  lam_plastu = 0.34;
       otherwise,  lam_plastu = 0.20;
   end;
end