radiokette.m

function dy = radiokette(t, y)
% rechte Seite der DGL bei radioaktiver Zerfallskette

% Parameter lambda1, lambda2 fest
lambda1 = 0.1;
lambda2 = 0.03;

dy1 = -lambda1*y(1);
dy2 = lambda1*y(1)-lambda2*y(2);
dy = [dy1; dy2];      % Spaltenvektor