exercise11.m

function exercise11()
  % Lösung von Aufgabe 11 mit fzero
  mdot = 5.668e-3;        % in kg/s
  Q12dot = 3225;          % in W
  t1 = 70;                % in °C
  M = 44.01e-3;           % in kg/mol, für CO2

  % Daten einlesen
  data = load("cps.dat");
  tD = data(:,1);         % in °C
  cpD = data(:,3)/M;      % in J/(kg K)
  
  % zu lösende Gleichung als Funktion von t2
  f = @(t2) Q12dot - mdot*cp(t1, t2, tD, cpD)*(t2 - t1);

  % Intervall für die Lösung finden und lösen
  epsi = 1e-5;
  cpmin = min(cpD);
  t2 = fzero(f, [t1 + epsi, t1 + Q12dot/(mdot*cpmin)])
end

%---------------------------------------------------------------------
function y = cp(t1, t2, tD, cpD)
  % berechnet die mittlere spezifische Wärmekapazität von t1 bis t2
  % t1        Starttemperatur in °C
  % t2        Endtemperatur in °C
  % tD        Tabellenwerte für t
  % cpD       Tabellenwerte für cp
  % y         Wärmekapazität in J/(kg K)
  
  cp_t1 = interp1(tD, cpD, t1);
  cp_t2 = interp1(tD, cpD, t2);
  y = (t2*cp_t2 - t1*cp_t1)/(t2 - t1);
end