%I (obiekt)
TwN = 20;                 %temperatura wewnetrzna
Twz = 20;                 %temperatura wewnętrzna-zadana
TzewN = -20;              %temperatura zewnetrzna
TgzN = 90;                %TgzN/TgpN = 90/70
TgpN = 70;      
Cpp = 1000;               %cieplo wlasciwe powietrza
Cpw = 4175;               %cieplo wlasciwe wody
rop = 1.2;                %gestosc powietrza
row = 1000;               %gestosc wody
Vp = 3 * 3 * 3;           %objetosc pomieszczenia(powietrza)
Vw = 1 * 1 * 0.1;         %objetosc grzejnika(wody)
QN = 5000;                %cieplo oddawane przez piec
QdN = 0;

T0=100; %opóźnienie

cvg=Cpw*row*Vw;
cvw=Cpp*rop*Vp;

TgsrN = (TgzN + TgpN) / 2;
TgpN = 2 * TgsrN - TgzN;
FN = QN / (2 * Cpw * row * (TgzN - TgsrN));
kg = QN / (TgsrN - TwN);
kstr = (QN + QdN) / (TwN - TzewN);

%II (warunki początkowe)
Tzew0=TzewN;
Tgz0=TgzN;
F0=FN;
Qd0=QdN;

%III (zmiany)
dTzew0=0;
dTgz0=0;
dF0=0;
dQd0=200;

Tw0 = Twz;
Tgsr0 = ( Twz*(kg + kstr) - kstr*Tzew0 + Qd0) / kg;

