turbulent.m

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