%% Dane
Pgen = 20000;
m = 25;
P1 = 100;
t1 = 500;
P2 = 6;
P3 = 0.07;
Wu = 10;
mi = 0.5;
etha_w1 = 0.8;
etha_w2 = 0.85;
etha_m = 0.98;
etha_g = 0.99;
etha_k = 0.8;
%% Entlpie
i1 = XSteam('h_pT',P1, t1);
s1 = XSteam('s_pT', P1, t1);
i2t = XSteam('h_ps', P2, s1);
i2 = i1-etha_w1*(i1-i2t);
s2 = XSteam('s_ph', P2, i2);
i3t = XSteam('h_ps', P3, s2);
i3 = i2-etha_w2*(i2-i3t);
i4 = XSteam('hL_p', P3);
i5 = XSteam('hL_p', P2);
%% Wyznaczenie mu i mk
mk = (Pgen/(etha_m*etha_g*(i2-i3)))-m*((i1-i2)/(i2-i3));
mu = m-mk;
%% Moc oddana do sieci cieplnej
Qc = mu*(i2-i5);
%% Bilans energetyczny wymiennika W
i6 = (mk*i4+mu*i5)/m;
%% Zapotrzebowanie na paliwo
B = (m*(i1-i6))/(etha_k*Wu);
%% Sprawność termodynamiczna
%% Wspójczynnik skojarzenia
sigma = Pgen/Qc;