laminar.m

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