clear all
TzewN = -20;
TwewN = 20;
TpN = 10;
QgN = 1000;
kolor = 'gbcrmk'

Kcw = QgN/(1.25*TwewN - 0.25*TpN - TzewN);
Kcwp = 0.25*Kcw;
Kcp = Kcwp*(TwewN - TpN)/(TpN - TzewN);

Tzew = [TzewN:10:30];
Qg = [0:200:QgN];

figure(1)
hold on;
grid on;
title('Charakterystyka Twew w zaleznosci od Qg','FontSize',12);
xlabel('Moc grzalki [W]','FontSize',16);
ylabel('Temperatura wewnetrzna [°C]','FontSize',16);
figure(2)
hold on;
grid on;
title('Charakterystyka Tp w zaleznosci od Qg','FontSize',12);
xlabel('Moc grzalki [W]','FontSize',16);
ylabel('Temperatura poddasza [°C]','FontSize',16);
figure(3)
hold on;
grid on;
title('Charakterystyka Twew w zaleznosci od Tzew','FontSize',12);
xlabel('Temperatura zewnetrzna [°C]','FontSize',16);
ylabel('Temperatura wewnetrzna [°C]','FontSize',16);
figure(4)
hold on;
grid on;
title('Charakterystyka Tp w zaleznosci od Tzew','FontSize',12);
xlabel('Temperatura zewnetrzna [°C]','FontSize',16);
ylabel('Temperatura poddasza [°C]','FontSize',16);

%Sprawdzenie poprawnosci wzorow  10 i 20
Tp1 = (QgN + Kcw * (TzewN) + 5 *(Kcp*(TzewN)))/(5*Kcp + Kcw)
Twew = (QgN + Kcw * (TzewN) + 0.25*Kcw*Tp1)/(1.25*Kcw)

%Rysowanie wykresow
for i = 1:length(Qg)
  Tp = (Qg(i) + Kcw.*Tzew+5*(Kcp.*Tzew))/(5*Kcp + Kcw);
  Twew = (Qg(i) + Kcw.*Tzew + 0.25*Kcw.*Tp)/(1.25*Kcw);
  figure(1)
  plot(Qg,Twew,kolor(mod(i, length(kolor))+1))
  figure(2)
  plot(Qg,Tp,kolor(mod(i, length(kolor))+1))
end

for i = 1:length(Tzew)
  Tp = (Qg + Kcw.*Tzew(i)+5*(Kcp.*Tzew(i)))/(5*Kcp + Kcw);
  Twew = (Qg + Kcw.*Tzew(i) + 0.25*Kcw.*Tp)/(1.25*Kcw);
  figure(3)
  plot(Tzew,Twew,kolor(mod(i, length(kolor))+1))
  figure(4)
  plot(Tzew,Tp,kolor(mod(i, length(kolor))+1))
end

%Zaznaczanie punktow na wykresie
for i = 1:length(Qg)
  Tp = (Qg(i) + Kcw.*Tzew+5*(Kcp.*Tzew))/(5*Kcp + Kcw);
  Twew = (Qg(i) + Kcw.*Tzew + 0.25*Kcw.*Tp)/(1.25*Kcw);
  figure(1)
  scatter(Qg,Twew,kolor(mod(i, length(kolor))+1),'filled')
  figure(2)
  scatter(Qg,Tp,kolor(mod(i, length(kolor))+1),'filled')
end

for i = 1:length(Tzew)
  Tp = (Qg + Kcw.*Tzew(i)+5*(Kcp.*Tzew(i)))/(5*Kcp + Kcw);
  Twew = (Qg + Kcw.*Tzew(i) + 0.25*Kcw.*Tp)/(1.25*Kcw);
  figure(3)
  scatter(Tzew,Twew,kolor(mod(i, length(kolor))+1),'filled')
  figure(4)
  scatter(Tzew,Tp,kolor(mod(i, length(kolor))+1),'filled')
end