function dy = testode(t,y)

% Freier Fall mit Stokescher Reibung

% W. Schweizer
% Simulation physikalischer Systeme - Computational Physics mit MATLAB -
% Vers. 1.1
% 02.2016

g =9.81;                 % Erdbechleunigung
k = g/350^2;             % Luftreibung gewichtet mit K?rpergewicht (80 kg)
dy = -g + k.*y.^2;  