freifall.m

function dy = freifall(t, y)
% rechte Seite der DGL beim freien Fall
g = 9.81;    % in m/s^2

dy1 = y(2);
dy2 = -g;
dy = [dy1 dy2]';