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