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