clc;
clear all;

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

Kcw=qN/((TwN-TzN)+0.25*(TwN-TpN));
Kcp=(0.25*qN/((TwN-TzN)+0.25*(TwN-TpN)))*((TwN-TpN)/(TpN-TzN));
Kwcp=0.25*Kcw;

kolory='rgbcm';
figure;

%rysowanie wykresu Twew(q)
q=[500:1:2000];
i=1;
for Tz=-20 : 10 : 20;
    Tw=(q/Kcw)+Tz-((Kcp*q)/(Kcw*(Kcw+5*Kcp)));
    subplot(2,2,1);
    hold on;
grid on;
    plot(q, Tw, kolory(mod(i,length(kolory))+1));
    title('Wykres Twew(q)');
    xlabel('q - Moc grzalki [W]');
    ylabel('Twew - temperatura wnetrza [*C]');

    i=i+1;
end
legend('Tzew= -20', 'Tzew= -10', 'Tzew=0', 'Tzew= 10', 'Tzew= 20');

%rysowanie wykresu Tp(q)
q=[500:1:2000];
i=1;
for Tz=-20 : 10 : 20;
    Tp=(q+Kcw*Tz+5*Kcp*Tz)/(Kcw+5*Kcp);
    subplot(2,2,2);
hold on;
grid on;
    plot(q, Tp,(kolory(mod(i,length(kolory))+1)));
  title('Wykres Tp(q)');
      xlabel('q - Moc grzalki [W]');
    ylabel('Tp - temperatura poddasza [*C]');
      i=i+1;
end
legend('Tzew= -20', 'Tzew= -10', 'Tzew=0', 'Tzew= 10', 'Tzew= 20');

%rysowanie wykresu Twew(Tzew)
Tz=[-20:1:20];
i=1;
for q=500:500:2000;
    Tw=(q/Kcw)+Tz-((Kcp*q)/(Kcw*(Kcw+5*Kcp)));
     subplot(2,2,3);
     hold on;
    grid on;
    plot(Tz, Tw,(kolory(mod(i,length(kolory))+1)));
       title('Wykres Twew(Tzew)');
       xlabel('Tzew - temperatura zewnetrzna [*C]');
    ylabel('Twew - temperatura wnetrza[*C]');
      i=i+1;
end
legend('q= 500', 'q=1000', 'q=1500', 'q=2000');

%rysowanie wykresu Tp(Tzew)
Tz=[-20:1:20];
i=1;
for q=500:500:2000;
   Tp=(q+Kcw*Tz+5*Kcp*Tz)/(Kcw+5*Kcp);
   subplot(2,2,4);
   hold on;
grid on;
    plot(Tz, Tp, (kolory(mod(i,length(kolory))+1)));
       title('Wykres Tp(Tzew)');
              xlabel('Tzew - temperatura zewnetrzna [*C]');
    ylabel('Twew - temperatura poddasza[*C]');
      i=i+1;
end
legend('q= 500', 'q=1000', 'q=1500', 'q=2000');