function ex18()
n = 3000/60;   % Drehzahl in 1/s
p1 = 1e5;
p2 = 18.379e5;
p3 = 43.979e5;
p4 = 2.393e5;
V1 = 1e-3;
V2 = 0.125e-3;
V3 = V2;
V4 = V1;
kappa = 1.40;

% Funktionen 1-2, 3-4
c1 = p1*V1^kappa;
func12 = @(V) c1./V.^kappa;
c3 = p3*V3^kappa;
func34 = @(V) c3./V.^kappa;

% Flächen unter 12 und 34
A12 = integral(func12, V2, V1)
A34 = integral(func34, V3, V4)

% Arbeit und Leistung in J bzw. W
Wa = A34 - A12
Pa = n/2 * Wa

% Funktion für b
al = 30;
m = 1.7;
fx = @(V) (V - V2)/V2;
fy = @(x) 1 + (al*m*x.^m - 1).*exp(-al*x.^m);
fp = @(V) p2*(1 - fy(fx(V))) + fy(fx(V)).*p3.*(V3./V).^kappa;

% analog Flächen etc.
A24 = integral(fp, V2, V1)
Wb = A24 - A12
Pb = n/2 * Wb
