clear; close all;
% Imie Nazwisko Indeks
% Grupa nr. 1, Wtorek, 7:30
% LAB 02

A= 0.3;
B= 0.7;

qgN= 5000; %[w]   N -> Nominalne
TzewN= -20; %[*C]
TwewN= 20; %[*C]
TpN= 10; %[*C]

Kcw= (A*qgN)/(TwewN-TzewN);
Kcp= (B*qgN)/(TpN-TzewN);
Kcwp2=  (B*qgN)/(TwewN-TpN); % NIE UŻYWANE
Kcwp= 0.25*Kcw;

Tzew= -20 : 5 : 20;
qg= 0 : 500 : 5000;

% Twew = ((Kcwp+Kcp)*(qg+Kcw*Tzew)+Kcwp*Kcp*Tzew) / ((Kcw+Kcwp)*(Kcwp+Kcp)-Kcwp^2);
% Tp = (Kcwp*(qg+Kcw*Tzew)+(Kcw+Kcwp)*Kcp*Tzew) / ((Kcw+Kcwp)*(Kcwp+Kcp)-Kcwp^2);

figure; % creates a new figure window

% w1 = "Kcwp = (B*qgN)/(TwewN-TpN)";
% w2 = 'Kcwp =0.25*Kcw';
% Opis = sprintf("T_{wew}(T_{zew}) \\n %s",w2);
lista=["co "; "kto "; "gdzie"];
opis1 = lista(1,:); opis2 = lista(1,:);
kolor = [[1 0 0]; [0.9 0.1 0.1]; [0.8 0.2 .2]; [0.7 0.3 0.3]; [0.6 0.4 0.4];
       [0.5 0.5 0.5]; [0.4 0.6 0.6]; [0.3 0.7 0.7]; [0.2 0.8 0.8]; [0.1 0.9 0.9]; [0 1 1]]; %kolory RGB

for j = 1 :length(qg)
    subplot(2,2,1) % Twew(Tzew)
    title("T_{wew}(T_{zew})"); xlabel("T_{zew} [-20;20] [°C]"); ylabel("T_{wew}");
    plot(Tzew, ((Kcwp+Kcp)*(qg(j)+Kcw*Tzew)+Kcwp*Kcp*Tzew) / ((Kcw+Kcwp)*(Kcwp+Kcp)-Kcwp^2), 'Color', kolor(j,:));
    hold on
    grid on;
    
    subplot(2,2,2) % Tp(Tzew)
    title("T_p(T_{zew})"); xlabel("T_{zew} [-20;20] [°C]"); ylabel("T_p");
    plot(Tzew, (Kcwp*(qg(j)+Kcw*Tzew)+(Kcw+Kcwp)*Kcp*Tzew) / ((Kcw+Kcwp)*(Kcwp+Kcp)-Kcwp^2), 'Color', kolor(j,:));
    hold on;
    grid on;
    
    opis1(j,:)=sprintf('q_g=%3d', qg(j));
end
legend(opis1, 'Location', 'southeast');
subplot(2,2,1); plot(TzewN, TwewN, "r*");% Twew(Tzew)
subplot(2,2,2); plot(TzewN, TpN, "r*"); % Tp(Tzew)

for j = 1 : length(Tzew)
    subplot(2,2,3) % Twew(qg)
    title("T_{wew}(q_g)"); xlabel("p_g [0;5000] [W]"); ylabel("T_{wew}");
    plot(qg, ((Kcwp+Kcp)*(qg+Kcw*Tzew(j))+Kcwp*Kcp*Tzew(j)) / ((Kcw+Kcwp)*(Kcwp+Kcp)-Kcwp^2), 'Color', kolor(11-j,:));
    hold on;
    grid on;
    
    subplot(2,2,4) % Tp(qg)
    title("T_p(q_g)"); xlabel("p_g [0;5000] [W]"); ylabel("T_p");
    plot(qg, (Kcwp*(qg+Kcw*Tzew(j))+(Kcw+Kcwp)*Kcp*Tzew(j)) / ((Kcw+Kcwp)*(Kcwp+Kcp)-Kcwp^2), 'Color', kolor(11-j,:));
    hold on;
    grid on;
    
    opis2(j,:)=sprintf('Tzew=%3d', Tzew(j));
end
legend(opis2, 'Location', 'southeast');
subplot(2,2,3); plot(qgN, TwewN, "b*"); % Twew(qg)
subplot(2,2,4); plot(qgN, TpN, "r*"); % Tp(qg)