%
% Diagramm Carnot-Wirkungsgrad
%
% Menner 2014
%

Titel = "Carnot-Wirkungsgrad";

Text_eta = "Wirkungsgrad η";

T_W = 150:10:1600; % 140 Werte
Text_T_W = "Temperatur TH in Kelvin";
 
T_K1 = 173
T_K2 = 273
T_K3 = 373;
Text_T_K1 = "TK1 = 173 K";
Text_T_K2 = "TK2 = 273 K";
Text_T_K3 = "TK3 = 373 K";
 
% T_W_multi = [T_W; T_W; T_W];
 
% eins_multi = ones;
 
% eta = 1-T_K/T_W
 
eta_1 = 1 - T_K1./T_W;
eta_2 = 1 - T_K2./T_W;
eta_3 = 1 - T_K3./T_W;


eta = [eta_1; eta_2; eta_3];
size(eta)
ylim([0 1])

plot([273 273], ylim, "color", "k", "linewidth", 1.); % 0°C Linie

hold("on");

plot([T_W; T_W; T_W]', eta', "linewidth", 1.)

ylim([0 1]);
xlim([min(T_W), max(T_W)]);
title(Titel, "fontsize", 14);

xlabel(Text_T_W, "fontsize", 12);
ylabel(Text_eta, "fontsize", 12);

legend_set = legend("T = 0 °C", Text_T_K1, Text_T_K2, Text_T_K3, "location", "southeast");

set(legend_set, "fontsize", 12);

% grid();

set (gca(), "linewidth", 1.5);

printf("Weiter mit beliebiger Taste\n");
pause;
 
print("carnot_wirkungsgrad.svg", "-dsvg", "-S1000,800")