%% Wykres
P = [10:-1:1];
% i1= XSteam('h_pT', 10, 500);


for k=1:10
if k==1
i(k) = XSteam('h_pT', P(k), 500);
s(k) = XSteam('s_pT', P(k), 500);
else
iact=i(k-1);
sact=s(k-1);
i(k)= XSteam('h_ps', P(k), sact);
i(k) = iact-0.8*(iact-i(k));

s(k) = XSteam('s_ph', P(k), i(k));

end
end

plot(s,i,'-*r');grid on;
xlim([7.5 8.5]);
ylim([2000 4000]);