hold on;
grid on;
color = 'rgybc';
title('Zaleznosc temperatury wewnatrz domu od temepratury zewnetrznej');
xlabel('Tzew [*C]');
ylabel('Twew [*C]');
W=0.25;                                               %wspolczynnik
Twew=20;                                              %temperatura wewnątrz
Tzew=-20;                                             %temperatura na zewnatrz
Tp=10;                                                %temperatura poddasza
Qg=1000;                                              %moc grzejnika elektrycznego
Kcw=Qg/((Twew-Tzew)+W*(Twew-Tp));                     %ciepło tracone przez ściany
Kcwp=W*Kcw;                                           %cieplo tracone przez sufit
Kcp=Kcwp*(Twew-Tp)/(Tp-Tzew);                         %ciepło tracone przez dach
Tzew=[-20:1:20];
Twew=(Qg+Tzew*(Kcw+( (Kcwp*Kcp)/(Kcwp+Kcp) ) ) )/( Kcw+Kcwp-((Kcwp*Kcwp)/(Kcwp+Kcp)));
%plot(Twew,Tzew); % pokaz charakterystyke
Qg=[500,1000,1500,2000]; % różne wartości na mocy na grzejniku

for i=1:1:length(Qg); % dla różnych mocy grzejnika
        Twew=(Qg(i)+Tzew*(Kcw+( (Kcwp*Kcp)/(Kcwp+Kcp) ) ) )/( Kcw+Kcwp-((Kcwp*Kcwp)/(Kcwp+Kcp)));

        plot(Tzew,Twew,color(mod(i,length(color))+1),'LineWidth',4);

end