function dy = turbulent(t, y) % rechte Seite der DGL beim freien Fall mit turbulenter Reibung g = 9.81; % in m/s^2 m = 75; % in kg b2 = 0.1; % in kg/m dy1 = y(2); dy2 = -g - b2/m*y(2).^2.*sign(y(2)); dy = [dy1 dy2]';