N0 = 10000; % Anfangsmenge lambda = 0.1; % Zerfallskonstante h = 0.1; % Schrittweite n = 0:1000; % 1000 Schritte t = h*n; % Zeitwerte N = N0*(1-lambda*h).^n; % Werte fuer N(t) plot(t, N);
plot(t, y(:,1), t, y(:,2)) legend('N_1(t)', 'N_2(t)', 0);
F = cos(t); plot(t, y(:,1), t, y(:,2), t, F) legend('Auslenkung', 'Geschwindigkeit', 'Kraft', 'Location', 'NorthWest');