%
% 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")