function dy = laminar(t, y) % rechte Seite der DGL beim freien Fall mit laminarer Reibung g = 9.81; % in m/s^2 m = 75; % in kg b1 = 5; % in kg/s dy1 = y(2); dy2 = -g - b1/m*y(2); dy = [dy1 dy2]';