clear all


QgN=1000;
TzN=-20;
TwN=20;
TpN=10;

Kcw=QgN/(1.25*TwN-TzN-0.25*TpN);
Kcwp=0.25*Kcw;
Kcp=Kcwp*(TwN-TpN)/(TpN-TzN);

Qg=500:250:2500;
Tzew=-20:5:20;
kolor='rgcbrgcbrgcbrgcb'


for i=1:length(Qg);
    figure(1);
    grid on;
    hold on;

    Tp=Qg(i)*(1/(5*Kcp+Kcw))+Tzew;
    Twew=(Qg(i)+Kcw*Tzew+Kcwp*Tp)/(Kcw+Kcwp)

    plot(Tzew, Tp, kolor(i));
    xlabel('Tzew');
    ylabel('Tp');
    legend('500W','750W','1000W','1250W','1500W','1750W', '2000W', '2250W', '2500W');

    figure(2);
    grid on;
    hold on;
    plot(Tzew, Twew,kolor(i));
    xlabel('Tzew');
    ylabel('Twew');
    legend('500W','750W','1000W','1250W','1500W','1750W', '2000W', '2250W', '2500W');
end

for i=1:length(Tzew);
    figure(3);
    grid on;
    hold on;

    Tp=Qg*(1/(5*Kcp+Kcw))+Tzew(i);
    Twew=(Qg+Kcw*Tzew(i)+Kcwp*Tp)/(Kcw+Kcwp)

    plot(Qg, Tp, kolor(i));
    xlabel('Qg');
    ylabel('Tp');
    legend('-20*C','-15*C','-10*C','-5*C','0*C','5*C', '10*C', '15*C', '20*C');

    figure(4);
    grid on;
    hold on;
    plot(Qg, Twew, kolor(i));
    xlabel('Qg');
    ylabel('Twew');
    legend('-20*C','-15*C','-10*C','-5*C','0*C','5*C', '10*C', '15*C', '20*C');
end