function m = supersize(p)

% Fuer ein gegebenes p in der GDG,
%     m' = (e/k) - (p/k)*m,  m(0) = m0
% berechne den Zeitverlauf m der Masse ueber 30 Tage

% Energiezufluss, pro Tag
    e = 5000;
% Dichte, kcal/kg
    k = 7800;
% Energieabfluss-Parameter, exakter Wert zu bestimmen
%   p = 21.6;

% zu loesende GDG ist m' = dm(t,m)
    dm = @(t,m) (e/k) - (p/k)*m;

% anfaengliche Masse
    m0 = 84;
% Zeitpunkte in der die Loesung ausgegeben wird
    tpte = (0:30)';
% Optionen
    opts = odeset('RelTol',1.0e-6);
% Loesung der GDG
    [t,m] = ode45(dm,tpte,m0,opts);

end
