>> 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');