% s-funkcia pre vytvorenie dynamickeho MM rurkoveho vym. tepla
% pri jeho rozdeleni na 4 casti( dz=1/4 )
% vstupna velicina je teplota na vstupe do vymennika t0
% vystupna velicina je teplota na konci vymennika
% je to bezrozmerovy model (cela dlzka =1)
function[sys,x0]=rvtsf(t,x,u,flag)
if flag==0,
sys=[4,0,1,1,0,0];
%x0=[339.0965  335.4977  332.1799  329.1211];
x0=[338.3104  334.0606  330.2095  326.7195];
elseif abs(flag)==1,
dz=5;
q=.002;ro=998;cp=4.182;al=1.1;%al=.9;
d=.05;z=20;tch=293;t0=343;
tau=d*ro*cp/(4*al);v=4*q/(pi*d^2);
k1=v/dz;k2=k1+1/tau;
% volba vstupnej veliciny
t0=u;
sys(1)=-k2*x(1)+tch/tau+k1*t0;
sys(2)=-k2*x(2)+tch/tau+k1*x(1);
sys(3)=-k2*x(3)+tch/tau+k1*x(2);
sys(4)=-k2*x(4)+tch/tau+k1*x(3);
sys=sys;
elseif flag==3,
   sys=x(4);
else
sys=[];
end