%I (obiekt)
TwN = 20;                 %temperatura wewnetrzna
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)
Vk = 1 * 1 * 1;           %objętość kotła
QN = 5000;                %cieplo oddawane przez grzejnik
QdN = 0;
PkN = 3 * QN;             %moc kotła

T0=1000; %opóźnienie
Tdelay = 1000;

cvg=Cpw*row*Vw;           %pojemność pieplna grzejnika
cvw=Cpp*rop*Vp;           %pojemność pieplna pokoju
cvk=Cpw*row*Vk;           %pojemność pieplna kotła

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

%II (warunki początkowe)
Tzew0=TzewN;
Tgz0=TgzN;
Tkz0=TkzN;
Tkp0=TkpN;

F0 = FN;
Qd0=QdN;
Pk0 = PkN;

%III (zmiany)
dTzew0=0;
dTgz0=0;

dF1=0;
dF2=0;
dF3=0;

dQd1=0;
dQd2=500;
dQd3=0;

dPk0=0;

Tgsr0 = ( 2 * Cpw*row*F0*Tgz0*(kg + kstr) + kg*(kstr*Tzew0 + Qd0) ) / ( kg * kstr + 2 * Cpw * row * F0 * (kg + kstr) );
Tw0 = (Qd0 + kg * Tgsr0 + kstr * Tzew0) / (kg + kstr);
